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ABSTRACT 


Two  particular  source  regions,  Lake  Superior  and  the 
Nevada  Test  Site  are  documented  in  detail  using  available 
geological  and  geophysical  observations.  Each  is  discussed 
individually  in  relation  to  its  setting  in  the  prevailing 
regional  crust  and  upper  mantle  structure  of  the  mid-continent 
and  the  Basin  and  Range  respectively.  Comparison  of  the  two 
source  areas  indicates  that  they  are  markedly  different  from 
one  another  not  only  at  shallow  depths  but  throughout  the 
upper  mantle  as  well.  Location  accuracies  for  events  at  both 
sites  are  quite  comparable,  however,  a  large  bibliography 
documents  the  extensive  literature  relavant  to  Lake  Superior 
and  the  Nevada  Test  Sites. 
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INTRODUCTION 


Between  July  8  and  July  30,  1963  a  series  of  78  chemical 
explosions  was  detonated  along  three  shot  lines  in  Lake  Superior 
(Figure  1)  as  part  of  a  large  cooperative  undertaking  designed 
to  study  intensively  the  crust  and  upper  mantle  beneath  a 
limited  geographical  region  by  explosion  seismology.  Because 
crustal  structure  deduced  from  explosion  seismology  has  long 
been  confused  by  the  subjective  judgements  of  the  individual 
interpreting  the  data,  a  number  of  groups  from  institutions  in 
the  United  States  and  Canada,  collectively  referred  to  as  the 
North  American  Seismic  Group  (for  a  list  of  participants,  see 
Table  1)  agreed  to  record  these  shots  in  order  to  secure  a 
number  of  independent  interpretations.  In  addition,  the 
participants  agreed  that  the  first-arrival  data  be  made 
available  in  one  place  (Steinhart,  1964)  and  to  provide  access 
to  the  seismograms  obtained. 

The  Lake  Superior  region  was  chosen  because  previous  in¬ 
vestigators  (Steinhart  and  Meyer,  1961)  had  suggested  that  the 
shield  areas  of  the  world  are  especially  uniform  and  structurally 
simple;  other  results  from  several  previous  studies  of  crustal 
structure  in  the  area  of  the  southern  Canadian  Shield  surround¬ 
ing  Lake  Superior  (Slichter,  1951;  Tatel  and  Tuve,  1952; 

Hodgson,  1953  a,b)  lent  support  to  this  contention.  However, 
the  expected  simple  structure  was  not  found.  Rather,  the 
crust  under  Lake  Superior  is  probably  as  structurally  complex 
as  anywhere  in  North  America.  To  date,  the  pooled  first  arrival 
data  obtained  during  the  1963  Lake  Superior  Experiment  have 
been  interpreted  by  three  groups:  Smith  et  al,  (1966a, b)  and 
Berry  and  West  (1966a, b)  used  the  time  term  approach  (Scheidigger 
and  Willmore,  1957)  whereas  OfBrien  (1968)  used  the  method  of 

matching  structure  to  reduce  the  data  and  deduce  a  velocity- 
depth  model.  As  is  to  be  expected,  the  three  models  differ 
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significantly,  although  all  suggest  that  the  seismic  velocities 
in  the  crust  are  extraordinarily  high  (P-wave  velocities  great¬ 
er  than  6.6  km/sec  at  depths  of  only  10  km)  and  that  in  an 
overall  sense  the  structure  beneath  the  Lake  is  synclinal.  In 
the  discussions  which  follow,  we  examine  the  validity  of  the 
assumptions  underlying  these  three  interpretations  and  attempt 
to  determine  what  effect,  if  any,  they  might  have  on  the  de¬ 
rived  models. 

The  one  ton  charges  detonated  during  the  1963  Lake  Superior 
Seismic  Experiment  were  recorded  at  surprisingly  large  distances 
by  stations  in  the  LRSM  and  VELA  observatory  network  (Steinhart, 
1964;  Mansfield  and  Evernden,  1966a, b;  Bancroft,  1966).  In 
October  of  1964  the  U.S.  Geological  Survey  detonated  a  series 
of  ten  shots  ranging  in  size  from  2  to  20  thousand  lbs  of  TNT  in  an 
experiment  designed  to  investigate  the  unusually  low  attenua¬ 
tion  of  P-waves  in  the  Interior  Plains  of  the  United  States 
that  was  observed  during  the  1963  experiment  and  to  determine 
whether  an  upper-mai  cle  "shadow  zone"  for  P-waves  exists 
beneath  the  plains  (Roller  and  Jackson,  1966a, b).  The  U.  S. 
Geological  Survey  recorded  these  ten  shots  at  35  stations  along 
a  profile  extending  from  the  Bayfield  Peninsula  to  Denver* 

Colorado.  Roller  and  Jackson  (1966  a,b)  analyzed  the  travel-time 
and  amplitude  data  recorded  along  this  profile  and  found  that 
the  P-wave  velocity  increased  with  depth  to  100  km  at  which 
depth  they  found  a  high-velocity  (8.5  km/sec)  layer.  They  also 
noted  that  the  amplitude  of  Pn  attenuated  as  1/A,  in  contrast 

3 

with  the  attenuation  of  1/A  commonly  observed  in  the  western 
United  States. 

The  successful  observation  of  the  1963  and  1964  shots 
at  large  distances,  plus  the  -nusually  low  attenuation  and  high 
Pn  velocities  compared  with  the  western  United  States  that  were 

-2- 


t 


observed  during  the  Lake  Superior  and  other  seismic  experi¬ 
ments  in  the  eastern  United  States  led  to  a  third  seismic  ex¬ 
periment  at  Lake  Superior  in  1966.  Named  Project  Early  Rise, 
this  experiment  consisted  of  38  shots  each  of  10,680  lbs.  of 
Dupont  Nitramon  WW(EL)  which  were  detonated  at  a  common  shot 
point  off  the  Keweenaw  Peninsula  during  a  twenty-one  day  period 
in  July  1966  and  recorded  along  extensive  profiles  which  ra¬ 
diated  outwards  from  the  shores  of  the  Lake  (Figure  2).  The 
recording  profiles  were  located  so  that  the  data  would  yield 
information  on  the  structure  of  the  crust  and  upper  mantle 
beneath  the  central  portion  of  the  North  American  continent. 

Ten  members  of  the  North  American  Seismic  Group  (Table  2) 
occupied  the  stations  on  these  profiles.  Warren  et  al,  (1967) 
have  published  a  compendium  of  first  arrival  data  obtained 
during  Project  Early  Rise. 

Lewis  and  Meyer  (1968),  Green  and  Hales  (1968),  and  Mereu 
and  Hunter  (1969)  have  interpreted  some  of  the  data  collected 
during  Project  Early  Rise  and  their  results  indicate  that  the 
structure  of  the  upper  mantle  beneath  the  Canadian  Shield  is 
quite  different  from  that  beneath  the  western  United  States  as 
found  by  Niazi  and  Anderson  (1965),  Johnson  (1967)  and  Archambeau 
et  al,  (1969).  Thus  in  comparing  Lake  Superior  as  a  seismic  • 
source  region  with  the  Nevada  Test  Site  we  must  examine  the 
regional  setting  of  the  two  areas.  To  accomplish  this,  we  have 
attempted  to  summarize  what  is  currently  known  about  the  physical 
properties  of  the  crust  and  upper  mantle  beneath  the  two  areas 
and  derived  geophysical  models  for  the  crust  beneath  the  Early 
Rise  shot  point  and  Yucca  Flat,  Nevada  Test  Site. 


In  constructing  possible  crustal  models  for  the  Early  Rise 
shot  point  and  the  Yucca  Flat  area  of  the  Nevada  Test  Site  we 


they  are  known.  We  should  point  out,  however,  that  we  have 
been  unable  to  apply  any  tests  for  uniqueness  to  these  models. 
Nevertheless,  in  view  of  the  constraints  placed  on  the  possible 
choices  by  the  varied  data  which  we  consider,  we  believe  that 
the  models  are  sufficiently  representative  of  the  crust  and 
upper  mantle  in  the  two  areas  to  serve  as  a  basis  for  compan- 
son  of  the  two  source  regions. 
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PHYSICAL  SETTING  OF  LAKE  SUPERIOR  EXPERIMENTS 


Geologic  Environment 

The  Relationship  *>f  Lake  Superior  to  the  Canadian  Shield 


The  Canadian  Shield  forms  a  large  stable  area  of  more 
than  two  million  square  miles  in  the  northeastern  part  of  North 
America  (Figure  3).  A  shield  is  typically  an  area  of  Precambrian 
igneous  and  metamorphic  rocks  which  have  been  eroded  to  a  con¬ 
vex  surface  with  little  relief  and  overlain  at  its  margins  by 
much  younger  gently  dipping  sedimentary  rocks  which  increase  in 
thickness  away  from  the  shield.  Where  these  sediments  attain 
considerable  thickness,  the  rocks  of  the  shield  form  the  basement. 
These  latter  areas  which  have  remained  relatively  unaffected  by 
later  orogenies  are  generally  referred  to  as  stable  continental 
interiors,  or  in  the  case  of  the  United  States,  the  Interior 
Lowlands  (King,  1959,  p.  23). 

The  Canadian  shield  departs  from  the  ideal  case  in  several 
important  ways  as  a  result  of  Phanerozoic  deformations  (King 
1959,  p.  11).  Beginning  in  Cretaceous  time  (King,  ibid.)  the 
northeastern  part  of  the  shield  underwent  severe  deformation 
which,  continuing  through  the  Tertiary,  resulted  in  the  dissec¬ 
tion  and  uplift  of  the  present  day  northeastern  margin  of  the 
shield  to  form  the  Torngat  Mountains  of  Labrador  and  the 
mountainous  peaks  of  Baffin  Land.  The  relationship  of  Greenland 
to  the  Canadian  Shield  is  problematical.  However  if  Greenland 
once  constituted  part  of  the  shield,  it  is  likely  that  it  was 
separated  from  the  major  part  of  the  shield  at  the  same  time 
that  the  '’mountains”  formed.  During  the  Pleistocene  and 
Recent,  the  Canadian  Shield  was  subjected  to  glacial  load¬ 
ing  which  resulted  in  its  being  depressed  beneath  the  center 
of  the  ice  caps.  Of  the  three  ice  caps  which  originally 
covered  the  shield,  only  the  ice  cap  on  Greenland  remains. 
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Glacial  erosion  considerably  altered  the  topography  of 
the  shield,  denuding  it  of  much  of  its  sedimentary  cover  and 
forming  the  many  small  lakes  of  the  shield's  interior  and  the 
much  larger  lakes  (the  Great  Lakes)  along  its  margins.  The 
Canadian  Shield  conforms  more  closely  to  the  definition  of  a 
shield  area  along  its  southern,  south-east  and  south-west 
borders  (King,  1959,  p.  10).  Here  the  shield  passes  beneath 
the  sedimentary  rocks  of  the  Interior  Lowlands  whose  scarps 
face  the  interior  of  the  shield. 

Seismological  studies  of  the  crust  and  upper  mantle  be¬ 
neath  the  Canadian  Shield  reflect  this  long  period  of  relative 
stability.  Brune  and  Dorman  (1963)  studied  phase  velocities 
for  fundamental  mode  Love  and  Rayleigh  waves  along  extended 
paths  across  the  Canadian  Shield  supplemented  by  group  veloc¬ 
ity  data  and  travel  times  for  body  waves.  Figure  4  shows 
the  model  CANSD  which  best"  fits  their  data.  Brune  and  Dorman 
found  that  the  crust  was  relatively  uniform  over  the  shield 
and  could  be  approximated  by  a  three- layer  crust  some  35  km 
thick.  Shear  wave  velocities  in  the  upper  mantle  indicate  a 
low-velocity  zone  about  200  km  thick  beginning  at  a  depth 
of  115  km  where  the  shear  wave  velocity  decreases  by  0.2  km/ 
sec.  CANSD  shows  the  P-wave  velocity  slowly  increasing  with 
depth  from  8.1  km/sec  at  the  Moho  to  8.3  km/sec  at  315  km 
where  there  is  a  sharp  increase  to  8.7  km/sec. 

Recently  Barr  (1967)  has  published  the  results  of  a 
study  of  P-wave  data  from  chemical  explosions  in  Hudson  Bay 
and  Lake  Superior.  Figure  4  also  shows  his  proposed  model  for 
P-wave  velocities  in  the  upper  mantle  beneath  the  Canadian 
Shield.  The  discontinuities  in  P-wave  velocity  are  determined 
from  the  positions  of  cusps  in  the  observed  travel-time  data. 
The  depths  of  the  discontinuities  in  P-wave  velocity  are  in 
close  agreement  with  the  depths  of  the  discontinuities  in 
S-wave  velocity  found  by  Brune  and  Dorman  (1963).  Barr  (1^67) 
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finds  no  evidence  for  a  P-wave  low  velocity  zone  and  concludes 
that  the  P-wave  velocity  structure  in  the  upper  mantle  is 
comparatively  uniform  across  the  Canadian  Shield. 

Heat  flow  measurements  in  continental  shield  areas,  in¬ 
cluding  the  Canadian  Shield,  reflect  their  geologic  stability. 

In  a  review  of  heat  flow  data  Lee  and  Uveda  (1965)  concluded 

2 

that  the  global  mean  heat  flow  is  1.6  +_  10%  pcal/cm-sec  at 
95%  confidence  leve.1 .  They  found  no  significant  difference 
between  the  mean  oceanic  and  mean  continental  heat  flow  but 
concluded  that  in  general,  regional  heat  flow  variations  great- 

r\ 

er  than  0.2  ycal/cm-sec  are  significant.  For  the  Canadian 
Shield,  Lee  and  Uyeda  found  a  mean  heat  flow  of  0.9  +  0.1  P 
cal/cm-sec  for  145  measurements  made  in  Lake  Superior  as  re¬ 
ported  by  Steinhart  et  alj  (1968).  The  low  heat  flow  in  the 
Canadian  Shield  is  not  due  solely  to  the  effects  of  Pleistocene 
glaciation  since  the  mean  for  seven  determinations  of  heat  flow 
in  the  Australian  Shield,  which  was  not  glaciated,  is  1.0  +_ 

0.1  (5)  pcal/cm-sec  (Lee  and  Uyeda,  1965).  Rather,  low  heat 
flow  values,  0.9  +  0.2  pcal/cm-sec,  are  characteristic  of 
shield  areas  throughout  the  world.  Considering  that  shield 
areas  have  an  abundance  of  radiogenic  rocks  in  the  crust 
(granites,  gneisses,  etc.)  which  contribute  to  the  measured 
heat  flow,  the  mantle  heat  flow  under  the  shields  must  be 
very  low  indeed  compared  to  oceanic  and  tectonically  active 
areas . 

Although  the  Canadian  Shield  has  behaved  as  a  single 
entity  since  Cambrian  times  it  is  composed  of  a  series  of 
provinces  with  different  ages,  structural  histories  and  lithol¬ 
ogies.  The  shield  is  presently  divided  into  seven  geologic 
provinces  (Figure  3).  Radiometric  dating  (Stockwell  1962, 

1964)  has  shown  that  the  rocks  of  the  Superior  Province  were 
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involved  in  the  Kenoran  orogeny  2.5  +_  0.1  b.y.  ago,  those 
of  the  Churchill  Province  in  the  Hudsonian  orogeny  1.7  +_ 

0.1  b.y.  ago,  and  those  of  the  Grenville  Province  in  the 
Grenville  orogeny  0.9  +_  0.06  b.y.  ago. 

Data  obtained  in  recent  seismic  refraction  surveys 
suggests  that  the  gross  crustal  structure  and  thickness  also 
varies  among  the  three  major  provinces.  Mereu  and  Hunter 
(1969)  in  an  analysis  of  data  obtained  during  Project  Early 
Rise  (Warren  et  al  ,  1967)  conclude  that  the  crust  is  30  - 

35  km  thick  under  the  Superior  Province  to  the  west  of  Hudson 
Bay  whereas  under  the  Churchill  Province  to  the  north  it  is 
greater  than  40  km  thick.  The  indicated  increase  in  crustal 
thickness  takes  place  across  a  region  about  200  km  wide  to 
the  north  of  the  Nelson  River  and  correlates  with  a  large  linear 
gravity  anomaly  (Innes,  1960).  Mereu  and  Hunter  (1967)  suggest 
that  a  model  consisting  of  a  single  6.3  km/sec  layer  for  the 
crust,  a  P  velocity  of  8.05  km/sec  and  a  8.43  km/sec  layer 
at  a  depth  of  84  O  km  is  required  to  fit  their  data. 

The  crustal  thickness  of  40  km  found  by  Mereu  and  Hunter 
(1969)  to  the  west  of  Churchill  agrees  with  the  earlier  re¬ 
sults  of  Hunter  and  Mereu  (1967),  Ruffman  and  Keen  (1967), 
and  Hobson  et  al,  (1967).  These  authors  interpreted  the 
travel-times  recorded  during  the  1966  Hudson  Bay  Experiment 
(Figure  5)  and  though  specific  details  in  interpretation  differ 
among  the  three  papers,  all  conclude  that  the  crust  is  approxi¬ 
mately  40  km  thick  at  Churchill  (Figure  5),  decreases  to  25  - 
30  km  150  km  west  of  the  Ottawa  Islands  and  then  increases 
rapidly  to  40  km  again  under  Povungnituk.  The  data  from  the 
N-S  profile  suggest  that  the  crust  is  25  -  30  km  thick  under 
Chesterfield  Inlet  and  increases  in  thickness  to  35  -  40km  in 
the  center  of  the  bay.  Th^  rapid  12  km  increase  in  crustal 
thickness  to  the  west  and  north  of  Ottawa  Island  corresponds 
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to  the  projection  of  the  Superior-Churchill  province  boundary 
eastwards  beneath  the  Hudson  Bay.  The  data  also  indicate  a 
single  layered  crust  with  a  velocity  of  6.3  km/sec  covered  by 
a  thin  veneer  of  Phanerozoic  sediments  which  occupy  the  sedi¬ 
mentary  basin  which  forms  the  Hudson  Bay.  This  sedimentary 
basin  is  not  reflected  in  the  Moho  topography  which  consists 
of  a  N-S  trending  trough  superimposed  upon  an  E-W  oriented 

saddle.  All  three  papers  report  a  P  velocity  of  8.2  to  8.3 

n 

km/sec  beneath  the  Hudson  Bay. 

Thus,  geophysical  studies  to  date  indicate  that  the  upper 
mantle  beneath  the  Canadian  Shield  is  laterally  relatively  uni¬ 
form.  The  crustal  structure  however  is  not  lateraly  uniform. 
The  various  provinces  recognized  by  geologists  appear  to  have 
differences  in  gross  crustal  structure  which  extend  to  the 
base  of  the  crust.  The  greater  part  of  Lake  Superior  is  sit¬ 
uated  in  one  of  these  provinces,  the  Southern  Province,  along 
the  southern  margin  of  the  shield.  There  is  every  indication 
that  crustal  structure  beneath  the  Lake  is  as  complicated  as 
anywhere  else  in  the  Canadian  Shield. 

Geology  of  the  Lake  Superior  Basin 

The  geology  of  the  Lake  Superior  region  has  been  under 
investigation  for  more  than  a  century.  Excellent  accounts  of 
the  earlier  work  can  be  found  in  Lane  (1911),  Van  Hise  and 
Leith  (1911),  and  Butler  and  Burbank  (1929).  More  up-to-date 
accounts  which  take  into  consideration  recent  geophysical  data 
can  be  found  in  White  (1966ar  1966b)  and  Halls  (1966).  These 
geologic  data  and  their  interpretation  provide  the  background 
necessary  to  interpret  the  results  obtained  during  the  1963-4 
Lake  Superior  Experiments  and  Project  Early  Rise  of  1966. 
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Rocks  outcropping  in  the  immediate  vicinity  of  Lake 
Superior  can  be  divided  into  four  groups  in  order  of  in¬ 
creasing  age: 

1.  Recent  glacial  deposits 

2.  Paleozoic  sediments 

3.  Keweenawan  rocks 

4.  pre-Keweenawan  rocks 

Figure  6  shows  the  areal  distribution  of  pre-glacial  rocks  in 
the  region.  Table  3  summarizes  the  pertinent  stratigraphic 
data.  Cores  taken  in  the  Lake  itself  (Zumberge  and  Cast,  1961; 
Zumberge,  1962)  suggest  that  the  Cambrian  rocks  which  outcrop 
along  the  southeast  shore  of  Lake  Superior  extend  beneath  parts 
of  the  lake. 

The  oldest  rocks  outcropping  in  the  Lake  Superior  region 
are  the  Archean  rocks  of  the  Ontarian  and  Timinskamian  (Halls, 
1966)  or  Timiskaming  (Stockwell,  1964)  and  thus  are  older  than 
2.4  b.y.  These  rocks  then  formed  the  basement  of  an  E-W  trend¬ 
ing  geosyncline  which  occupied  much  of  the  western  part  of  the 
present  day  lake  and  probably  extended  eastwards  along  the 
northern  shores  of  Lake  Huron,  The  lack  of  Animikean  or  Huron- 
ian  rocks,  the  structural  trend,  and  the  Kenoran  dates  of  the 
rocks  exposed  along  the  eastern  shore  of  Lake  Superior  suggest 
that  the  Penokean  orogeny  did  not  extend  into  western  Ontario, 
Deposition  of  Animikean  rocks  and  their  equivalents  in  this 
trough  terminated  with  the  Penokean  orogeny  1,8  -  1.6  b.y,  ago 
which  is  equated  by  Stockwell  to  the  Hudsonian  orogeny  in  the 
Churchill  province  to  the  north.  The  Animikean  rocks  along 
the  northern  shore  of  Lake  Superior  were  only  slightly  defor¬ 
med  by  the  Penokean  orogeny,  the  deformation  being  more  inten¬ 
sive  in  central  Minnesota,  northern  Michigan,  and  Wisconsin. 

After  a  short  period  of  erosion  a  large  structural  basin 


deve^ped  in  which  lower  Keweenawan  sediments  accumulated. 

Thef  i  are  overlain  by  a  very  thick  series  of  basic  lava  1  low  , 
the  Middle  Keweenawan  lavas.  The  sources  for  the  lavas  have 
not  been  recognized.  Subsidence  of  the  Lake  ^jperior  basin 
kept  pace  with  the  eruption  of  the  lavas  (White,  1966b),  an* 
continued  after  the  volcanism  had  stopped  and  the  deppsition 
of  the  Upper  Keweenawan  s'diments  began.  The  Keweenawan  lavas 
are  intruded  by  the  Duluth  Gabbro  which  has  a  radiometric  age 
of  1.0  b.y .  The  lavas  are  thus  younger  than  1.6  b.y.  and 
older  than  1.0  b.y.  The  age  of  the  Upper  Keweenawan  sedi¬ 
ments  is  uncertain.  They  are  younger  than  the  lavas,  but 
sedimentation  in  the  Lake  Superior  basin  continued  into  the 
Cambrian  (Hamblin,  1961)  with  only  minor  interruptions, 
which  suggests  that  they  were  deposited  between  1.  b.y.  and 

0,6  b.y.  ago. 

Lake  Superior  itself  is  a  large  (o3,  *0  km  )  temper ate 

fresh-water  lake  with  a  mean  depth  of  ISO  m  and  a  maximum 
depth  of  400  m.  (Hart  and  Tilton,  1966),  which  occupies  the 
greater  part  of  the  Keweenawan  basin.  Figure  7  shows  the 
configuration  of  the  bottom  topography.  In  the  west  the 
lake  bottom  is  strongly  asymetric  being  deepest  along  the 
Minnesota  shore,  whereas  in  the  eastern  sector  the  lake 
bottom  is  made  up  of  a  series  of  ridges  and  troughs  orientated 
approximately  N-S.  These  latter  features  are  probably  the  re¬ 
sult  of  glacial  scouring  (Weld  and  Ostenso,  1966). 

At  the  present  day  the  Keweenawan  rocks  form  a  large  syn¬ 
cline  first  recognized  by  Irving  (1883),  most  of  which  is 
covered  by  the  lake  (Figure  6).  In  the  western  half  of 
the  lake  the  inferred  synclinal  axis  trends  HE.  Off  the 
Keweenaw  Peninsula  it  trends  E-W,  and  in  the  eastern  half  of 
the  lake  it  trends  HW-SE.  There  is  strong  evidence  from 
gravity,  aeromagnetic,  borehole  and  seismic  data  (Lyons,  -  ’* 
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Theil,  1956;  Craddock  et  al,  1963;  Zietz,  1965;  Cohen,  1966; 
Cohen  and  Meyer,  1966),  that  the  Keweenawan  ro cV"L  continue 
southwards  from  the  southwest  end  of  Lake  Superior  in  a  narrow 
trough  as  far  as  Kansas  and  that  they  are  closely  associated 
with  the  mid-continent  gravity  high  (Woollard  and  Joesting, 
196M),  Gravity  dcta  obtained  by  Bacon  (1957)  tentatively  sug¬ 
gest  that  the  Keweenawan  rocks  persist  southward  from  the 
eastern  end  of  the  Lake  also. 

Although  the  gross  structure  of  the  Lake  Superior  region 
has  been  known  fcr  some  time,  the  details  of  the  stricture  are 
still  poorly  known.  White  (1966b)  has  interpreted  the  struc¬ 
ture  between  Isle  Roy ale  and  the  Keweenaw  Peninsula  from  geo¬ 
logic  data.  His  interpretation  is  shown  in  Figure  8.  White 
gives  the  permissible  range  of  thickness  for  the  sedimentary 
rocks  as  3  to  7  km.  The  inferred  thickness  of  the  sedimentary 
section  is  dependent  upon  how  the  dip  of  the  upper  surface  of 
the  Keweenawan  lavas  is  extrapolated;  the  most  probable  thick¬ 
ness  is  5km.  The  total  thickness  of  the  lavas  is  unknown* 
but  White  (1966b)  estimates  them  to  be  6.5  km  thick  cn  Keweenaw 
Point  and  m.s  km  thick  on  Isle  Royale  based  on  Lane's  (1911) 
stratigraphic  correlation  between  the  two  outcrops,  although 
White  points  out  that  they  may  be  as  thick  as  9.5km, 

Geophysical  Evidence  on  Crustal  Structure  Under  Lake  Superior 

Seismic  refraction,  gravity  and  aeromagnetic  surveys 
have  been  conducted  in  the  vicinity  of  Lake  Superior  in  re¬ 
cent  years.  Since  1963  extensive  seismic  refraction,  under¬ 
water  gravity,  magnetic  and  heat  flow  surveys  have  been  made 
in  the  lake  itself.  The  results  of  much  of  this  work  are 
published  in  Monograph  10  of  the  American  Geophysical  Union 
(Steinhart  and  Smith,  1966)  and  the  Carnegie  Institution  of 
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Washington  Year  Books  for  1963  through  1967.  Because  the 
seismic  data  are  of  primary  concern  in  this  report,  we  shall 
examine  them  first  and  attempt  to  discriminate  among  those 
crustal  models  which  fit  the  seismic  data  by  using  the  remain¬ 
ing  geophysical  and  geological  data. 

Seismic  Measurements 

Ref notion  Studies.  -  The  1963  Lake  Superior  Seismic 
Experiment  consisted  primarily  of  a  series  of  one  ton  shots 
firei  along  a  line  from  Duluth  to  Otter  Cove  (Tigure  1),  which 
were  recorded  by  a  number  of  stations  (Figure  9)  situated  at 
either  end  of  thi"*  profile.  In  addition,  data  were  obtained 
from  two  secondary  profiles,  the  Bayfield  Peninsula-Sibley 
Peninsula  profile  also  referred  to  as  the  Gambler  Line  in 
some  publications,  and  the  Keweenaw  Peninsula-Hichipicoten 
profile  or  South  Line.  Steinhart  (1964)  has  publishea  shot 
locations  and  origin  times,  station  locations,  and  first  arrival 
data  from  this  experiment. 

When  the  first  arrival  data  for  stations  at  either  end 
of  the  profile  are  plotted  as  reversed  travel- time  distance 
plots  (for  example,  see  Smith  et  al,  1966a;  Berry  and  West, 
1966a)  it  is  clear  that  a  simple  dipping-layer  model  is  in¬ 
adequate  to  explain  the  observations.  In  an  attempt  to  over¬ 
come  these  difficulties,  Berry  and  West  (1966a, b)  and  Smith 
et  al,  (1966a, b)  have  interpreted  the  first  arrival  data  using 
the  time-term  method,  originally  described  by  Scheidegger  and 
Willmore  (1957). 

The  velocity-depth  crustal  models  obtained  by  Berry  and 
West  and  r.ralth  et  al,  are  shown  in  figure  10,  Both  models  are 
for  the  main  profile  orientated  approximate ly  E-W  (figure  1). 
Each  model  has  a  two  layer  crust,  although  the  thickness  and 
velocity  of  the  two  layers  varies  between  the  two  models.  The 
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models  a'*o  in  agreement  to  the  extent  that  they  show  a  con¬ 
siderable  thickening  (some  15  -  20  km)  of  the  crust  to  the 
east  of  the  Apostle  Islands.  The  Berry  and  West  model  (Berry 
and  West,  1966a,b),  has  an  upper  layer  about  3  -  5  km  thick 
under  the  shores  of  the  lake  with  an  assumed  velocity  of  5.5 
km/sec  which  attains  a  maximum  thickness  of  12  km  west  of  the 
Apostle  Islands.  The  lower  crustal  layer  has  an  assumed 
velocity  of  6.6  km/sec.  The  depth  to  the  Moho  is  30  km  at 
the  western  extremity  of  the  profile  decreasing  to  26  km  off 
the  western  shore  of  the  lake  and  increasing  to  a  maximum  of 
57  km  approximately  70  km  west  of  Otter  Cove*,  to  the  east  of 
this  point  the  crustal  thickness  decreases  abruptly  to  u8  km. 
Smith  et  al,  (1966a)  interpret  this  sharp  change  as  a  fault. 

The  upper  layer  of  Smith  et  alj's  model  (Figure  10)  shows 
considerable  lateral  variation  in  both  velocity  and  thickness. 
The  thickness  is  roughly  constant  at  5  km  under  the  lake,  but 
increases  rapidly  to  more  than  10  km  to  the  east  of  Otter  Cove. 
The  velocity  in  this  layer  (assumed)  varies  laterally  from  m.8 
km/sec  to  6.2  km/sec.  The  lower  crustal  layer  everywhere  has 
a  velocity  of  6.8  km/sec.  Depths  to  the  Moho  are  similar  to 
Berry  and  Wes^s  model.  They  range  from  35  Van  under  station 
Zipgun,  to  28  km  under  the  western  shore  of  the  lake,  in¬ 
creasing  rapidly  to  63  km  some  70  km  west  of  Otter  Cove.  Again 
there  is  a  sharp  decrease  in  depth  at  this  point  of  around  10  km 

and  the  depth  to  the  Hoho  continues  to  decrease  eastwards  to 

• 

30  km  under  Foylet.  The  Pn  velocity  of  8.07  km/sec  closely 
agrees  with  the  8.1  km/sec  found  by  Berry  and  West. 

The  Lake  Superior  Seismic  Experiment,  as  conducted,  does 
not  meet  the  theoretical  requirements  of  the  time-term  method 
rigorously.  In  particular,  the  Pn  time-termn  for  shots  close 
to  either  end  of  the  lake  were  obtained  mainly  from  stations 
located  on  only  one  side  of  the  shot.  The  structure  of  the 


crust-mantle  interface  from  the  Pn  time-term  analysis  is 
thus  considerably  exaggerated  and  displaced  westward  by 
approximately  70  km  (O'Brien,  1968)  since  there  was  a  pre¬ 
ponderance  of  recording  stations  along  the  western  shore  of 
the  lake.  O'Brien  (1968)  has  reinterpreted  the  data  from  the 
1963  Lake  Superior  Experiment  using  the  method  of  matched  pro¬ 
files.  (O'Brien's  model,  again  for  the  main  E-W  profile,  is 
also  shown  in  Figure  10). 

The  upper  layer  in  O'Brien  's  model  is  6  -  8  km  thick  and 

represents  the  near  surface  rocks.  The  velocity  in  this  layer 
varies  laterally  from  4.7  km/sec  to  4.9  km/sec,  and  was  de¬ 
termined  rather  than  assumed.  This  surface  layer  is  under¬ 
lain  by  a  single  layer,  the  upper-surface  velocity  of  which 
varies  laterally  from  6.6  km/sec  to  6.7  km/sec  and  in  which 
the  velocity  increases  slightly  with  depth  to  6  75  km/sec  at 
the  base  of  the  crust.  The  depth  to  the  Moho  ranges  from  45  km 
under  the  western  shore  of  the  lake  to  35  km  off  the  Bayfield 
Peninsula,  increases  to  a  maximum  of  54  km  approximately  130  km 
west  of  Otter  Cove  and  then  decreases  to  42  km  under  the  eastern 
shore  of  the  lake.  The  weighted  mean  crustal  velocity  for  this 
model  is  6.75  km/sec  which  agrees  well  with  the  6.85  km/sec 
which  O'Brien  obtained  from  t  -x  analyses  of  wide-angle  re¬ 
flections  from  the  Moho,  the  6.8  km/sec  used  by  Smith  et  al, 
to  convert  the  Pn  time-terms  to  depths,  and  the  6.9  km/sec 
layer  at  a  depth  of  8  km  under  the  western  half  of  Lake 
Superior  reported  by  Anzoleaga  et  al,  (1969). 

O'Brien  also  discusses  a  three  layered  crustal  model,  the 
thira  layer  of  which  has  a  velocity  of  around  7.0  km/sec.  This 
model  (Figure  11)  is  based  on  intra-crustal  reflections  re¬ 
corded  at  Otter  Cove,  Tom,  and  Woodsman.  The  depth  to  the  re¬ 
flector  io  15  km  near  the  eastern  shore  and  25  km  in  the  east- 
central  sector  of  the  lake.  In  addition,  O'Brien  attributes 
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late  arrivals  recorded  at  stations  X-ray  and  Tyro  in  the 
western  sector  of  the  lake  to  reflections.  However,  these 
arrivals  are  more  scattered  and  could  not  be  correlated  over 
any  distance.  The  records  at  the  stations  Yankee  and  Victor 
did  not  snow  distinct  reflections,  but  the  amplitudes  of  tns 
first  arrival  were  10  times  those  recorded  at  .stations  Tyro, 
Woody,  and  X-ray.  Similar  large  amplitude >  were  also  found 
by  Roller  and  Jackson  (1966a)  to  the  southwest  of  the  lake. 

0  Brien  suggests  that  these  are  due  to  reinforcement  of  the 
first  arrival  by  the  reflection  on  account  of  the  small 

travel-time  difference  between  these  two  phases  at  this  dis¬ 
tance  . 


The  small  amplitude  of  the  reflections  in  the  eastern 
section  of  the  lake  suggests  that  they  may  be  due  to  a  second- 
order  discontinuity  in  velocity.  Until  additional  evidence 
is  forthcoming,  we  prefer  the  two  layer  model.  The  wide 
angle  reflection  data  does  indicare,  according  to  01 Brien 
(1968),  that  the  P-wave  velocity  increases  with  depth  be¬ 
neath  the  lake  and  suggests  that  there  are  no  low  velocity 
layers  present  within  the  crust. 


Over  the  period  1964  to  1967  the  University  of  Wisconsin 
continued  shallow  seismic  refraction  profiling  in  the  western 
half  of  Lake  Superior.  Anzoleaga  et  al,  (1969)  have  interpre¬ 
ted  the  shallow  structure  from  4  profiles  along  the  line  of 
the  1963  Knife  River,  Minnesota  to  Otter  Cove,  Ontario  (Figure 
12)  profile  using  the  gravity  data  of  Wold  (1968)  as  a  con¬ 
straint.  Figure  13  shows  the  model  derived  by  Anzoleaga  et  al, 
(1969).  The  shallow  structure  above  the  upper  refractor  of 
Smith  et  al,  (1966a,  1966b)  is  divided  into  three  layers  with 
P-wave  velocities  of  3.3  -  3.5  km/sec,  4.5  -  4.7  km/sec  and 

5.4  -  5.6  km/sec  respectively  which  pinch  out  under  the  Minnesota 
shore  of  the  lake  against  a  6.5  km/sec  dipping  layer.  The 
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structure  shows  a  thickening  of  all  three  layers  centered  some 
150  km  from  the  Knife  River  station  to  form  a  basin,  the  posi¬ 
tion  of  which  is  coincident  with  the  gravity  low  mid-way  be¬ 
tween  the  Bayfield  Peninsula  and  Isle  Royale  (Figure  17). 
Anzoleago  et  al,  (1969)  find  a  velocity  of  6.9  km/sec  for  the 
upper  refractor,  which  is  in  agreement  with  the  velocity  of 

6.8  km/sec  found  by  Smith  et  al,  (1966a,  1966b)  and  the  6.75 
6.85  km/sec  reported  by  O’Brien  (1968).  From  its  spatial 
configuration,  the  6.5  km/sec  layer  probably  represents  the 
Duluth  Gabbro.  The  remaining  three  layers  above  the  6.8 

6.9  km/sec  refractor  are  probably  the  Upper  Keweenawan  sedi¬ 
ments  and  the  Keweenawan  lavas . 

O’Brien's  model  gives  a  two-d: mensional  view  of  crustal 
structure  along  an  E— W  line  through  Lake  Superior.  In  an 
attempt  to  define  the  structure  of  the  Lake  Superior  basin  in 
three  dimensions,  we  have  plotted  fence-diagrams,  Figure  14, 
showing  crustal  structure  in  the  vicinity  of  Lake  Superior 
based  on  the  results  of  earlier  workers  in  this  area  and  the 
data  obtained  during  Project  Early  Rise.  To  the  south  of 
the  lake  the  work  of  Steinhart  et  al,  (1961)  indicates  total 
crustal  thicknesses  ranging  from  36  -  42  km  and  PR  velocities 
of  8.1  km/sec.  The  crustal  models  obtained  by  Steinhart  et  al 
(1961)  suggest  that  the  arrivals  can  be  explained  by  a  two- 
layer  crust  consisting  of  a  thin  layer  with  a  P-wave  velocity 
of  4,6  -  5.5  km/sec  overlying  a  thicker  layer  with  a  velocity 
of  6.3  -  6.4  km/sec.  However,  there  were  indications  in  North 
Central  Wisconsin  that  the  deeper  layer  should  be  considered 
as  two  layers  with  velocities  of  6.1  and  6,5  km/sec,  and  for 
the  Keweenan  Profile,  Steinhart  et  al,  (1961)  suggest  that  two 
layers  with  velocities  of  6.4  and  6.7  km/sec  are  indicated. 

To  the  southwest  of  Duluth,  Cohen  (1966)  and  Cohen  and  Meyer 
(1966)  found  crustal  thicknesses  of  46  km  along  the  axis  of 
the  Midcontinental  Gravity  High,  with  a  6.5  -  6.6/sec  layer 
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occurring  at  a  depth  of  1.8  km  and  an  apparent  velocity  of 

а. l  km/sec.  R.  P.  Meyer  (personal  communication,  1969)  has 
re-examined  this  data  and  believes  that  velocities  as  high  as 

б. 9  km/sec  may  be  appropriate  for  the  major  crustal  layer. 
Along  the  gravity  low  flanking  the  eastern  margin  of  the  high, 
Cohen  (1966)  and  Cohen  and  Meyer  (1966)  found  velocities  sim¬ 
ilar  to  those  found  along  the  axis  of  the  high,  however  the 
6.5  -  6.6  km/sec  layer  occurs  at  a  depth  of  4.6  km  and  the 
crustal  thickness  is  only  42  km.  In  an  analysis  of  travel¬ 
time  data  from  Early  Rise,  Mereu  and  Hunter  (1969)  found  that 
to  the  north  and  east  of  Lake  Superior  the  crust  is  30  -  35 

km  thick  under  the  Superior  Province  of  the  Canadian  Shield, 
thickening  to  greater  than  40  km  under  the  Churchill  Province. 
They  found  mean  velocities  along  a  profile  north  of  Lake 
Superior  for  Pg  =  6.30  km/sec  and  Pr  =  8.05  km/sec. 

Thus,  the  Lake  Superior  crustal  syncline,  on  the  basis 
of  the  refraction  data  presented  above,  is  cJ early  delineated. 
Figure  15  is  a  contour  map  of  the  Moho  depth  in  the  vicinity 
of  Lake  Superior  ba^ed  solely  on  the  refraction  data  dis¬ 
cussed  above.  The  high  mean  crustal  velocities,  6.7  -  6.9 
km/sec,  found  under  Lake  Superior  have  not  been  reported  from 
areas  outside  the  lake  with  one  exception.  Recent  examination 
of  refraction  data  recorded  along  a  path  entirely  within  the 
Midcontinental  Gravity  High  to  the  southwest  of  the  lake 
suggests  that  the  high  velocities  of  6.9  km/sec  at  depths  of 
8  km  beneath  the  lake  reported  by  Anzoleaga  et  al,  (1969) 
may  be  present  .along  axis  of  the  Midcontinental  Gravity  High 
(R.P.  Meyer,  personal  communication,  1969). 

Surface  Wave  Dispersion  Observations.  -  We  have  been  un¬ 
able  to  find  any  published  observations  of  surface  wave  dis¬ 
persion  across  Lake  Superior.  Brune  and  Dorman  (1963)  have 
published  curves  showing  Love  and  Rayleigh  dispersion  across 


the  Canadian  Shield  and  derived  a  model  (Pigure  4;  based 
on  their  observations.  However,  none  of  the  paths  studied  tv 
Brune  and  Dorman  cross  the  lake.  Additional  data  are  provided 
by  Pilant  (1967)  who  studied  Rayleigh  wave  dispersion  in  North 
America.  Pilant  measured  phase  and  group  velocities  using  the 
WW3SS  network  and  prepared  a  series  of  contour  maps  showing 
observed  fundamental  mode  phase  velocities  over  North  America 
as  a  function  of  period.  We  have  taken  the  values  a* 

Superior  and  plotted  them  together  with  the  theoretical  dis¬ 
persion  curve  predicted  by  Brune  and  Dorman's  model  A..  E  • 
Figure  16  shows  the  results. 

Pilant's  observed  values  fall  systematically  below  the 
theoretical  curve  at  periods  shorter  than  50  seconds.  However, 
we  can  attach  no  great  significance  to  this  since  the  deviation 
from  the  theoretical  curve  is  less  than  the  contour  interval 
of  Pilant's  data. 

Because  the  mean  crustal  thickness,  the  mean  crustal  P- 
wave  velocity  and  hence  by  inference  the  mean  shear  wave 
velocity  and  density  are  different  beneath  Lake  Superior  from 
those  incorporated  in  the  model  CANSD,  we  do  not  anticipate 
that  the  CANSD  phase  velocity  dispersion  curve  shown  in  Figure 
16  will  be  representative  of  the  Lake  Superior  region.  This 
can  be  seen  from  the  theoretical  dispersion  curve  also  shown 
in  Figure  16  which  was  computed  for  a  composite  Lake  Superior 
model  (Table  4 )  derived  from  seismic  refraction  and  gravity 
observations . 

Gravity  Data 

During  September  1964  the  Dominion  Observatory  of  Canada 
conducted  a  reconnaissance  underwater  gravity  survey  in  Lake 
Superior  (Weber  and  Goodacre,  1966).  Additional  data  along 
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the  south  shore  of  the  lake  and  in  Keweenaw  Bay  were  collected 
in  August  1964  by  the  University  of  Wisconsin  (Wold  and  Ostenso, 
1966).  Weber  and  Goodacre  have  computed  Bouguer  anomaly  values 
from  these  data  and  presented  a  contour  map  showing  Bouguer 
anomalies  in  the  Lake  Superior  region  at  10  mgal  intervals. 
w,e  anomalies  over  the  lake  were  joined  to  existing  gravity 
maps  cf  the  areas  surrounding  the  lake,  (Innes,  1960;  Woollard 
and  Joesting,  1964;  McConnell,  1966).  The  resulting  map 

(Weber  and  goodacre,  1966,  Figure  2)  covers  the  area  45  -  50°iJ, 

8  3  -  9  3°W . 

'  ince  1  the  University  of  Wisconsin  has  completed  a 
gravity  survey  of  Lake  Superior  on  a  5-mile  grid,  occupying 
some  1600  stations  (Wold,  1968).  Although  there  are  some 
minor  differences  in  detail  between  Wold’s  map  and  the  earlier 
results  of  Weber  and  Goodacre  (1966),  the  overall  features  are 
similar  (R.  Wold,  personal  communication,  1969), 

In  or  Jer  to  examine  the  relationship  between  the  gravity 
rield  at  Lake  Superior  and  the  surrounding  area  on  a  regional 
t asit  ,  we  have  compiled  a  Bouguer  anomaly  map  of  the  area  38  - 

,  /h  -  104°W  (Figure  17).  The  base  map  scale  was  1:2.5x10^ 
to  correspond  to  scales  of  geological,  tectonic,  and  gravity  maps 
of  the  United  States  published  by  the  U.  S.  Geological  Survey. 

The  contour  interval  is  10  mgal.  Th«  majority  of  the  data  are 
from  the  sources  quoted  by  Weber  and  Goodacre  (1966)  with  addi¬ 
tional  data  from  a  series  of  maps  compiled  by  the  Dominion 

Observatories  Branch,  Department  of  Mines  and  Technical  Surveys, 
Ottawa,  Canada. 

The  most  conspicuous  feature  on  the  gravity  map  is  the 
Midcontinental  Gravity  High  which  runs  northeast  from  Kansas 
through  Iowa  and  appears  to  terminate  south  of  the  Bayfield 
Peninsula  in  Wisconsin.  It  is  composed  of  four  elongate  areas 
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of  strong  positive  Bouguer  anomalies  flanked  by  regions  of 
pronounced  negative  anomalies.  Black  (1955)  has  interpreted 
these  anomalies  as  indicating  "an  echelon  rift  pattern  de¬ 
veloped  in  the  crust  during  late  Precambrian  time  and 
genetically  related  to  the  development  of  the  Lake  Superior 
geosyncline  during  the  Proterozoic’' .  However,  the  relative 
lows  which  separate  the  highs  could  equally  be  due  to  thinning 
of  the  perturbing  body  since  the  magnetic  anomalies  associated 
with  the  high  do  not  show  signs  of  appreciable  offset. 

Seismic  refraction  studies  (Cohen,  1966;  Cohen  and  Meyer, 
1966)  on  the  Midcontinental  Gravity  High  indicate  that  the 
crust  is  somewhat  thicker  (46  km)  than  in  the  surrounding 
areas  (42  km).  They  also  indicate  high  velocity  rocks  at 
shallow  depths  under  the  high  and  somewhat  deeper  under  the 
lows.  Cohen  and  Meyer  (1966)  found  rocks  with  P-wave  velocity 
of  6.5  km/sec  under  the  high.  However,  recent  data  suggest 
that  the  6.9  km/sec  rocks  found  beneath  Lake  Superior  (Anzoleaga 
et  al,  1969)  may  be  present  at  shallow  depths  under  the  high 
(R.  P.  Meyer,  personal  communication,  1969).  The  gravity  lows 
correlate  well  with  the  thick  accumulations  (up  to  5.0  km)  of 
low  velocity  sediments  found  by  Cohen  and  Meyer  (1966). 

The  seismic  studies  support  the  hypothesis  postulated  by 
Ihiel  (1956)  and  Craddock  et  al,  (1963)  that  the  positive  anom- 
lies  are  due  to  thick  accumulations  of  Keweenawan  volcanic 
rocks  in  a  synclinal  trough  which  has  been  uplifted  to  form 
a  horst  and  that  the  negative  anomalies  flanking  the  highs  are 
due  to  thick  accumulations  of  low  density  sediments. 

To  the  southwest  of  the  Bayfield  Peninsula  the  Bouguer 
anomalies  change  from  +60  mgal  to  -80  mgal  over  a  distance  of 
100  km.  In  this  region  the  synclinal  axis  of  the  Lake  Superior 
syncline  swings  E-v7  before  turning  NE-SW  and  passing  under  the 
lake  between  the  Bayfield  and  Keweenaw  peninsulas.  The  reversed 
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faults  which  bound  the  horst,  the  Douglas  Fault  to  the  north 

and  the  Lake  Owen  fault  to  the  south,  also  swing  E-W  in  this 
area  (Figure  18). 

The  large  negative  gravity  anomaly  over  the  Bayfield 
Peninsula  correlates  with  a  large  accumulation  of  Upper 

Keweenawan  sediments  whose  thickness  has  been  estimated  by 
White  (1966)  as  7  -  9  km.  To  the  east  of  the  Bayfield  Pen¬ 
insula  there  is  a  region  of  broad  positive  anomalies  which 
extends  from  the  Porcupine  Mountains  and  the  Keweenaw  Penin¬ 
sula  to  Isle  Royale.  To  the  west  of  the  Bayfield  Peninsula 

there  is  a  large  positive  anomaly  centered  over  the  Minnesota 
shore. 

Apart  from  three  small  closed  highs  centered  over 
Michipicoten  Island,  over  Marathon,  and  to  the  south  of  Isle 
‘ t*  Ignace,  the  remainder  of  the  anomalies  over  Lake  Superior 
are  nagative.  The  mean  Bouguer  anomaly  over  the  lake  is  -25  mgal 
as  compared  to  the  average  value  of  -45  mgal  over  the  surrounding 
highlands  (Weber  and  Goodacre,  1966). 

The  gravity  high  over  the  Minnesota  shore  of  Lake  Superior 
correlates  with  the  outcrop  of  the  North  Shore  Volcanics  which 
are  intruded  by  the  Duluth  Gabbro,  a  lopolithic  body  which 
from  geologic  evidence  thins  lakewards  from  ^he  shore.  If 
the  6.5  km/sec  layer  in  Figure  13  represents  the  Duluth  Gabbro, 
then  the  seismic  refraction  data  of  Anzoleaga  et  al,  (1969) 
are  consistent  with  this  interpretation.  According  to  White 
(1966a)  th  -  Keweenawan  rocks  to  the  east  of  Beaver  Bay  complex 
are  separated  from  those  under  the  rest  of  Lake  Superior  by  a 
narrow  ridge  in  the  pre-Keweenawan  surface  over  which  the  lavas 
thin.  However,  the  gravity  data  suggest  that  the  structure  is 
more  complicated  than  this. 

The  closed  circular  low  NW  of  the  Bayfield  corresponds  to 
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a  relatively  thick  sedimentary  basin  filled  by  low  velocity 
sediments  (Ansoleaga  nr  al.,  1969)  occupying  a  50  km  wide  de¬ 
pression  in  the  upper  sun  ace  of  the  Keeweenawan  lavas  which 
themselves  appear  to  occupy  a  trough  in  the  pre-Keweenawan 
surface . 

The  origin  of  the  positive  anomaly  between  the  Keweenaw 
Peninsula  and  Isle  Royale  is  more  uncertain.  There  are  steep 
gradients  in  the  Bouguer  anomalies  over  the  Keweenaw  fault 
and  the  Isle  Royale  fault.  These  faults  (Figure  18),  are 
reversed  or  thrust  faults  which,  by  analogy  with  the  structure 
southeast  of  Ashland,  suggest  that  the  lavas  occupy  a  synclinal 
trough  which  has  been  uplifted  to  form  a  horst.  However, 
O'Brien's  (1968)  interpretation  of  the  structure  along  the 
1963  main  line  profile  shows  a  decrease  in  thickness  of  the 
upper  layer  of  his  model  through  this  region.  This  is  some¬ 
what  surprising  since  the  profile  is  coincident  with  the 
inferred  synclinal  axis  at  this  point,  although  uplift  and 
erosion  of  the  syncline  as  suggested  above  could  account  for 
such  a  result. 


The  gravity  data  from  the  eastern  sector  if  the  lake  in¬ 
dicate  a  quite  different  structure  from  that  in  the  western 
sector.  The  broad  negative  anomalies  suggest  thick  accumula¬ 
tions  of  sediments.  However,  the  mean  anomaly  over  the  lake 
is  positive  with  respect  to  the  surrounding  areas  of  the  shi'*!^ 
(-25  mgals  vs.  -45  mgals)  which  together  with  the  known  out¬ 
crops  of  Keweenawan  lavas  on  Michipicoten  Island  and  the 
eastern  shore,  suggest  that  lavas  are  present  at  depth  beneath 
eastern  Lake  Superior.  This  interpretation  is  consistent 
with  the  seismic  data  discussed  earlier  and  the  aeromagnetic 
data  discussed  below. 


Aeromagnctlc  Data 


l 
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Wold  and  Oitenso  (1966)  and  Hfnze  al,  (1966)  hava 
compiled  a  total  magnetic  intanrity  map  oi  tht  Laka  Superior 
region  from  *ome  22,000  km  of  flight  linai  ipaced  at  9.6  km 
intervals.  Their  map  is  reproduced  as  Figure  19.  From  this 
data  Wold  and  Ostenso  (1966)  have  interpreted  the  geologic 
structure  under  the  western  half  of  the  lake  and  Hinze  at  al. 
(1966)  have  interpreted  the  structure  under  the  ea«iarn  sec¬ 
tion.  Together  their  interpretations  supplement  the  evidence 
from  seismic  and  gravity  data  above. 

In  the  western  section  of  the  lake  the  magnetic  anomalies 
indicate  a  broad  closed  low  which  extends  from  the  Bayfield 
Peninsula  to  the  north  of  the  Keweenaw  Peninsula.  Flanking 
-v,e  magnetic  low  to  the  north  and  south  there  are  a  series  of 
magnetic  highs.  Along  the  south  shore  of  the  lake  these  highs 
correlate  with  the  Keweenawan  lavas  in  the  upthrust  block  be¬ 
tween  the  Douglas  and  Lake  Owen  faults  and  the  Keweenawan 
lavas  exposed  to  the  north  of  the  Keweenaw  fault  on  the 
Keweenaw  Peninsula.  North  of  the  low,  the  positive  anomalies 
correlate  with  exposures  of  the  Duluth  Gabbro  and  the  North 
Shore  Volcanics.  The  positive  ar^maly  over  Isle  Royale  corre¬ 
lates  with  the  exposures  of  lakeward  dipping  Keweenawan  lavas. 

The  magnetic  low  is  evident  in  the  N-S  aeromagnetic  pro¬ 
files  across  the  western  half  of  the  lake  (Figure  20)  which 
are  remarkably  smooth.  The  broad  low  correlates  with  the  in¬ 
ferred  synclinal  axis  and  the  smoothness  of  the  negative 
anomalies  suggest  that  the  trough  is  occupied  by  a  thicl  sec¬ 
tion  of  sediments.  The  trend  of  the  anomalies  in  the  vicinity 
of  the  Bayfield  Peninsula  is  consistent  with  the  interpre¬ 
tation  of  White  (1966a)  who  suggested  that  the  horizon  which 
fcivts  rise  to  these  particular  anomalies  is  high  within  the 
sequence  of  Keweenawan  lavas  and  outlines  the  nose  of  a  NE 
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plunging  *ynclin«.  The  r*§ult§  of  thU  d»tAll*d  §urv*y 
support  Vfhitt’s  (1966*)  poitulat#  th*t  th*  l*v*>  thin  ovar 
«  11-3  ridg#  in  th#  pr*-K«w##n*w*n  *urf*ct  batvaan  th#  Bayfitld 
P#nin*ul*  #nd  Ba*v#r  B*y  which  Whit#  h#d  «ug£#it«d  on  tht  b#*i# 
w*  data  coll#ct«d  from  *l«v#n  i #connai#i*nc#  profil#*  *croi* 
th«  lak#. 

To  tht  «**t  of  th#  K#w*«niiw  Ptnintula,  th#  ch*r*ct*r  of 
th#  #noGs#li%5  chang#*  cmrhtdly .  Th#  anomali##  in  th.i  At**** 
fora  *  t«ri«s  of  clo§#d  highs  #nd  lows  (Figur#  21).  Hownvtr , 
th#  m*gn«tlc  low  in  th#  wtit^n  p*rt  of  th*  la>,«  c#n  b*  tr*c«d 
east  of  th*  K«w#«n*w  P#ninsul*  whtrt  it  swings  *pproxia*t*ly 
H-S  and  «xt*nds  into  th#  Horthtrn  P#ninsul*  of  Michigan  suggest¬ 
ing  that  th#  s«dia«ntary  basin  does  also.  Th*  magnetic  anomaly 
minimum,  B,  of  Figure  21  also  correlates  with  a  Bougu«r  gravity 
anomaly  minimum  (Bacon,  1967),  further  supporting  this  inter¬ 
pretation  of  th*  magnetic  anomaly  (Hinz#  *t  al,  1966 . ) 

The  positive  anomaly,  C,  in  Figure  21  correlates  with 
the  exposures  of  Keweenawan  volcanics  along  the  Keweenaw 
Peninsula  and  its  axis  also  swings  through  E-W  to  join  with 
the  Il-S  trending  axis  of  anomaly  D  which  is  a  continuation 
of  anomaly  C  (Hinze  et  al,  1966)  extending  south  at  least  as 
far  as  the  northern  shore  of  Lake  Michigan.  The  configuration 
or  anomaly  D  indicates  that  the  volcanics  dip  eastwards  and 
are  more  deeply  buried  to  the  south. 

Along  the  northern  and  eastern  shores  of  Lake  Superior 
there  are  a  series  of  six  positive  anomalies.  Anomaly  E  is 
a  continuation  of  the  anomaly  associated  with  the  Keweenawan 
lavas  outcropping  on  Isle  Royale.  The  configuration  of  the 
Anomaly  E  indicates  that  the  lavas  dip  SW  (Hinze  et  al,  19bO 
into  the  Keweenawan  syncline.  This  anomaly  is  separated  from 
anomaly  G  by  the  Isle  Royale  thrust  fault.  The  configuration 
of  anomaly  G  indicates  that  the  lavas  dip  toward  the  axis  of 


tht  basin  (Hint#  «?  «1,  1966).  Anomaly  r  correlates  with 

th«  outcrop  o I  Kawatnawan  lavas  on  Hichipicoten  Island,  which 
suggests  that  th#  source  of  the  anomaly  is  a  sequence  of  basic 
volcanics  dipping  westward  into  the  basin  (Hinzt  et  al.,  1966). 
To  the  north  anomaly  F  merges  into  another  magnetic  high,  anom- 
F  passes  into  the  E-W  trending  anomaly  H,  which  Hinze  et  alj 
1966  interpret  as  being  due  to  an  uplifted  block  of  volcanics 
overlain  by  a  thin  veneer  of  clastic  sediments.  The  southern- 
*©st  anomaly  I  correlates  with  the  Keweenawan  volcanics  of 
Haimainse  Point,  which  dip  at  30°  towards  the  lake  (Thompson, 
1953).  From  the  magnetic  data,  the  volcanics  appear  to  ex¬ 
tend  beneath  Wnitefish  Point.  Anomaly  J  is  similar  in  shape 
to  I  but  has  no  expression  in  the  surface  geology. 

The  aeromagnetic  profiles  also  indicate  the  extent  of 
faulting  within  Lake  Superior.  We  have  combined  Wold  and 
Ostenso*s  interpretation  of  the  structure  of  the  western  half 
of  Lake  Superior  with  Hinze  et  al.'s  interpretation  of  the 
structure  of  the  Lake  Superior  Basin.  Firstly,  in  the  western 
half  of  the  lake  the  magnetic  data  clearly  delineate  the  thrust 
faults  which  bound  the  upthrusted  central  block  of  the  syncline. 
Secondly,  although  the  aeromagnetic  data  indicate  that  the  syn¬ 
cline  is  continuous  under  the  eastern  part  of  the  lake,  they 
suggest  that  the  structure  is  somewhat  different.  The  thrust 
faults  which  parallel  the  synclinal  axis  in  the  western  part 
of  the  lake  give  way  to  a  series  of  faults  which  cut  across 
the  structure  east  of  the  Keweenaw  Peninsula.  And  thirdly, 
the  magnetic  data  suggest  that  either  the  Douglas  Fault  dies 
out  beyond  central  Bayfield  County  or  that  it  continues  at 
great  depth  and  has  small  displacements  (Patenaude,  1966). 

There  is  no  evidence  from  the  aeromagnetic  data  to  suggest 
that  the  Isle  Royale  Fault  is  an  extension  of  the  Douglas 
Fault . 


-26- 


Heat  Flow  Observations 


The  Carnegie  Institution  of  Washington  began  a  regional 
heat  flow  survey  of  Lake  Superior  in  1966  which  they  completed 
in  1967.  The  results  of  this  survey  are  shown  in  Figure  22 
which  is  taken  from  Steinhart  et  al. ,  (1968).  Although  the 
moan  value  of  the  measurements  (1.0  +  0.2  ycal/cm2-sac )  is  in 
good  agreement  with  the  mean  heat  flow  over  the  Canadian 
Shield,  the  range  of  the  observations  (0.5  to  1.5  y cal/cm2-sec) 
is  quite  large;  rapid  variations  in  heat  flow  taking  place 
over  a  few  tens  of  kilometers.  Steinhart  et  al.,  (1968)  con¬ 
clude  that  these  variations  are  not  due  to  measurement  errors 
which  only  account  for  deviations  of  10%  nor  are  they  due  to 
errors  related  to  topography,  etc.,  which  contribute  less  than 
10%  to  the  observed  variations.  If  these  conclusions  are 
correct,  then  the  variations  in  heat  flow  must  reflect  varia¬ 
tions  in  structure  beneath  the  lake. 

Steinhart  et  al.,  (1968)  divided  the  heat  flow  patterns 

in  Lake  Superior  into  three  anomalies  (see  Figure  22).  The 

first  of  these,  the  "north  shore  low"  with  heat  flow  values 

2 

of  0.5  -  0.9  ycal/cm  -sec  correlates  with  known  exposures  of 

Keweenawan  lavas  'tnd  with  near  surface  occurrences  of  the  lavas 

as  inferred  from  gravity  and  magnetic  data.  In  the  eastern 

and  central  part  of  the  lake  the  "central  high"  with  heat  flow 

2 

values  of  1.0  -  1.5  ycal/cm  -sec  correlates  with  the  gravity 
and  magnetic  anomalies  associated  with  the  Upper  Keweenawan 
sediments.  In  the  western  part  of  the  lake  the  heat  flow 
high  correlates  only  with  the  magnetics  (Steinhart  et  al, 

1968).  The  third  anomaly,  the  "central  low"  is  characterized 
by  heat  flow  values  between  0.8  and  1.0  ycal/cm  -sec.  Al¬ 
though  there  are  lavas  exposed  along  the  Keweenaw  Peninsula 
where  the  heat  flow  is  low  there  are  no  gravity  or  magnetic 
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highs  which  correlate  with  the  closed  heat  flow  low  to  the 
northeast  of  the  peninsula. 

The  broad  correlation  of  heat  flow  anomalies  in  Lake 
Superior  with  the  near  surface  geology  suggests  that  the 
lateral  variations  in  lithology  of  the  upper  10  km  of  the 
crust  may  be  sufficient  to  explain  the  observed  differences 
in  heat  flow.  Steinhart  et  al,  (1968)  state  that  although 
such  a  model  could  qualitatively  explain  the  observations, 
quantitatively  it  would  require  that  the  Oronto  group  of 
sediments  have  the  equivalent  radioactivity  of  a  granite. 
Despite  this  objection,  Steinhart  et  al.  (1968)  ascribe 
the  neat  flow  anomalies  to  variations  _n  composition  of  the 
upper  10  km  of  the  crust,  claiming  that  variations  in  total 
crustal  thickness  along  are  insufficient  to  account  for  the 
wide  range  of  the  observations. 

The  fact  that  the  brc \d  heat  flow  anomalies  seem  to 
correlate  with  the  synclinal  structure  of  Lake  Superior, 
with  high  values  occurring  along  the  axis  of  the  syncline 
and  low  values  on  the  flanks  can  be  explained  simply  by  the 
geometry  of  the  synclinal  structure  in  the  crust.  One  can 
show  that  the  synclinal  axis  should  exhibit  the  highest 
heat  flow  values  assumming  a  laterally  uniform  temperature 
at  depth  (the  Moho,  say)  below  the  syncline.  The  maximum 
lateral  variation  in  heat  flow  depends  on  the  "wave-length" 
ox  the  syncline  in  cross-section  and  the  "amplitude"  of  its 
undulation.  The  same  effect  would  explain  the  correlation 
o..  smaller  anomaly  highs  with  topographic  lows  in  the  lake 
bottom  and  vice  versa  but  not  the  observed  magnitude  of 
the  anomalies. 
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To  determine  the  seismicity  of  the  Lake  Superior  region, 
we  searched  the  Preliminary  Determination  of  Epicenters  cards 
published  by  the  U.  S.  Coast  and  Geodetic  Survey  for  events 
occurring  in  the  area  45  -  52°N,  80  -  100°W  for  the  period 
January  1960  through  January  1968.  We  found  only  one  small 
event  located  at  48.1°N  00.0°W.  However,  in  the  literature 
we  found  reports  of  two  other  significant  earthquakes  (July 
26,  1905  and  May  26,  1906)  which  occurred  before  the  in¬ 
stallation  of  a  sensitive  seismic  recording  network  in  the 
United  States  (Eppley,  1965).  Both  events  were  located  near 
Lake  Superior  on  the  Keweenaw  Peninsula.  The  July  26,  190 > 
event  was  felt  over  the  entire  Keweenaw  Peninsula  and  caused 
damage  at  several  places.  The  May  26,  1906  event  caused  dam¬ 
age  and  significant  elevation  changes  at  the  Atlantic  Mine  on 
Keweenaw  Peninsula;  i4*  affected  an  area  about  160  km  in  diam¬ 
eter.  In  addition  there  are  two  earthquakes  located  some¬ 
where  in  central  Minnesota  (Figure  23)  which  may  also  be 
associated  with  Lake  Superior  tectonics.  The  best  documented 
event  occurred  on  3  September  1917  and  caused  damage  at 
Brainard  and  Staples,  Minnesota.  The  other  in  the  same  gener¬ 
al  region,  occurred  during  186'  but  little  is  known  about  - t. 
Both  events  are  located  in  the  vicinity  of  the  Midcontinental 
Gravity  High  where  it  passes,  under  the  southwest  tip  of  Lake 
Superior . 

Thus  in  recent  years  the  Lake  Superior  region  as  been  rel¬ 
atively  aseismic.  However,  the  observed  faulting  discussed 
above  suggests  that  the  area  must  have  been  tectonically  active 
in  the  geologic  past. 


A  Geophysical  Model  for  the  Crust  Beneath  the  Early  Rise  Shot 
Point 


During  Project  Early  Rise,  38  shots  each  consisting  of 
4,844  kg  (0.005  KT)  of  DuPont  Nitramon  WW(EL)  (Borcherdt 
et  al,  1967)  were  detonated  on  the  bottom  of  Lake  Superior 
during  July,  1966.  Borcherdt  et  al.,  (1967)  located  the  ac¬ 
tual  shot  points  from  the  travel-times  of  hydro-acoustic  waves 
recorded  by  six  hydrophone  stations  around  Lake  Superior.  Fig¬ 
ure  24  shows  the  computed  shot  points.  The  shots  fall  within 
a  radius  of  0.7  km  from  the  mean  shot  point  (47.55°N,  88.94°W) 
with  the  exception  of  four  shots.  Of  these  four,  three  are 
within  approximately  1  km  of  47.55°N,  88.94°W;  the  fourth, 
shot  #1  is  located  some  3.3  km  to  the  east  of  Ihe  mean  shot 
point.  The  coordinates  of  the  38  shot  points  and  the  origin 
times  are  included  in  Appendix  I. 

The  University  of  Wisconsin  conducted  continuous  sub¬ 
bottom  profiling  surveys  in  western  Lake  Superior  (Figure  25) 
during  1966.  Figure  26  is  a  profile  through  the  Early  Rise 
shot  point  (Personal  communication,  R.  J.  Wold,  1968)  obtained 
using  a  salt  water  bag  sparker  system  which  excites  80  -  150  Hz 
acoustic  waves.  The  profile  shows  that  the  lake  bottom  in  the 
immediate  vicinity  of  the  Early  Rise  shot  point  is  nearly 
horizontal  with  little  or  no  relief.  The  two  way  travel- 
times  of  the  acoustic  waves  reflt^ted  from  the  bottom  at  the 
Early  Rise  shot  point  is  0.275  sec. 

To  convert  the  two-way  travel- times  to  a  water  depth,  we 
need  to  know  the  velocity  of  sound  in  the  lake.  Browning  et  al, 
(1968)  have  published  velocity-depth  profiles  obtained  at  three 
locations  in  eastern  Lake  Superior  during  September  1967.  These 
profiles  show  that  the  velocity  is  almost  constant  at  1422  m/sec 
for  depths  greater  than  50  m.  Between  50  m  and  the  surface 
the  velocity  increases  approximately  linearly  to  1445  m/sec. 
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Therefore,  considering  the  lake  as  a  two  layered  model,  a 
50  m  thick  layer  with  velocity  of  1434  m/sec  overlying  a  lower 
layer  with  a  velocity  of  1422  m/sec,  we  calculated  the  depth 
of  the  lake  beneath  the  Early  Rise  shot  point  to  be  196  m. 

This  result  agrees  closely  with  the  depth  of  approximately 
190  m  reported  by  Borcherdt  et  a]  ,  (1967). 

The  University  of  Wisconsin  (Anzoleaga  et  a],  1969)  has 

also  conducted  shallow  seismic  refraction  profiling  through 
the  Early  Rise  shot  point  (Figure  12).  Figure  27  shows 
the  first  arrivals  from  the  reversed  profile.  The  close 
agreement  between  the  apparent  velocities  of  the  four  sets 
of  arrivals  recorded  from  the  eastern  shot  point  with  those 
from  the  western  shot  point  suggests  that  the  shallow  layers 
are  very  nearly  horizontal.  Figure  13  shows  Anzoleaga  et  a] 
interpretation  of  the  shallow  structure  beneath  the  western 
half  of  Lake  Superior. 

T  e  velocity  of  6.9  km/sec  for  the  first  major  crustal 
layer  is  somewhat  higher  than  the  velocity  reported  by  Smith 
et  al,  (3966a, b),  Berry  and  West  (1966ab)  and  O’Brien  (1968). 
Since  it  is  unlikely  that  the  P-wave  velocity  is  constant 
in  the  major  crustal  layer,  but  rather  increases  with  depth, 
we  will  adopt  6.8km/sec  for  the  velocity  of  this  layer,  a 
value  which  represents  a  mean  of  the  reported  measurements. 

The  choice  of  a  thickness  for  the  major  crustal  layer  also 
presents  some  difficulties.  O’Brien's  (1968)  model  (Figure 
10)  indicates  that  the  crust  is  approximately  40  km  thick, 

100  km  to  the  west  of  the  Early  Rise  shot  point  and  50km 
thick,  100  km  east  of  the  shot  point.  To  the  west  of  the  shot 
point,  the  Moho  dips  7°E.  However,  to  the  east  of  the  shot 
point  the  dip  is  highly  irregular.  Since  we  are  attempting 
to  construct  a  plane-parallel  layered  model,  we  will  adopt  a 
total  crustal  thickness  of  48  km,  implying  that  the  major 
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crustal  layer  is  40km  thick  beneath  the  Early  Rise  shot  point. 
The  velocity  in  the  Lake  Superior  region  is  well  defined  at 
8.1  km/ sec. 

The  shear  wave  velocities  for  the  crust  beneath  Lake 
Superior  are  less  well  known.  Lewis  and  Meyer  (1968)  report 
shear  wave  arrivals  with  apparent  velocities  of  3.65  km/sec 
and  3.80  km/sec  recorded  on  the  Keweenaw  Peninsula  during  the 
1903  Lake  Superior  Seismic  Experiment  from  shots  fired  in  the 
western  half  of  the  lake.  In  their  model,  Lewis  and  Meyer 
(1968)  assign  a  shear  wave  velocity  of  3,65  km/sec  to  the  upper 
10  km  of  the  crust.  However,  this  is  almost  certainly  too  high 
for  the  uppermost  layers  of  the  model  proposed  by  Anzoleaga 
et  al  (1969).  Therefore,  in  constructing  our  models 
(Table  4  )  we  have  calculated  the  shear  wave  velocities  in 
layers  1,  2,  and  3  from  the  compressional  wave  velocities 
assuming  a  value  of  Poisson’s  ratio  a  -  0.25.  Since  Poisson's 
ratio  may  be  somewhat  greater  than  this  in  the  shallow  layers,’ 
th°ir  true  shear  velocities  may  be  lower  than  computed.  The 
velocity  of  Sn  is  4.72  km/sec,  taken  from  Lewis  and  Meyer's 
model.  This  velocity  for  S^  is  in  agreement  with  Brune  and 
Dorman's  (1963)  observations  in  the  Canadian  Shield. 

To  complete  the  model  we  need  to  know  the  density  distri¬ 
bution.  Anzoleaga  et  al,  (1969)  constructed  their  model  (Fig¬ 
ure  13)  from  the  P-wave  arrival  data,  and  constrained  their 
results  to  fit  the  gravity  data  of  Wold  (1968).  The  densities 

of  layers  1,  2,  3  are  similar  to  the  densities  of  Jacobsville 

3  3 

Sandstone,  2.37  gm/cm  ;  Oronto  Group,  2.62  gm/cm  *,  and  mafic 

3 

rocks,  2.92  gm/cm  ;  obtained  empirically  by  White  (1966)  which 
are  in  general  agreement  with  the  results  of  earlier  workers 
(e.g.  Thiel,  1956).  The  density  of  the  6.9  km/sec  layer  is 
3.10  gm/cm  (Anzoleaga  et  a]  ,  1969).  When  the  gravity  anom¬ 
aly  of  the  crustal  model  ERLS111  is  calculated  with  respect 
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to  CANSD,  a  density  of  3.10  gm/cm  for  the  40  km  thick  layer 

gives  an  anomaly  of  +221  mgals  which  does  not  agree  with  the 

observed  anomaly  of  approximately  +60  mgals  with  respect  to 

the  Canadian  Shield  near  Lake  Superior.  That  is,  from  the 

Bouguer  gravity  map  (Figure  17),  the  Bouguer  anomaly  at  the 

Early  Rise  Shot  Point  is  2  0  mgals  while  the  mean  value  in  the 

area  surrounding  Lake  Superior  is  a  -60  mgal.  If  the  density 

3 

of  the  40  km  thick  layer  is  reduced  to  3.0  gm/cm  (or  alter¬ 
nately  if  the  density  in  some  40  km  of  the  shield  model  is  in¬ 
creased  by  0.1  gm/cm3)  then  the  anomaly  relative  to  CANSD  is 
+57  mgals  in  close  agreement  with  the  observed  anomaly.  Thus, 
the  relative  anomaly  between  the  observed  Bouguer  gravity  at 
the  Early  Rise  shot  point  and  the  area  surrounding  Lake  Supei  ior 
is  in  good  agreement  with  the  Bouguer  anomaly  between  the  new 
model  ERLS112  and  CANSD  assuming  that  the  mean  density  of  the 
40  km  thick  layer  is  3.0  gm/cm  .  Throughout  we  have  adopted 

Lewis  and  Meyer’s  (1968)  value  of  p  for  the  upper  mantle,  3. 

3 

3.42  gm/cm  . 

In  addition  to  models  ERLS111  and  ERLS112,  we  have  con¬ 
structed  a  third  model,  ERLS121  assuming  a  Poisson’s  ratio  of 

3 

0.25  throughout  the  crust  and  a  density  of  3.10  gm/cm  for  the 
6.8  km/sec  layer.  Of  the  three  models,  ERLS112  is  probably 
the  most  representative.  Figure  16  shows  the  theoretical 
fundamental  mode  Rayleigh  wave  dispersion  waves  computed 
from  these  models. 


DISCUSSION 


The  geophysical  investigations  of  the  properties  of  the 
each's  crust  beneath  Lake  Superior  over  the  past  decade  have 
shown  that  Irving's  (1883)  hypothesis  that  the  crust  beneath 
Lake  Superior  has  been  downwarped  to  form  a  broad  syncline  is 
correct.  However,  the  same  investigations  have  also  shown 
that  crustal  structure  is  considerably  more  complex  than 
was  originally  supposed.  Neverthelf ;s  there  is  a  strong 
correlation  between  the  inferred  geology,  structure  and  geo¬ 
physical  properties  of  the  shallow  part  of  the  earth's  crust 
in  the  Lake  Superior  region.  For  the  deeper  part,  the  results 
from  the  1963  Lake  Superior  Seismic  Experiment  suggest  that 
the  structure  of  the  lower  crust  is  considerably  more  complex 
than  was  originally  supposed.  In  view  of  the  high  seismic 
velocities  for  the  lower  crust  observed  by  all  investigators, 
the  deeper  structure  of  the  Lake  Superior  region  probably  is 
not  representative  of  the  Canadian  Shield  as  a  whole. 

In  the  western  half  of  Lake  Superior,  the  central 
part  of  the  syncline  has  been  upthrusted  to  form  a  horst. 

This  portion  of  the  structure  appears  to  be  genetically  re¬ 
lated  to  the  Midcontinent a?  Gravity  High.  There  is  an  ex¬ 
cellent  correlation  between  the  surface  geology,  gravity, 
magnetics,  heat  flow  and  the  shallow  seismic  refraction  re¬ 
sults.  All  suggest  that  the  Keweenawan  lavas  and  sediments 
are  contained  in  a  synclinal  trough  approximately  10  km  deep, 
the  axis  of  which  is  orientated  in  a  NE  -  SW  direction.  The* 
results  of  Anzoleaga  et  al  ,  (1969)  suggest  that  many  of  the 
gravity  and  magnetic  anomalies  in  the  region  particularly  those 
where  large  gradients  are  involved,  can  be  explained  by  lateral 
variations  in  structure  of  the  upper  10  km  of  the  crust. 
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To  the  east  of  the  Keweenaw  Peninsula  the  shallow 

,  ,  nattern  of  the  anomalies  change. 

:rustal  structure  an  .  from  NE_SW  west  of  the  penin- 

'he  aXiS  °f  f*  Sytn^finit  '  The  positive  gravity  anomalies 

““ » *  ,«U" -hi°h fo™ ,h* . 

,rth  steep  flank  gl  ^  vioinity  0f  the  Keweenaw  Penxnsula. 
gravity  Hig  J  the  apparent  disappearance  of  the  high 

in  crustal  structure 

the  eastern  section  of  the  lake. 

,  +v,_  lake  there  is  less  shallow 
in  the  eastern  section  o^he^.  Qf  ^  anomalies  change. 

seismic  refraction  ata  an  ^  sufficient  to  indicate  that  the 

However  the  refrac  ion  ^  thickness  and  velocity. 

near  surface  lay  ^  ^  heat  floW)  magnetic  and 

These  variations  g  structure  of  the  shallow  part 

gravity  data  although  the  magnetic  data 

itdii:;usor:r -  -  ::nrlting- 

If  the  magnetic  highs  in  the  southeaster n^par^,^  ^ 

superior  ^^^^^Keweenawan  sediments,  then  the 

and  the  magneti  persists  to  the  south- 

n •  .1  ctructure  probably  persib 

shallow  syncli  #  continuance  is  reasonable 

rl^l/lro^^n: tentatively  supported  by  the  gravity 
unT,k  jf  Bacon  (19  57  ). 

f  the  deeper  structure  beneath  Lake  Superior 
The  variation  16  we  show  a  contour  map 

is  less  well  known.  8  Superior  region  based 

of  the  depth  to  the  Moho  m  the  ™*J***  c  refraction 

solely  on  a  compilation  of  structure  depends  heavily 

data.  This  interprets  ion  crustal  structure 

on  the  results  of  O’Brien  (1968)  for 
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beneath  the  main  line  of  the  1963  Lake  Superior  Seismic  Experi- 
ment.  We  chose  to  incorporate  O’Brien’s  structure  rather 
than  either  Berry  and  West's  (1966a, b)  model  or  Smith  et  ala’s 
(1966a, b)  model  because  the  latter  data  do  not  fulfill  the 
requirements  of  the  time-term  method  which  was  used  to  de¬ 
rive  them.  Interpretation  of  the  Lake  Superior  travel-time 
data  by  this  method  results  i  1  a  significant  exaggeration  of 
the  Moho  structure.  Additional  support  of  the  choice  of  the 
O'Brien  (1968)  structure  is  that  it  satisfactorily  explains 
the  observations  of  Bednarek  and  Meyer  (1966).  Bednarek  and 
Meyer  (1966)  plotted  reduced  travel-times  to  the  WWSSN  station 
MDS  (t-A/8.1)  seconds  minus  the  shot  point  time-term  reported 
by  Berry  and  West  (1966a)  as  a  function  of  distance  from  the 
Zipgun  station  and  found  that  the  difference  remained  constant 
at  about  3.4  seconds  for  stations  less  than  375  km  from  Zip- 
gun.  For  the  more  distant  stations,  i.e.,  the  stations  in  the 
western  half  of  the  lake  where  the  most  rapid  variations  in 
crustal  thickness  occur,  the  differerces  increase  with  range 
attaining  a  maximum  value  of  5.0  seconds  at  the  western  end 
of  the  lake.  If  instead  of  the  time-term  for  the  shot  point, 
the  delay  time  for  the  ray  offset  to  west  of  the  shot  point 
is  subtracted  from  the  reduced  travel  time,  the  difference 
remains  constant  and  is  independent  of  range  from  Zipgun 
(O'Brien,  1968). 

Figure  15  shows  that  the  crust  beneath  Lake  Superior 
attains  a  thickness  greater  than  50  km.  However,  the  figure 
also  shows  that  these  unusually  large  crustal  thicknesses  are 
probabDy  confined  to  the  eastern  section  of  Lake  Superior. 

The  crust  thins  rapidly  northwest  of  Lake  Superior,  and  has 
a  normal  thickness  for  t ha  Canadian  Shield  (  *35  km)  to  the 
northeast  of  the  lake  in  Ontario.  Earlier  refraction  work 
(Steinhart  et  al,  1961)  indicate  that  the  crust  thins  to  the 
south  of  Lake  Superior  also. 
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The  unusually  high  velocity  of  the  major  crustal  layer 
reported  by  several  investigators  (Smith  et  al,  1966a, b; 
O'Brien,  1968;  Anzoleaga  et  al^  1969)  appears  to  be  largely 
confined  to  the  crust  beneath  Lake  Superior.  Seismic  re¬ 
fraction  profiles  in  the  area  around  the  Lake  indicate  that 
the  corresponding  layer  has  a  velocity  of  6.3  to  6.4  km/sec, 
which  is  similar  to  values  found  further  north  in  the  shield. 
The  only  evidence  for  the  existence  elsewhere  in  thi  region 
of  velocities  as  high  as  6.9  km/sec  at  depths  as  shallow  as 
8  km  deep  in  the  crust  is  data  obtained  during  a  recent  re¬ 
examination  of  the  travel- time^.  from  shots  in  Lake  Superior 
recorded  along  a  profile  entirely  within  the  Midcontinental 
Gravity  High  (R.P.  Meyer,  personal  communication  1969). 

Thus  the  occurrence  of  the  unusually  high  velocities  appears 
to  correlate  with  the  occurrence  of  the  Keweenawan  lavas. 
However,  Halls  (1969)  has  examined  compressional  wave  veloci¬ 
ties  in  the  Portage  Lake  lavas  and  demonstrated  that  sample 
velocities  in  the  slightly  metamorphosed  members  at  2  Kb 
pressure  are  insufficient  to  explain  the  high  velocity  (6.6  - 
6.9  km/sec)  of  the  upper  refractor.  But,  Halls  found  that 
amphihole-bearing  volcanics,  which  represent  a  higher  metamor- 
phic  grade,  taken  from  near  the  base  of  the  Keweenawan  extru¬ 
sive  sequence,  had  sample  velocities  of  6.7  km/sec  at  2Kb  and 
thus  possibly  represent  the  top  surface  of  the  upper  refractor. 
This  suggests  that  possible  sources  for  these  lavas  are  magmas 
which  pervaded  the  deeper  layers  of  the  crust.  However,  more 
conclusive  data  will  be  required  to  validate  this  hypothesis. 

Based  on  all  the  geophysical  and  geologic  data  discussed 
above,  we  have  constructed  a  composite  plane-layered  model 
ERLS112  (Table  4  )  which  we  believe  is  representative  of  the 
crustal  structure  beneath  the  Early  Rise  shot  point  as  it  is 
presently  known.  Clearly  this  model,  as  any  plane-layered 
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crustal  model  must,  represents  a  major  simplification  of  the 
true  crustal  structure  in  what  is  clearly  a  structurally 
complex  area.  Fur'  lermore,  the  shear  velocities  and  densi¬ 
ties  in  this  model  are  indirectly  inferred  estimates  rather 
than  established  quantities.  Nevertheless,  Model  ERLS112  will 
serve  as  a  useful  basis  for  comparison  of  the  structure  beneath 
the  Early  Rise  shot  point  and  that  at  Yucca  Flat,  NTS. 
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EPICENTER  LOCATIONS,  TRAVEL-TIMES,  AND  UPPER  MANTLE  STRUCTURE 

FOR  LAKE  SUPERIOR  REGION 


In  this  section  of  the  report  we  deal  with  the  observed 
travel-times  from  shots  in  Lake  Superior  which  were  recorded 
by  stations  at  regional  and  teleseismic  (A>1000  km)  distances. 
Because  the  1963  Lake  Superior  Seismic  Experiment  was  primarily 
intended  to  establish  th*  crustal  structure  beneath  the  lake 
itself,  there  was  no  extensive  station  coverage  at  large  dis¬ 
tances;  consequently,  we  will  limit  our  discussion  of  travel- 
times  to  data  obtained  from  the  October  1964  series  of  shots 
and  project  Early  Rise  in  1966  for  which  many  recording  sta¬ 
tions  were  available*  The  October  1964  series  of  shots  affords 
an  excellent  opportunity  to  assess  our  ability  to  detect  and 
locate  events  in  Lake  Superior  since  five  of  the  ten  shots 
were  fired  in  essentially  the  same  location,  whereas  the  re¬ 
maining  five  were  detonated  at  neighboring  points.  During 
project  Early  Rise,  38  shots  were  detonated  at  essentially  a 
common  location,  the  recordings  being  made  over  long  profiles 
which  for  the  most  part  radiated  from  the  shot  point  (Figure  2). 
Inversion  of  the  travel-times  and  amplitude  data  recorded 
along  these  profiles  yields  a  series  of  velocity-depth  struc¬ 
tures  for  the  upper  mantle  for  eastern  North  America. 


Location  of  the  October  1964  Shots  Using  the  LRSM  and  Observa- 
tory  Networks 


In  October  1964,  the  U.  S. Geological  Survey  fired  a 
series  of  ten  shots  in  western  Lake  Superior.  The  shots  ranged 
in  five  sizes  rrom  907kg  (0.001  KT)to  9070  kg  (0.01  KT),one 
shot  of  each  size  being  detonated  at  a  common  location.  The 
locations,  origin  times,  and  shot  sizes  were  originally  re¬ 
ported  by  Roller  and  Jackson  (1966).  As  reported,  the  loca- 
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tions  for  shots  46A  and  47A  are  in  error.  The  correct  locatioi  s 
determined  from  hydrophone  data  (W.  H.  Jackson,  personal  com¬ 
munication,  1969)  are  given  in  Table  A-2  of  Appendix  I,  together 
with  the  water  depth  at  each  shot  point  (J.  C.  Roller,  personal 
communication,  1969)  as  determined  by  the  firing  ship's  depth 
finder . 

We  determined  the  arrival-times  and  amplitudes  at  those  LRSM 
stations  and  observatories  which  recorded  these  shots ;  these  ob¬ 
servations  are  tabulated  in  Tables  A- 3  through  A-12  of  Appendix  I 
(Basic  Observational  Data).  The  arrival- times  were  obtained 
by  visually  cross-correlating  the  trace  from  each  shot  with 
the  trace  from  shot  49A  at  each  station.  For  the  smaller  shots, 
we  were  not  able  to  discern  any  arrival  at  many  of  the  more  dis¬ 
tant  stations.  However,  for  completeness,  these  stations  are 
included  in  Tables  A-3  through  A-12  of  Appendix  I.  Figure  28 
shows  the  distribution  of  LRSM  stations  which  recorded  the 
October  1964  shots. 

From  the  travel-times  in  Tables  A-3  through  A-12  of  Appendix  I, 
we  located  the  10  shots  with  the  computer  program  SHIFT  (Chiburis, 
1968).  These  locations  are  given  in  Table  6.  This  program  minimizes 
the  relative  travel-time  anomaly  errors  for  a  set  of  stations 
rather  than  travel-time  residual  errors  which  are  used  in  the  pro¬ 
gram  LOCATE  (Flinn,  1963;1965).  The  travel-time  anomaly  at  a 
station  i  relative  to  station  j  may  be  defined  as  (Chiburis,  1968): 


Ti-Tj’Hi  +  Hj  +  ei  “  ej 


(1) 


where  T  is  the  observed  travel— time,  H  is  the  expected  travel- 
time  computed  from  the  Herrin  68  travel-time  tables,  and  e  is 
an  ellipticity  and  elevation  correction.  The  anomaly  error 
is  then 
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dA  •/.  =  A  •  ,  .  -  A.,.  (2) 

l/j  l/j  l/j 

where  A.^,.  is  the  previously  determined  relative  anomaly  for 

the  region  in  which  the  event  occurs,  and  the  quantity  to  be 
minimized  becomes 


E •  /  •  *  dA •  /  •  dA • / • 

1/ D  -l/D  1/] 


where  dAj^  is  the  theoretical  relative  anomaly  given  by 


dA!  , . 
1/] 


3  A .  ,  .  3  A .  ,  . 

— ill  dA  +  d<fr  +  c. 


Cj  is  the  mean  error  at  the  reference  station  and  A  and  <£  are 
the  longitude  and  latitude  of  the  trial  epicenter  respectively. 

Figure  29  shows  the  epicenters  of  the  10  shots  as  computed 
from  the  whole  LRSM  and  observatory  network  without  the  anomaly 
corrections  The  computed  epicenters  for  all  10  shots  lie 

to  the  east  of  the  respective  shot  points.  The  distance  be¬ 
tween  the  epicenter  and  shot  point  is  outside  the  95%  confi¬ 
dence  ellipse  for  the  computed  epicenter. 

Since  there  are  no  values  of  A^^  currently  available  for 

the  Lake  Superior  region,  we  re-computed  the  epicenters  for 

nine  of  the  shots  using  travel-time  anomalies  from  shot  49A 

relative  to  station  WMO  (Table  5)  as  estimates  of  A^  (Tabic  6). 

The  recomputed  epicenters  (Figure  29)  with  the  exception  of  46A, 

still  lie  to  the  east  of  the  respective  shot  points.  However, 

the  mean  location  error  is  now  only  2.1  km,  and  the  shot  point 

falls  within  the  95%  confidence  ellipse  in  all  cases.  The  mean 

2 

area  of  the  95%  confidence  ellipses  is  now  31.2  km  compared  tc 
465.0  kin  without  anomalies  for  the  same  network  of  stations 
(Table  6).  It  is  of  interest  to  note  that  Figure  29  suggests  that, 
knowing  the  location  of  shot  49A  and  hence  the  displacement  vector 
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of  the  epicenter  computed  without  anomaly  corrections,  we  could 
hti^c  improved  the  location  of  the  remaining  nine  shots  by  simply 
graphically  subtracting  the  displacement  vector  for  snot  4SA 
from  their  computed  locations.  This  is  a  consequence  o,  hav¬ 
ing  basically  the  same  set  of  receiving  stations,  hence  the 
same  bias,  in  each  case. 

.he  +rave^-times  to  the  ten  Arizona  stations  (see  in- 
ser  4 ,  Figure  28)  are  systematically  late  with  respect  to  the 
Herrin  68  travel-times .  Since  these  stations  are  located 
over  a  small  distance  and  azimuth  range,  they  may  be  expected 
to  Lias  the  locations.  We  recomputed  the  epicenters  for  all 
ten  shots  again  using  the  LRSM  and  observatory  network,  but 
Wi' h  only  two  of  the  stations  shown  in  the  insert  to  Figure 
28 ’  iFu  and  FN-UT.  The  results  are  shown  in  Table  7.  Without 
anomalies,  the  mean  location  error  is  reduced  to  12.6  km  and 
the  azimuth  of  the  computed  epicenter  from  the  actual  shot 
P°int  increased  by  around  5°.  The  mean  area  of  the  96%  con¬ 
fidence  ellipse  decreases  slightly  despite  the  smaller  number 
of  stations  used  in  the  location.  However,  the  actual  shot 
point  still  falls  outside  the  95%  confidence  ellipse  for  the 
computed  epicenter.  With  anomaly  corrections,  these  locations 
are  essentially  the  same  as  those  obtained  using  the  whole 
network. 

Figure  oC  shows  that  there  are  considerable  differences 
m  travel-times  between  stations  in  the  eastern  and  western 
United  States  with  respect  to  the  Herrin  68  travel-time  tables. 
The  eastern  stations  beyond  1000  km  recorded  early  arrivals, 
whereas  the. western  stations  beyond  2000  km  are  systematically 
late  by  about  two  seconds.  Similar  variations  in  travel-times 
were  reported  by  Mansfield  and  Evernden  (1966)  for  the  1963 
Lake  Superior  Seismic  Experiment.  To  investigate  the  effect 
of  these  regional  differences  in  travel-times  on  the  computed 
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epicenters,  we  divided  the  LRSM  and  observatory  network  into 
two  groups  of  stations.  One  group  contained  all  stations  east 
of  the  Rocky  Mountains  (Figure  28);  the  other  contained  all 
stations  within  and  to  the  west  of  the  Rocky  Mountains.  The 
locations  computed  using  these  two  networks  with  no  anomaly 
corrections  are  shown  in  Figure  31  and  Table  8. 

The  results  for  the  eastern  network  alone  show  consider¬ 
able  improvement  over  those  obtained  from  the  whole  network. 

The  mean  displacement  of  the  computed  epicenters  is  only  8.5  km. 
The  mean  area  of  confidence  ellipses  is  somewhat  larger,  prob¬ 
ably  because  of  the  smaller  number  of  stations  used  in  computing 
the  locations.  In  all  cases,  the  actual  shot  point  lies  with¬ 
in  the  95%  confidence  ellipse  for  the  computed  epicenter.  The 
results  for  the  western  stations  show  the  computed  epicenters 
considerably  displaced  to  the  southwest.  This  is  not  what 
one  vould  expect  from  a  consideration  of  the  observed  travel- 
times,  which  are  late  with  respect  to  Herrin  68  times.  How¬ 
ever,  since  all  stations  are  located  in  one  quadrant,  the  pro¬ 
gram  adjusted  the  origin  time  until  the  slope  of  the  travel¬ 
time  curve  fitted  the  slope  of  the  theoretical  travel-time 
curve,  with  the  result  that  the  estimated  origin  times  are 
late  by  some  9  seconds;  thus  allowing  the  computed  epicenters 

to  move  southwest.  The  95%  confidence  ellipses  have  a  mean 

2 

area  of  approximately  12,000  km  and  are  elongated  in  a 
NE-SW  direction  with  axial  ratios  of  around  4:1.  This  re¬ 
flects  the  narrow  distance  and  azimuth  aperture  of  the  network 
(Flinn,  1965). 

Spence. and  Alexander  (1968,  1969)  have  developed  an  al¬ 
ternative  method  for  relative  event  location  based  on  a  scheme 
used  by  Gutenberg  (Richter,  1958;  Appendix  X)  to  improve  the 
location  of  shallow  teleseisms.  We  define  the  travel-time 
difference  to  a  given  station  from  two  seismic  events  as: 
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when  is  the  observed  travel-time  from  the  reference  and  T? 
is  the  travel-time  from  the  event  whose  location  we  wish  to 
determine.  Then  for  two  surface  events  (Figure  32)  we  may  ex¬ 
press  t*  at  a  distant  station  as, 

.  •  sin  i  t  ~  sin  i  „  ,  ,  . 

t  -  _ o  6*r  _  _ _o.  5  cos  (ct-<J>) 

v^  v 

o  o 

where  i  is  the  take-off  angle  at  the  source,  v  the  P— wave 

_  o 

velocity  at  the  source,  6  is  the  distance  vector  between  ref¬ 
erence  event  and  the  unknown  event,  whose  length  6  and  azimuth 
a  we  wish  to  determine,  and  r  is  the  unit  vector  along  the 

azimuth  from  the  reference  event  to  the  station.  But  s^n  ^o 

vq 

is  simply  p(A),  the  seismic  ray  parameter.  Thus  by  plotting 

d  =pfA)vs  ^  we  can  estimate  6  and  a  from  the  amplitude  and 
intercepts  of  the  resulting  waves  (Figure  33). 

Selecting  shot  49A  as  the  reference  event,  we  formed  the 
quantity  t*  by  an  approximation  to  p(A)  given  by  the  first 
difference  of  the  Herrin  68  travel-time  tables  about  A  for 
each  station,  to  form  d (0 ) .  The  values  of  d  for  each  shot- 
pair  were  then  fitted  by  least-squares  to  a  smooth  curve  from 
which  we  calculated  6  and  cx»  Figure  34  shows  the  relative  lo¬ 
cations  of  the  1964  series  of  shots  computed  by  this  method. 

The  mean  location  error  (see  Table  9)  is  2.0  km.  Both  the 
locations  and  the  local  ion  errors  are  similar  to  those  com¬ 
puted  from  the  travel-time  anomalies. 

The  results  of  locating  the  ten  shots  using  program 
oHIFT,  show  that  utilizing  the  whole  LRSM  and  observatory  net¬ 
work,  we  are  able  to  locate  events  in  Lake  Superior  to  within 
15.6  km  of  the  actual  locations.  The  results  also  show  that 
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for  the  whole  network  there  is  an  apparent  "location  bias," 
since  the  computed  epicenters  lie  on  almost  a  common  azimuth 
f^om  the  actual  shot  points.  This  may  be  due  to  the  geometry 
of  the  array  since  there  are  only  two  stations  between  325° 
and  100°  in  azimuth  from  the  shot  points.  Removing  the  west¬ 
ern  stations  from  the  LRSM  and  observatory  network  improves 
the  location  in  that  it  reduces  the  distance  of  the  computed 
epicenter  from  the  shot  point.  However,  the  area  of  the  con¬ 
fidence  ellipse  for  the  computed  epicenter  increases.  This 
is  due  primarily  to  the  smaller  number  of  stations  in  the 
eastern  network  since  the  azimuthal  aperture  of  the  network 
remains  fixed.  The  accuracy  and  precision  of  the  locations 
from  the  western  network  of  stations  is  poor.  This  is 
largely  due  to  the  narrow  azimuthal  aperture  (29°). 

Using  travel-time  anomalies  obtained  from  one  of  the 
shots,  we  were  able  to  locate  (with  Program  SHIFT)  the  re¬ 
mainder  to  within  2.1  km  with  a  high  degree  of  confidence. 

The  relative  locations  obtained  using  the  Spence  and  / lexander 
approach  were  quite  comparable  in  accuracy  (average  2.0  km) 
which  is  expected  because  both  methods  account  for  variations 
in  propagation  velocity  to  individual  stations.  The  2  km 
accuracy,  therefore,  probably  represents  the  practical  limit 
on  location  of  events  at  Lake  Superior  using  this  network. 
Moreover,  in  view  of  the  known  structural  complexity  of  Lake 
Superior  the  travel-time  anomalies  derived  using  shot  49A  as 
the  master  event  probal ly  cannot  be  used  to  locate  events  to 
this  accuracy  if  they  are  more  than  about  25  km  from  shot  49A. 

Travel-times  and  Upper  Mantle  Structures  from  Events  in  Lake 
Superior 

Figure  2  shows  the  configuration  of  the  seismic  refraction 
profiles  which  were  used  to  record  the  Early  Rise  shots. 


Warren  et  al,  (1967)  and  Iyer  et  al,  (1969)  have  published 

xininary  "travel— times  and  amplitudes  of  the  first  arrivals 
recorded  during  Project  Early  Rise.  In  addition,  Warren  et  al, 
(1967)  show  record  sections  for  many  of  the  profiles.  Some 
detailed  analyses  of  this  data  have  been  published  by  Lewis 
and  Meyer  (1968),  Green  and  Hales  (1967a, b;  1968),  Barr  (1967), 
Mereu  and  Hunter  (1969)  and  Willis  (1968).  Roller  and  Jackson 
(1966a, b)  have  reported  on  travel-times  for  the  October  1964 
series  of  shots  in  Lake  Superior.  Iyer  et  al.  (1969)  pub¬ 
lished  an  analysis  of  all  the  first-arrival  data  from  Project 
Ear xy  Rise  and  derived  a  P-wave  velocity  model  from  the  com¬ 
bined  data.  After  extensive  smoothing,  they  inverted  the  re¬ 
sulting  travel— time  data  by  the  usual  Wiechert— Herglotz  method 
( Bullen ,  1963).  By  treating  the  data  in  this  manner  they 
effectively  averaged  out  any  regional  differences  in  travel- 
times.  oince  it  is  these  regional  differences  in  travel- 
times  and  their  implications  which  we  wish  to  examine,  we  will 
simply  show  the  reduced  travel- times  from  Project  Early  Rise 
(Figure  35)  together  with  their  theoretical  travel-times 
calculated  from  the  model  NCER  1  shown  in  Figure  36. 

Lewis  and  Meyer  (1968)  investigated  the  upper  mantle  to 
the  west  of  Lake  Superior  using  P  and  S  wave  arrivals  recorded 
along  a  profile  extending  from  400  km  to  2500  km  due  west  of 
i  “•-*  Lake  (Profile  1  in  Figure  2).  Figure  37  shows  a  record 
section  for  Profile  1  showing  the  P-wave  arrivals  compared 
with  the  arrivals  computed  from  the  model  obtained  by  Lewis 
and  Meyer  (Figure  39),  Figure  38  is  a  similar  figure  showing 
the  S-wave  arrivals  compared  with  the  theoretical  travel- 
times  computed  from  their  model.  An  interesting  feature  of 
Lewis  and  Meyer's  model  (Figure  39)  is  the  rapid  increase  in 
P-wave  velocity  at  a  depth  of  126  km  which  is  accompanied  by 
a  sharp  decrease  in  shear  wave  velocity.  At  150  km  depth  there 
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is  a  zone  40  km  thick  in  which  the  P-wave  velocity  decreases 
by  0.04  km/sec.  This  is  the  region  in  which  the  shear  waves 
show  a  pronounced  low  velocity  zone.  Both  the  P-wave  and 
S-wave  velocities  increase  abruptly  at  450km  depth.  Lewis  and 
Meyer  (1968)  also  calculated  two  distributions  of  density  with 
depth  from  the  two  relationship  between  P-wave  velocity  and 
density  given  by  Birch  (1964). 

To  the  southwest  of  Lake  Superior,  Green  and  Hales 
(1967a, b;  1968)  investigated  travel-times  along  two  radial  pro¬ 
files  (Profiles  2  and  3  in  Figure  2)  and  an  arcuate  profile 
(Profile  4).  Figure  40  shows  the  travel-times  of  the  first 
arrivals  and  the  theoretical  travel-time  curves  computed  from 
the  two  models  shown  in  the  insert.  The  travel-times  from 
Early  Rise  Model  1  fit  the  observed  data  to  within  0.8  sec. 
Early  Rise  Model  2  is  an  alternative  model  with  a  low  velocity 
zone  between  130  and  160  km  depth.  The  theoretical  travel- 
times  from  this  model,  shown  by  a  dashed  line  in  Figure  38, 
also  fit  the  observations  well,  accounting  for  the  offset  in 
the  observed  times  at  1500  km.  However,  Green  and  Hales 
(1968)  emphasize  that  this  offset  could  be  accounted  for  by 
reading  error  or  lateral  variations  (in  crustal  thickness  or 
mean  crustal  velocity),  and  prefer  Early  Rise  Model  1,  con¬ 
cluding  that  the  low  velocity  zone  for  P-waves  is  either  very 
thin  or  occurs  at  a  depth  greater  than  150  Jem  beneath  the 
central  United  States.  The  travel-times  along  a  profile  north¬ 
east  from  the  GNOME  explosion  near  Carlsbad,  New  Mexico  (Romney 
et  al,  1962;  Lehmann,  1964),  which  essentially  reverses  Pro¬ 
file  2,  agree  well  with  the  travel-times  calculated  from 
Early  Rise  Model  1. 

Profile  5  in  Figure  2  very  nearly  reverses  an  earlier  pro¬ 
file  occupied  by  the  University  of  Michigan  to  record  shots 
fired  in  the  Atlantic  during  the  1962  North  Carolina  Experiment 


(Meyer  et  al,  1965),  the  CHASE  series  of  underwater  explosions 
-r.  the  Atlantic,  and  the  East  Coast  On-Shore  Off-Shore  Experi¬ 
ment  (Hales  et  al;  1967).  Figure  41  is  a  plot  of  the  travel 
times  of  the  first  arrivals  from  the  Early  Rise  shots  recorded 
al°ng  Profile  5.  Willis  (1968)  has  inverted  this  data  to  ob¬ 
tain  the  following  P-wave  velocity  model  for  the  Lake  Superior 
end  of  the  profile. 


Velocity  km/sec 

Depth  km 

Thickness  km 
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8.  3 

50.4 
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8.4 

72.0 

OO 

the  southern  end  of 

the  profile. 

Willis  (1968)  found  a 

^8  km  thick  crust  at  the  continental  margin  and  the  higher 
velocity  layer  (8.4  km/sec)  at  a  depth  of  99  km. 

Also  to  the  southwest  of  Lake  Superior,  Roller  and  Jackson 
(1967  a,  b)  analyzed  the  travel-times  from  the  October  1964 
series  of  shots  along  a  profile  extending  from  the  Bayfield 
Peninsula  to  Denver,  Colorado,  which  was  extended  to  central 
Arizona  by  incorporating  several  of  the  LRSM  stations.  Figure 
42  shows  Roller  and  Jackson’s  results.  Roller  and  Jackson 
(1967a, b)  fit  two  linear  segments  to  the  mantle  first-arrival 
data  representing  waves  with  apparent  velocities  of  8.2  km/sec 
and  8. 5  km/sec,  the  crossover  occurring  at  a  distance  of 
1100  km  from  the  shot  point. 

Barr  (1967)  studied  the  travel-times  for  the  Early  Rise 
shots  and  the  Hudson  Bay  se:.smic  experiments  in  196  5  recorded 
at  the  standard  seismograph  stations  maintained  across  Canada 
by  the  Dominion  Observatory.  The  P-wave  trav  -times  were 
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referenced  to  a  standard  crust  by  subtracting  the  difference 
between  the  time-term  for  the  Early  Rise  shot  point  (4.15  sec 
obtained  from  Berry  and  West’s  (1966)  data)  and  the  mean  time- 
term  from  the  Hudson  Bay  shots  (3.56  sec.  derived  from  th f. 
data  obtained  by  Hobson  et  al,  (1967))  from  the  Early  Rise 
travel-times.  This  reduced  the  Early  Rise  times  by  0.59  s;c. 

The  corrected  travel-times  from  the  Lake  Superior  shots  to¬ 
gether  with  those  from  the  Hudson  Bay  experiment  are  shown  in 
Figure  43.  Figure  43  also  shows  the  theoretical  travel-times 
calculated  from  Barr’s  model.  In  constructing  this  model, 

Barr  (1967),  determined  the  positions  of  the  cusps  in  the  ob¬ 
served  travel-times  to  ensure  that  no  low  velocity  layers 
were  omitted  from  his  model.  Barr  (1967)  concludes  that  no 
low  velocity  zones  for  P-waves  in  the  upper  (350  km  of  the) 
mantle  are  indicated  by  his  data. 

The  8.42  km/ sec  layer  in  Barr’s  model  occurs  at  126  km 
depth.  However,  the  first  arrivals  at  three  of  his  stations 
are : 

BLC  T  =  8.27  +  at  310-960  km 

FFC  T  =  7.27  +  at  720-1390  km 

FBC  T  =  10.41  +  at  790-1340  km 

When  combined,  these  irrivals  yield 

T  =  6.74  +  J^"22  310-1390  km  distance  range 

One  can  readily  see  that  two  of  the  three  independent  station 
velocities,  BLC  and  FBC,  are  significantly  higher  than  the  vel¬ 
ocity  from  the  combined  data.  Barr  (1967)  attributes  these 
higher  velocities  to  local  structure  in  the  upper  mantle  beneath 
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.tie  Moho  under  the  Hudson  Bay,  since  the  immediate  sub- 
crustal  velocity  is  known  from  earlier  work  (Hobson  et  al, 
1967;  Hall  and  Brisbin,  1961;  Hodgson,  1953)  to  be  very 
nearly  8.2  km/sec. 

..ereu  and  Hunter  (1969)  have  examined  the  travel- 
times  to  stations  along  a  profile  north  from  Lake  Superior  to 
Churchill ,  Manitoba  (Profile  8  in  Figure  2).  The  record  sec¬ 
tion  obtained  by  Mereu  and  Hunter  (Figure  44)  shows  three 
distinct  branches  of  the  travel— time  curve  with  apparent  vel¬ 
ocities  of  (a)  6.3  km/sec  up  to  2.4°,  (  b)  8.10  km/sec  from 
,*4°  to  5.8C,  (c)  8.54  km/sec  from  5.8°  to  13.4°.  The 

branch  of  the  curve  with  apparent  velocity  of  8.54  km/sec 
corresponds  to  arrivals  refracted  from  a  layer  with  velocity 
of  8.43  km/sec  which,  when  corrected  for  earth  curvature, 
occurs  at  a  depth  of  84  +  3  km  (Mereu  and  Hunter.,  1969). 

±iic  possibility  of  a  low  velocity  zone  in  the  upper  mantle 
for  F-waves  cannot  be  ruled  out  on  the  basis  of  Mereu  and 
Hunter's  data  since  first  arrivals  in  the  distance  range 
1140  to  1350  km  have  very  small  amplitudes.  There  is  also 
some  suggestion  of  an  offset  in  the  travel— time  curve  occur¬ 
ring  in  this  distance  range. 

Brune  and  Dorman  (1963)  in  their  study  of  surface  wave 
dispersion  across  the  Canadian  Shield  used  the  travel-times 
reported  by  Hodgson  (1953a, b),  Lehmann  (1955),  and  Hall  and 
Brisbain  (1961)  to  supplement  their  surface-wave  dispersion 
data  in  determining  the  model  CANSD  (Figure  2).  Brune  and 
Dorman  (1963)  report  a  Pn  velocity  of  8.18  km/sec  at  about 
400  km  distance  increasing  to  8.33  km/sec  at  about  2000  km, 
and  a  Sn  velocity  of  4.72  km/ sec  over  this  range.  The  com- 
pressional  wave  velocities  in  the  lower  layers  of  CANSD  were 
assigned  rather  than  observed,. 
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PHYSICAL  SETTING  OF  THE  NEVADA  TESTS 


Geological  Environment 

The  Relationship  of  the  Nevada  Test  Site  to  the 

Basin  and  Range 

The  Basin  and  Range  province  covers  a  large  area  of  the 
south-western  United  States  (Figure  45).  It  has  a  character¬ 
istic  landscape  consisting  of  short,  sub-parallel  mountain 
ranges  and  intermontane  desert  basins  (King,  1959,  p.152). 

Much  of  the  area  (the  Great  Basin  of  King)  is  a  region  of  in¬ 
terior  drainage.  The  topography  of  the  Basin  and  Range  pro¬ 
vince  is  structurally  controlled,  consisting  of  a  series  of 
horst,  grabens  and  tilted  blocks  which  are  bounded  by  Cenozoic 
normal  faults  (Eardley,  1962;  Gilluly,  1963;  Thompson,  1966; 
Thompson  et  al,  1967)  and  has  been  evolving  since  the  time  of 
the  Cordilleran  orogeny. 

According  to  King  (1957),  during  early  Tertiary  time,  in 
the  eastern  section  of  the  Great  Basin  deposition  of  lacustrine 
and  fluvial  sediments  took  place  over  a  broad  flood  plain 
bounded  to  the  east  by  the  low  lying  hills  of  the  Rocky  Moun¬ 
tain  region  and  the  denuded  Sierra  Nevada  Batholith  to  the  west. 
During  the  middle  Tertiary  time,  the  western  part  of  the  Great 
Basin  was  covered  by  volcanic  rooks,  and  in  some  areas  was 
faulted  and  tilted  during  the  middle  Miocene.  In  the  late 
Miocene  and  early  Pliocene  deposition  of  terresirial  sediments 
took  place  and  fossil  plants  occurring  in  these  rocks  suggest 
that  the  whole  area  had  ->ly  moderate  relief  and  was  relatively 
low-lying.  Uplift  of  the  Sierra  Nevada  block  began  in  late 
Pliocene  and  continued  into  the  Pleistocene.  During  this  time 
the  Great  Basin  area  was  disrupted  by  extensive  block  faulting, 
wtiich  resulted  in  the  uplift  of  not  only  the  present  day  mountains 


but  the  basins  also.  Geophysical  evidence  suggests  that  this 
uplift  is  still  taking  place  today. 

The  Basin  and  Range  province  is  a  tectonically  active 
region.  Seismological  studies  have  shown  that  the  structure 
of  the  crust  and  particularly  the  upper  mantle  is  quite  differ¬ 
ent  from  that  found  in  non-tectonic  regions.  Hill  and  Pakiser 
(1966)  in  a  study  of  a  series  of  reversed  seismic  refraction 
profiles  between  Eureka,  Nevada  and  Boise,  Idaho  concluded  that 
the  crust  in  the  northern  Basin  and  Range  province  consists  of 
a  6.0  km/sec  layer  19  to  24  km  thick  overlying  a  6.7  km/sec 
layer  some  10  to  12  km  thick  and  that  both  the  crust  as  a 
whole  and  the  6.7  km/sec  layer  increase  in  thickness  at  the 
boundary  between  the  Basin  and  Range  province  and  the  western 
Snake  River  Plain.  They  also  found  a  Pn  velocity  of  7.9  km/sec 
obtained  by  Pakiser  (1963),  Stuart  et  al,  (1964),  Pakiser  and 
Steinhart  (1964),  and  Healy  and  Ploger  (1965).  In  addition  to 
the  northwards  thickening,  Hill  and  Pakiser  find  that  the  Moho 
shows  local  variations  in  depth  of  as  much  as  5  km  under  the 
Basin  and  Range  province,  which  are  probably  due  to  faulting. 

Hill  and  Pakiser's  (1966)  interpretation  of  crustal  struc¬ 
ture  in  the  Northern  Basin  and  Range  has  recently  been  verified 
by  the  work  of  Prodehl  (1969)  Prodehl  reinterpreted  the  data 
accumulated  by  the  U.  S.  Geological  Survey  during  the  years 
1961-1964.  He  found  that  the  mean  crustal  thickness  in  the 
northern  Basin  and  Range  province  is  30  km,  and  the  average 
P-wave  velocity  is  6.1  to  6.2  km/ sec  depths  of  15  or  20  km, 
there  being  good  evidence  for  an  intermediate  layer  in  the 
northern  but  not  in  the  southern  part  of  the  Great  Basin. 
Alexander  (1963)  also  found  evidence  for  an  intermediate  layer 
from  observations  of  surface  wave  dispersion  across  the  Basin 
and  Range  province,  although  the  intermediate  layer  in  his 
model  is  somewhat  deeper  and  has  a  higher  velocity. 


The  physiographic  delimitation  of  the  Basin  and  Range  is 
relected  in  the  spatial  distribution  of  upper  mantle  provinces 
in  the  western  United  States.  Archambeau  et  al.  (1969)  have 
demonstrated  that  a  correspondence  exists  between  patterns  of 
compressional  wave  velocity  at  the  top  of  the  mantle  (P  veloc¬ 
ity)  P-wave  travel-time  delays  and  the  various  physiographic 
provinces;  the  Basin  and  Range  province  being  characterized 
by  anomalously  low  P^  velocities  and  positive  P-wave  delays 
with  respect  to  the  Jef f reys-Bullen  travel-time  tables.  Its 
boundaries  are  defined  by  rapid  changes  in  Pn  velocity.  Their 

results  also  show  that  variations  in  P  velocities  are  mani- 

n 

festations  of  profound  differences  in  the  structure  of  the  upper 
mantle.  Figure  46  shows  Archambeau  et  al.'s  model  CIT  111P 
for  P-wave  velocity  structure  beneath  the  Basin  and  Range 
piovince.  The  model  shows  no  velocity  reversals  with  depth 
but  rather  the  low  velocity  zone  (7.7  km/sec)  for  P-waves  ex¬ 
tends  from  a  depth  of  150  km  to  the  base  of  the  crust.  This 
latter  zone,  on  the  basis  of  the  results  of  Archambeau  et  al. 
(1968)  is  where  major  differences  in  mantle  structure  beneath 
the  various  physiographic  provinces  of  the  western  United 
States  probably  occur. 

Shear  wave  velocity  models  also  indicate  the  presence  of 
a  low  velocity  zone  in  the  mantle.  Kovach  (1968),  and  Kovach 
and  Robinson  (1969)  have  studied  dt/dA  for  long  period  S-waves 
recorded  at  the  extended  TFO  array  in  Arizona.  The  important 
features  of  their  model  which  is  shown  in  Figure  47  are  (1)  a 
thin  lid  zone,  9  km  thick  with  a  shear  wave  velocity  of  4.5 
km/sec  overlying  a  low  velocity  zone;  (2)  a  change  in  velocity 
gradient  at  a  depth  of  160  km  followed  by  a  gradual  increase 
in  velocity  to  a  depth  of  360  km;  (3)  a  rapid  increase  in 
velocity  beginning  at  a  depth  of  620  km.  The  main  features 
of  Kovach  and  Robinson’s  model  are  consistent  with  the 
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previous  results  of  Toksoz  and  Anderson  (1965)  from  studies 
of  Love  wave  dispersion  in  tectonic  regions. 

Heat  flow  measurements  in  the  Basin  and  Range  province 

reported  by  Roy  et  al.  (1968a)  indicate  a  mean  heat  flow  of 
2 

2.0  ycal/cm  sec  (Blackwell,  1969).  More  recent  results  re¬ 
ported  by  Warren  et  al.  (1969)  suggest  that  the  mean  heat  flow 

r\ 

may  be  as  high  as  2.21  +  0.43  ycal/cm  -sec  in  the  Basin  and 
Range  province  south  of  36°N.  Both  of  these  values  are  sig¬ 
nificantly  higher  than  the  world-wide  mean  heat  flow  of  1.5 
ycal/cm  -sec  (Lee  and-~Uyeda,  1965  ).  Figure  48  summarizes 
the  heat  flow  measurements  made  in  the  western  United  States 
to  date.  The  high  heat  flow  over  the  Basin  and  Range  is 
clearly  delimited.  To  the  west,  the  mean  heat  flow  in  the 
Sierra  Nevada  is  somewhat  low  (  *0.9  ycal/cm2  sec,  mean  of 
7  values  published  by  Roy  et  al,  1968c)  ,  and  to  the  east 
in  the  Colorado  Plateau  and  Texas  Foreland,  the  heat  flow  is 
also  somewhat  low  (*1.1  -  1.2  ycal/cm2  sec,  Warren  et  al,  1969). 
Roy  and  Blackwell  (1966)  investigated  the  heat  flow  transition 
between  the  Basin  and  Range  province  and  the  Sierra  Nevada  and 
found  that  the  heat  flow  and  seismic  boundaries  are  coincident 
and  about  50  km  east  of  the  structural  and  physiographic 
boundary  (see  Figure  49).  When  the  contribution  to  the  surface 
heat  flow  from  sources  within  the  upper  10  km  of  the  crust  is 
removed  using  the  method  of  Roy  et  al,  (1968b, c),  the  contrast 
is  even  sharper.  This  suggests  that  the  broad  heat  flow  anom¬ 
aly  over  the  Basin  and  Range  province  is  not  related  to  excess 
heat  production  in  the  upper  crust,  but  rather  that  it  is 
due  to  anomalously  high  heat  flow  in  the  upper  mantle.  Warren 
et  al ,  (19  69  )  examined  the  heat  flow  in  the  southeastern  part 
of  the  Basin  and  Range  and  the  Texas  Foreland.  Again,  they 
found  a  sharp  transition  between  the  high(>2 . Oycal/cm2  sec) 
heat  flow  in  the  Basin  and  Range  and  the  normal  or  somewhat 
low  (*1.0  ycal/cm  sec)  heat  flow  in  the  Texas  Foreland 
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which  takes  place  at  150  km.  Warren  et  al,  (1969)  interpret 
the  high  heat  flow  in  terms  of  high  subcrustal  temperatures. 
This  interpretation  is  consistent  with  the  geomagnetic  data 
of  Schmuker  (1964)  and  the  magnetotelluric  data  of  Swift  (1967) 
which  indicate  a  zone  of  anomalously  high  electrical  conductiv¬ 
ity  under  southern  Arizona  and  New  Mexico.  Both  Schmucker 
(1964)  and  Swift  (1967)  have  interpreted  this  zone  to  mean  that 
the  higher  temperature  isotherms  are  shallower  in  the  Basin 
and  Range  than  in  the  Texas  Foreland. 

The  coexistence  of  high  heat  flow,  low  seismic  velocities 
and  the  low  values  of  Q  in  the  upper  mantle  (Archambeau  et  al, 
1969)  suggests  that  the  low  velocity  zone  may  be  due  to  partial 
melting  of  the  mantle  material.  The  absence  of  a  high  velocity 
lid  to  the  low  velocity  zone  for  P-waves  in  the  Basin  and  Range 
province  implies  that  the  zone  of  partial  melting  extends  to 
the  base  of  the  crust.  Archambeau  and  Davies  (1969)  have  ex¬ 
amined  the  geophysical  data  and  concluded  that  they  are  con¬ 
sistent  wiu’h  the  hypothesis  that  the  Basin  and  Range  province 
overlies  the  extension  of  the  East  Pacific  Rise  (Heezen,  1960; 
Menard,  1960;  Cook,  1966;  Thompson,  1966)  where  the  North 
American  crustal  plate  has  overridden  the  spreading  zone.  If 
this  hypothesis  is  substantially  correct,  then  we  should  ex¬ 
pect  observations  of  travel-times,  amplitudes,  and  the  asso¬ 
ciated  models  of  crustal  and  mantle  structure  to  appear  anoma¬ 
lous  when  compared  to  the  tectonically  inactive  or  stable  re¬ 
gions  of  North  America. 

In  the  following  sections  we  summarize  some  of  the  ob¬ 
servations  which  demonstrate  the  anomalous  nature  of  the  crust 
and  mantle  beneath  Basin  and  Range  province  in  much  the  same 
manner  as  we  did  above  for  Lake  Superior.  In  particular,  we 
shall  be  concerned  with  the  derivation  of  a  geophysical  model 
for  NTS  which  will  serve  as  a  basis  for  comparison  between  the 
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two  regions. 


The  Nevada  Test  Site  (NTS)  is  situated  in  the  Great  Basin 
section  of  the  Basin  and  Range  province  (Figure  50)  approx¬ 
imately  110  km  to  the  northwest  of  Las  Vegas.  Some  recent 
tests,  for  example  FAULTLESS,  have  been  conducted  somewhat  to 
the  north  of  the  main  site,  as  were  certain  earlier  events 
such  as  SHOAL.  However,  we  will  discuss,  principally,  the 
main  test  area  shown  in  Figure  50,  and  will  attempt  to  derive 
a  geophysical  model  for  the  Yucca  Flat  area. 

Geology  of  the  Nevada  Test  Site 

The  geology  of  the  Yucca  Flat  and  Frenchman  Flat  area  of 
the  Nevada  Test  Site  (Figure  51)  is  fairly ’well  known  from 
surface  observations  and  borehole  data.  In  1952  Johnson  and 
Hibbard  conducted  a  reconnaissance  survey  of  the  area  the  re¬ 
sults  of  which  were  published  in  1957  (Johnson  and  Hibbard,  1957). 
Over  the  period  1960-63  the  U.  S.  Geological  Survey  remapped 
the  Yucca  Flat  area  and  conducted  numerous  subsurface  investi¬ 
gations.  The  results  are  published  in  a  series  of  Technical 
Letters  and  it  is  from  these,  together  with  Memoir  110  of  the 
Geological  Society  of  America  (Eckel,  1968),  that  much  of  the 
information  presented  below  is  drawn. 

Physiographically ,  Yucca  Flat  is  a  broad,  relatively  flat 
intermontane  basin  bordered  to  the  north,  east  and  west  by 
mountains  of  moderate  relief,  and  separated  from  Frenchman  Flat 
to  the  south  by  a  somewhat  lower  ridge  of  pre-Quaternary  rocks. 
The  basin  is  an  area  of  inland  drainage,  a  small  playa  (Yucca 
Lake)  forming  during  the  wet  season  at  the  southern  margin  of 
the  flat. 

In  general,  the  mountains  are  composed  of  late  Pre- 
cambrian  and  Paleozoic  sedimentary  rocks  and  Cenozoic  volcanics 
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whereas  the  valleys  are  filled  with  detritus  derived  from  the 
surrounding  mountains.  The  rocks  exposed  at  the  surface  in 
the  Yucca  Flat  area  range  in  age  from  late-Precambrian  to  Re¬ 
cent.  They  may  be  conveniently  divided  into  three  main  groups, 
the  late-Precambrian-Permian ,  the  Cenozoic  volcanics ,  and  the 
Quaternary  and  Recent  valley  fill.  The  stratigraphy  of  the 
area  is  summarized  in  Table  10.  The  first  group  consists  of 
some  10.6  km  of  argillites,  quartzites,  limestones  and  dolo¬ 
mites  which  are  divided  into  eight  sequences  on  the  basis  of 
gross  lithology  and  age  (Barnes  et  al,  1963).  The  second  group 
consists  of  volcanic  rocks  of  Miocene  or  younger  age  (radio- 
metric  ages  obtained  fnom  the  ash  flow  tuffs  range  from  26  my 
to  7  my  (Ekren ,  1968))  locally  more  than  1.06  km  thick 
(Barnes  et  al,  1963).  (The  above  thicknesses  refer  to  the 
stratigraphic  section  not  the  present  vertical  extent  of  these 
units.)  These  rocks  are  composed  of  tuffs,  welded  tuffs, 
tuffaceous  sedimentary  rocks  and  lava  flows  which  were  de¬ 
posited  on  the  deeply  eroded  surface  of  the  Palezoic  rocks. 
Seven  principal  ash-flow  tuffs  have  been  recognized  (Ekren, 
1968)  which  range  in  lateral  extent  from  100  to  200  km,  and 
are  thought  to  have  originated  to  the  west  of  Yucca  Flat.  In 
the  valleys  both  the  volcanics  and  the  pre-Cenozoic  rocks  are 
overlain  by  rocks  of  the  third  group,  the  Quaternary  and  Re¬ 
cent  valley  fill.  These  rocks  are  typical  fluvial  deposits 
composed  of  detritus  from  the  local  mountains,  ranging  from 
coarse  agglomerates  to  fine  silts  and  clays  with  locally  occur¬ 
ring  evaporites.  Together  these  deposits  are  referred  to  as 
alluvium.  The  alluvium  attains  a  maximum  thickness  of  0.61  km 
(Barnes  et  al,  1963)  beneath  the  south-central  portion  of 
Yucca  Flat. 

Although  the  predominant  structural  pattern  of  the  Yucca 
Flat  area  is  the  result  of  block  faulting,  detailed  mapping 
has  shown  that  the  area  is  structurally  very  complex.  Barnes 


et  al.  (1963)  divided  the  area  into  an  eastern  and  western  block, 
on  the  basis  of  the  age  of  the  outcropping  rocks  and  the  orien¬ 
tation  of  the  major  folds  and  faults  along  the  line  A-A*  in 
Figure  52.  The  rocks  in  the  eastern  block  are  chiefly  of  late 
Precambrian  and  early  Paleozoic  age.  Major  folds  trend  north¬ 
west,  and  the  high-angle  faults  strike  north-northeast.  The 
i  aleozoic  rocks  in  this  block  are  often  repeated  by  thrust 
faulting  (Figure  52),  which  in  some  areas  places  late  Precambrian 
and  lower  Paleozoic  rocks  over  middle  and  upper  Paleozoic  strata. 
The  Yucca  Fault,  which  runs  north-south  through  '.ucca  Flat  is 
mostly  within  the  eastern  block.  Quaternary  movement  along  this 
fault  has  caused  the  eastern  side  to  be  down-thrown  by  as  much 
as  0.3  km  (Barnes  et  al,  1963). 

In  the  western  block,  most  of  the  outcropping  pre-Cenozoic 
rocks  are  of  late  Paleozoic  age.  Here  the  rocks  are  folded 
along  northeast  trending  axes,  the  younger  Devonian  to  Permian 
rocks  being  overthrust  by  Cambrian  and  Ordovician  rocks.  The 
bigh-angle  faults  in  this  block  are  less  regular  in  trend  as 
can  be  seen  from  Figure  52.  They  are  predominently  oriented 
north  and  north-northeast,  with  apparent  vertical  displacements 
of  up  to  0.45  km.  Ekren  (1968)  attributes  the  change  in 
strike  from  north  to  northeast  in  the  northern  part  of  NTS  to 
bending  of  the  faulted  ranges  caused  by  drag  along  the  Las 
Vegas  Valley  shear  zone  (Longwell,  1960)  as  several  of  the 
northeast  trending  faults  show  left-lateral  displacements  of 
up  to  4  km. 

Yucca  Flat  itself  lies  between  the  two  blocks.  There  are 
no  outcrops  of  the  basement.  However,  the  structure  and  con¬ 
figuration  of  the  Paleozoic  and  Tertiary  rocks  are  fairly  well- 
known  from  seismic,  gravity,  magnetic,  and  bore-hole  data. 
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Geophysical  Evidence  on  Crustal  Structure  under  the  Nevada 


Test  Site 

Many  detailed  geophysical  investigations  (summarized 
in  Hazlewood  et  al^  1963)  have  been  made  of  the  Yucca  Flat 
area  in  order  to  elucidate  the  shallow  structure  of  the  area. 

In  addition,  seismic  surveys  have  been  carried  out  over  much 
of  Nevada  and  the  surrounding  states.  Many  of  these  seismic 
profiles  use  events  at  NTS  as  sources  of  seismic  energy  and 
so  it  is  not  possible  to  obtain  directly  the  deep  crustal 
structure  immediately  beneath  the  test  site  itself  from  these 
profiles ,  and  so  we  shall  have  to  study  the  seismic  data  on 
a  regional  basis.  In  addition  to  the  seismic  data,  local 
gravity  data  (Healey  and  Miller,  1963;  Healey,  1968)  and 
regional  Bouguer  gravity  data  (Joesting  and  Woollard,  1964) 
aeromagnetic  data  and  regional  heat  flow  data  (Blackwell, 

1969  )  are  available.  Furthermore,  NTS  lies  in  the  area 
covered  by  the  Transcontinental  Geophysical  Survey  (USGS, 1968). 

Seismic  Measurements 

Refraction  Studies 

Shallow  seismic  refraction  studies  (Hazelwood  et  al,  1963) 
were  made  in  the  Yucca  Flat  area  to  determine  the  thickness  of 
the  alluvium  beneath  the  flat  and  the  configuration  of  the 
Paleozoic  and  Tertiary  rocks  which  form  the  basement.  fable  10 
summarizes  the  seismic  velocities  for  the  various  rock  types 
outcropping  in  the  Yucca  Flat  region  as  reported  by  Hazelwood 
et  al,  (1963).  The  large  variation  in  velocity  for  alluvium 
is  due  to  the  rapid  increase  in  velocity  with  depth.  The  re¬ 
sults  of  a  seismic  refraction  survey  in  Yucca  Flat  aided  in 
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the  delineation  of  a  buried  Paleozoic  ridge,  which  from  gravity 
data  (discussed  below)  and  bore  holes  is  known  to  extend  north- 
south  through  the  center  of  the  flat. 

The  gross  crustal  structure  in  the  vicinity  of  the  Nevada 
Test  Site  is  known  from  the  U.  S.  Geological  Survey’s  crustal 
refraction  work  reported  in  Hill  and  Pakiser  (1966);  Pakiser 
(1963);  Eaton  (  1963  );  Diluent  et  al,  (1961);  Gibbs  and  Holler, 
(1964,  1966);  Johnson  (1965),  and  Prodehl  (1969).  In  addition, 
Pakiser  and  Hill  (1963),  and  Ryall  and  Stuart  (1963)  used  the 
travel-times  of  seismic  phases  generated  by  events  at  NTS  to 
study  crustal  structure  in  the  western  United  States.  We  have 
plotted  the  results  of  these  surveys  as  a  fence  diagram  (Fig¬ 
ure  53). 


To  the  north  of  NTS,  Hill  and  Pakiser  (1966)  studied 
travel-times  along  a  reversed  profile  to  Boise,  Idaho.  Their 
results  indicate  that  the  crust  is  30  km  thick,  100  km  north 
of  NTS,  increasing  to  48  km  under  the  western  Snake  River 
Plain  680  km  to  the  north  of  the  test  site.  In  the  vicinity 
of  NTS,  Hill  and  Pakiser  (1966)  (Figure  17,  p.  412),  show 
the  crust  composed  of  approximately  1  km  of  low-velocity 
"sediments”  overlying  a  19  km  thick  layer  with  velocity  of 
6.0  km/sec,  a  10  km  thick  intermediate  layer  with  a  velocity 
of  6.7  km/sec,  and  a  Pn  velocity  of  7.9  km/sec.  Their  model 
indicates  considerable  local  variations  in  crustal  thickness 
which  Hill  and  Pakiser  (1966)  attribute  to  faulting  at  the 
base  of  the  crust.  The  model  also  shows  that  both  the  crust 
as  a  whole  and  the  6.7  km/ sec  layer  rapidly  increase  in  thick¬ 
ness  at  the  physiographic  boundary  between  the  northern  Basin 
and  Range  Province  and  the  western  Snake  River  Plain.  The  P 
velocity  found  by  Hill  and  Pakiser  (1966)  agrees  with  earlier 
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studies  of  similar  data  along  the  sad  orofile  by  Pakiser  and 
Hill  (  1963  )  and  Knopoff  and  Teng  (190..  The  crustal  thick¬ 
ness  of  30  km  at  a  point  100  km  to  the  north  of  NTS  agrees  with 
the  results  obtained  by  Diment  et  al3  (1961)  from  analyses 
of  travel-times  from  events  at  NTS  recorded  along  a  300  km 
line  southeast  from  NTS  to  Kingman,  Arizona. 

To  the  east  of  NTS,  Ryall  and  Stuart  (1963),  have  inter¬ 
preted  travel-time  data  recorded  along  a  profile  from  NTS  to 
Ordway,  Colorado.  Their  results  indicate  a  two  layered  crust 
with  velocities  of  6.0  km/sec  and  6.5  km/sec,  which  increases 
in  thickness  from  26  km  at  a  point  approximately  60  km  east 
of  NTS  to  42  km  thick  beneath  the  western  part  of  the  Colorado 
Plateau  400  km  distant  from  NTS.  Ryall  and  Stuart  (1965) 
found  an  apparent  Pn  velocity  of  7.6  km/sec.  The  2°  35'  E  dip 
in  the  Moho  is  calculated  on  the  assumption  that  the  true  P 
velocity  at  NTS  is  7.9  km/sec.  The  thickness  of  the  interme¬ 
diate  layer  in  this  model  is  only  3.8  km.  Therefore  it  could 
not  possibly  be  recorded  as  a  first  arrival;  the  main  evidence 
for  its  existence  comes  from  the  large  amplitude  second  arrivals 
recorded  in  the  distance  range  123  -  319  km. 

To  the  south  of  NTS,  Gibbs  and  Roller  (1964,  1966)  have 

determined  the  crustal  structure  between  NTS  and  Ludlow, 

California  from  seismograms  recorded  along  a  reversed  profile 

between  the  two  points.  For  events  at  NTS,  they  found  that  the 

first  arrivals  out  to  155  km  had  a  velocity  of  6.10  km/sec 

which  they  correlate  with  Pg.  In  the  distance  range  155  km 

to  265  km  south  of  NTS,  the  first  arrivals  had  a  velocity  of 

8.04  km/sec  (P  ).  South  of  Ludlow,  in  the  distance  range 

315  km  to  490  km,  the  first  arrivals  had  a  velocity  of  7.75 

km/sec  (P  ).  For  events  at  Ludlow  recorded  northwards  from 

n 

Ludlow  toward  NTS,  they  found  first  arrivals  with  velocities 
of  6.10  km/sec  (Pg),  6.80  km/sec  (P*),  and  7.76  km/sec  (PR). 
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Gibbs  and  Roller  (1964,  1966)  internret  the  high  apparent  vel¬ 
ocity  of  recorded  south  from  NTS  as  evidence  that  the  crust 
thins  southwards  and  that  the  velocity  of  7.7  5  km/ sec  re¬ 
corded  south  of  Ludlow  indicates  that  the  Moho  flattens  out 
beyond  315  km  from  NTS,  Their  preferred  crustal  model  (Model 
II)  shows  a  total  crustal  thickness  of  24  km  at  Ludlow  and  34 
Km  at  NTS  made  up  of  the  following  layers: 
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the  thickness  of  the  6.10  km/sec  layer  being  assumed.  The  dip 
of  the  Moho  is  1.7°  S  (Gibbs  and  Roller  1964;1936). 

Johnson  (1965)  has  interpreted  the  crustal  structure  be¬ 
tween  Mono  Lake,  California  and  Lake  Mead,  Nevada,  His  pro¬ 
file  crosses  the  NTS  -  Ludlow  profile  (Gibbs  and  Roller,  1964; 
1966)  approximately  80  km  south  of  NTS.  Figure  53  shows  the 
structure  obtained  by  Johnson  (1965),  and  the  structure  ob¬ 
tained  by  Gibbs  and  Roller  at  the  point  where  the  two  profiles 
intersect.  The  total  crustal  thickness  of  33  km  found  by 
Johnson  agrees  well  with  the  32.5  km  obtained  by  projecting  the 
crustal  thickness  of  34  km  at  NTS  found  by  Gibbs  and  Roller, 
However,  the  agreement  could  well  be  fortuitous  since  Johnson’s 
model  consists  of  a  1.9  km  tnick  layer  with  velocity  3.0  km/sec 
(assumed),  a  2 3  km  thick  layer  with  velocity  6.15  km/ sec,  an 
intermediate  layer  8  km  thick  with  a  velocity  of  7.10  km/ sec 
and  a  P^  velocity  of  7.80  km/sec,  whereas  the  Gibbs  ^\nd  Roller 
model  has  an  intermediate  layer  20  km  thick  (see  above). 
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Johnson’s  model  is  also  interesting  in  that  it  shows  sharp 
flexures  of  6  km  in  the  Moho  25  km  south  of  Mono  Lake  and 
flexures  of  2* 5- km  in  the  Moho  30  km  north  of  Lake  Mead, 
which  correlate  with  large  surface  faults  (Johnson,  1965). 
Between  100  km  and  215  km  SW  of  Mono  Lake  the  crust  thins 
from  40  km  to  30  km  in  thickness.  This  latter  feature  (a 
monocline)  has  no  surface  expression;  however  Johnson  (1965) 
found  that  the  Bouguer  gravity  anomalies  calculated  from  the 
refraction  model  using  the  empirical  compressional  wave  veloc¬ 
ity  vs.  density  relations  of  J.  E.  Nafe  and  C.  L.  Drake  (Talwani 
et  alj  1959)  fitted  the  observed  Bouguer  gravity  anomalies 
along  the  profile. 

Northwest  of  NTS  Eaton  (1963)  interpreted  the  structure 
between  Fallon,  Nevada  and  Eureka,  Nevada.  Eaton  found  no 
wave  refracted  from  the  top  of  the  intermediate  layer  as  a 
first  arrival  on  this  profile;  the  main  evidence  for  the  exis¬ 
tence  of  the  intermediate  layer  comes  from  the  tentative 
identification  of  the  reflected  phase  PjP.  Assuming  that  an 
intermediate  layer  is  present,  Eaton  (1963)  found  that  the 
crust  at  Fallon  is  22.3  km  thick  and  consists  of  a  17.2  km 
thick  layer  with  an  adopted  velocity  of  6.0  km/sec,  and  a 
6.9  km  thick  layer  with  velocity  6.6  km/ sec.  Beneath  Eureka, 
Eaton  calculated  that  the  crust  is  33.9  km  thick,  consisting 
of  a  22.3  km  thick  6.0  km/sec  layer  and  a  11.6  km  thick  inter¬ 
mediate  layer  with  velocity  6.6  km/sec.  The  velocity, 
when  corrected  for  the  2°E  dip  in  the  Moho,  is  7.82  km/sec. 
Assuming  that  the  intermediate  layer  is  absent,  Eaton  calcu¬ 
lated  that  the  crust  (6.0  km/ sec)  is  22.2  km  thick  at  Fallon 
and  31.5  km  thick  at  Eureka,  the  true  velocity  being  7.82 
km/sec  when  . corrected  for  the  1.7°  dip  of  the  Moho.  The  32  - 
34  km  of  crustal  thickness  at  Eureka  found  by  Eaton  (1963) 
agrees  well  with  33  -  36  km  crustal  thickness  in  the  vicinity 
of  Eureka  reported  by  Hill  and  Pakiser  (1966,  p.  401). 
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Recently  Prodehl  (1969)  has  reinterpreted  the  data  from 
46  seismic-refraction  profiles  recorded  by  the  U.  S.  Geological 
Survey  in  the  western  United  States  during  the  years  1961  to 
1966.  He  concludes  that  the  average  crustal  thickness  under 
the  Basin  and  Range  province  is  30  km  and  the  mean  P-wave 
velocity  of  the  upper  15  to  20  km  of  the  crust  is  6.1  to  6.2 
km.  Prodehl  finds  good  evidence  for  an  intermediate  layer  in 
the  northern  part  of  the  Great  Basin  but  not  in  the  southern 
part.  However  he  states  that  the  curvature  of  the  branch  of 
the  travel-time  curves  corresponding  to  the  phase  P^P >  ..hich 
is  a  dominant  phase  on  all  profiles,  is  explained  by  a  rela¬ 
tively  large  velocity  gradient  at  the  base  of  the  crust.  Thus, 
Prodehl  (1969)  concludes  that  it  is  possible,  but  somewhat  un¬ 
likely,  that  this  zone  of  steep  velocity  gradients  is  10  -  15 
km  thick  beneath  the  southern  part  of  the  Great  Basin. 

Surface  Wave  Dispersion  Observations 

Alexander  (1963)  studied  Love  and  Rayleigh  wave  group 
velocity  dispersion  across  the  Basin  and  Range  province  from 
recordings  of  two  events  along  the  Utah-ldaho  border  and  the 
aftershock  sequence  of  the  Montana  earthquake  of  August  1959. 

His  data  together  with  his  model  35CM2  derived  for  the  crustal 
and  upper  mantle  structure  of  the  Basin  and  Range  province  are 
shown  in  Figures  54,  55,  56,  57.  An  unusual  feature  of  model 
35CM2  is  the  25  km  thick  layer  with  a  S-wave  velocity  of  4.1 
km/sec  and  a  P-wave  velocity  of  7.7  km/ sec  at  25  km  depth 
which  Alexander  (1963)  designates  section  X.  A  thin  high 
shear  velocity  layer  at  a  depth  of  about  50  km  is  necessary 
to  bring  the  theoretical  group  velocities  of  the  fundamental 
mode  Love  and  Rayleigh  waves  at  the  longer  periods  into  agree¬ 
ment  with  the  observations.  Alexander  (1963,  p.  23)  considers 
that  the  overall  agreement  of  the  model  for  the  fundamental 
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mode  Love  and  Rayleigh  waves,  the  higher  Love  modes,  and  the 
second  higher  Rayleigh  mode  is  good,  and  that  the  shape  of 
the  theoretical  dispersion  curve  for  the  first  higher  Rayleigh 
mode  agrees  satisfactorily  with  the  observations.  The  level 
of  the  computed  higher-mode  curves  at  the  longer  periods  could 
be  brought  into  better  agreement  with  the  observations  if 
a  correction  for  earth  curvature  were  made  to  the  model. 

Parts  of  Model  35CM2  do  not  agree  well  with  the  recent 
P-wave  data  of  Prodehl  (1969)  or  with  some  of  the  work  of  the 
U.  S.  Geological  Survey.  These  later  results  indicate  a 
mean  crustal  thickness  of  30  km  in  the  Basin  and  Range  prov¬ 
ince.  Furthermore,  the  results  of  Hill  and  Pakiser  (1966) 
suggest  the  existence  of  an  intermediate  layer  in  the  crust 
of  the  Basin  and  Range  province  with  a  P-wave  velocity  of  6.7 
km/sec.  However,  the  shear  velocity  distribution  of  Model 
35CM2  which  is  well  constrained  by  the  surface  wave  dis¬ 
persion  data  could  be  put  together  with  this  P-velocity  struc¬ 
ture  since  surface  wave  dispersion  is  relatively  insensitive 
to  perturbations  in  P-wave  velocity.  We  defer  further  exam¬ 
ination  of  the  correspondence  between  models  obtained  from 
refraction  studies  and  the  model  obtained  from  the  observa¬ 
tion  of  surface  wave  dispersion  across  the  Basin  and  Range 
province  until  a  later  section  of  this  report. 

Gravity  Data 

The  results  of  gravity  investigations  at  Yucca  Flat  and 
contiguous  areas  are  reported  in  Hazelwood  et  al,  (1963)  and 
in  Healey  '(1968  ) .  These  studies  were  conducted  in  order  to 
determine  ’the  structure  of  the  pre-Cenozoic  rocks  which  lie 
beneath  Yucca  Flat.  Consequently,'  it  is  possible  to  determine 
the  thickness  and  configuration  of  the  Cenozoic  volcanic  rocks 
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and  alluvium  which  constitute  the  valley  fill  in  this  area. 

On  a  regional  basis,  the  Bouguer  gravity  map  of  the  United 
States  (Woollard  and  Joesting,  1964)  and  the  Bouguer  gravity 
map  from  the  Transcontinental  Geophysical  Survey  (U.  S.  Geo¬ 
logical  Survey,  1968)  give  data  for  the  Basin  and  Range. 
Figure  58  shows  the  regional  Bouguer  gravity  in  the  vicinity 

of  NTS,  whereas  Figure  59  is  a  detailed  map  of  the  Yucca  Flats 
area. 

Although  the  rocks  outcropping  at  NTS  have  a  wide  range 
of  lithologies,  they  can  be  classified  into  three  groups  on 
the  basis  of  bulk  densities.  The  first  group  consists  of 
the  Precambrian  and  Paleozoic  rocks  which  have  an  overall  den¬ 
sity  range  from  2.49  gm/cm3  for  quartzite  to  2.85  gm/cm3  for 
dolomite,  with  a  mean  density  of  approximately  2.67  gm/cm3 
(Hazelwood  et  al,  1963).  The  Mesozoic  granitic  rocks  can 
conveniently  be  included  in  this  group  also.  The  second  group 
consists  of  the  Cenozoic  lavas  and  welded  tuffs.  These  rocks 
range  in  density  from  1.8  to  2.8  gm/om3,  with  a  mean  bulk  den- 
sity  of  2.40  gm/cm  (Hazelwood  et  al,  1963).  The  third 
group  consists  of  the  non-welded  and  partially  welded  tuffs 
and  the  alluvium  which  taken  together  have  a  mean  density  of 
1.95  gm/cm  .  We  have  incorporated  the  densities  of  the  in¬ 
dividual  rock  types  as  listed  by  Healey  and  Miller  (Hazelwood 
2t  al,  1963)  in  Table  10. 

Figure  59  shows  in  detail  the  Bouguer  gravity  anomalies 
in  the  Yucca  Flat  area.  The  map  taken  from  Healey  (1963)  in¬ 
corporates  the  data  compiled  by  Healey  and  Miller  (1963)  from 
1400  gravity  stations  with  approximately  one  observation  per 
square  km.  The  Bouguer  gravity  has  a  maximum  value  of  -137 
mgals  at  Mine  Mountain  and  a  minimum  of  -185  mgals  over 
Emigrant  Valley.  There  is  a  regional  gradient  of  -0.6  mgal/km 
northwards  across  Yucca  Flat  (Healey  and  Miller,  1963)  which 
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correlates  with  the  trend  of  the  mean  elevation  above  sea- 
level  across  the  Yucca  Flat  area.  In  general,  the  gravity 
highs  correlate  with  the  outcrops  of  the  pre-Cenozoic  rocks 
which  form  the  mountains  which  surround  Yucca  Flat,  whereas 
the  gravity  lows  correlate  with  the  thick  deposits  of  alluv¬ 
ium  beneath  the  flats. 

A  conspicuous  feature  on  the  Bouguer  gravity  map  is  the 
northward  trending  high  which  occurs  over  the  western  half  of 
Yucca  Flat.  A  major  gravity  low  parallels  this  feature  under 
the  eastern  half  of  the  flat,  a  smaller  low  occurring  along 
the  western  side  of  the  flat.  The  gravity  high  is  due  to  a 
buried  ridge  or  horst  approximately  2.4  km  wide  and  about  18 
km  long  which  is  bounded  by  high-angle  faults  (Healey  and 
Miller,  1963)  in  the  pre-Cenozoic  surface.  Geologic  pro¬ 
files  across  the  ridge  (Barnes  et  alj.  1963)  show  the  pre- 
Cenozoic  erosion  surface  is  0.1a  -  U.24  km  below  the  present 
surface  of  the  flat.  The  gravity  low  to  the  east  of  the  high 
defines  a  trough  approximately  32  km  long  and  4.8  km  wide 
containing  up  to  1.4  km  of  Tertiary  tuff  and  Quaternary  alluv¬ 
ium.  This  value  was  determined  from  the  gravity  data  (Healey 
and  Miller,  1963;  Healey,  1968)  substantiated  by  the  seismic 
work  of  Diment  et  al,  (1959)  and  verified  by  drilling.  Table 
1  of  Healey  (1968)  shows  that  the  estimates  based  on  gravity 
of  the  thickness  of  the  valley  fill  at  the  19  drill  holes 
whose  locations  are  shown  in  Figure  59  are  accurate  to  within 
17  percent  indicating  that  the  density  contrasts  between  the 
various  rocks  in  the  area  are  fairly  well  known.  To  the  east 
of  Yucca  Flat,  the  steep  gradients  away  from  the  trough  to¬ 
wards  the  Painted  Ridge  and  Banded  Mountain  are  indicative 
of  faulting.  Healey  and  Miller  (1963)  postulate  that  two  high- 
angle  faults  may  be  present  along  the  eastern  margin  of  the 
trough.  The  gravity  low  to  the  west  of  the  ridge,  according 
to  Healey  and  Miller  (1963),  represents  a  graben-like  structure 
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let  down  between  the  Eleena  Rang-,  and  the  buried  ridge,  in  which 
up  to  1  km  of  Cenozoic  rocks  ha/e  accumulated.  The  Yucca  Fault 
(Barnes  et  al,  1963)  runs  along  the  eastern  flank  of  the  buried 
ridge  for  much  of  its  length. 


Aeromagnetic  Data 


Figure  60  shows  the  total  magnetic  intensity  in  the  re- 
gion  of  NTS.  The  figure  covers  the  same  area  as  the  Bouguer 
gravity  map  (Figure  58)  and  is  taken  from  the  magnetic  map 
compiled  by  Zietz  and  Kirby  as  part  of  the  Transcontinental 
Geographical  Survey  (U.  S.  Geological  Survey,  1968).  The  area 
of  intense  anomalies  both  positive  and  negative  in  the  north¬ 
west  part  of  NTS  corresponds  to  the  outcrop  of  Tertiary  lavas. 
Some  of  these  lavas  have  sufficient  remanent  intensity  (both 
normal  polarity  and  reversed  polarity)  to  account  for  the 
aeromagnetic  anomalies  (Bath,  1968). 


The  large  scale  aeromagnetic  map  of  the  Yucca  Flat  area 
has  been  compiled  by  the  U.  S.  Geological  Survey  from  data 
collected  in  1960  along  profiles  flown  at  2400  meters  with 
.8  km  spacing  (Hazelwood  et  al,  1963).  The  range  in  values 
of  the  total  magnetic  intensity  over  the  Yucca  Flat  area  is 
750  gammas.  There  are  two  prominent  positive  magnetic  anoma¬ 
lies  which  correlate  with  the  Climax  stock  and  the  Gold  Mea¬ 
dows  intrusive.  (These  anomalies  can  be  seen  on  Figure  60). 
Both  are  granitic  bodies  of  Mesozoic  age.  The  two  positive 
anomalies  are  separated  by  a  negative  anomaly  which  Hazelwood 
et  al,  (1963)  suggest  is  at  least  partially  due  to  remanent 
magnetization  in  the  volcanic  rocks  of  the  Rainer  Mesa. 
Allingham  (written  communication,  1967,  in  Hinrichs,  1968) 

confirms  that  this  negative  anomaly  is  caused  by  a  reversely 
polarized  ash-flow  tuff. 
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Over  Yucca  Flat  itself,  the  magnetic  intensity  is 
relatively  uniform  with  a  gentle  southward  gradient  of  about 
10  gammas/km  (Hazelwood  et  al,  1963).  This  suggests  that  the 
area  beneath  the  flat  is  uncomplicated  by  major  intrusives. 
There  is  no  expression  in  the  aeromagnetic  data 
shown  in  Figure  60,  of  the  north-south  trending  ridge  which 
was  found  from  the  interpretation  of  the  gravity  data  (Healey 
and  Miller,  1963).  However  a  detailed  aeromagnetic  map  of  the 
north-central  part  of  Yucca  Flat  published  by  Hinricks  (Figure 
6,  1968)  shows  a  steep  negative  gradient  along  the  eastern  or 
downthrown  side  of  the  Yucca  Fault. 

Heat  Flow  Data 


We  were  unable  to  find  any  published  reports  of  heat  f]  w 
determinations  with  the  Nevada  Test  Site  in  the  geophysical 
literature.  However,  Blackwell  (1969)  has  prepared  a  map  show¬ 
ing  the  distribution  of  heat  flow  determinations  in  the  western 
United  States  (Figure  48).  The  mean  heat  flow  for  the 
Basin  and  Range  province  from  Blackwell's  data  is  approximately 
2.0  ycal/cm2-sec  .  To  'the  west  of  the  Basin  and  Range 
province  in  the  Sierra  Nevada  heat  flow  is  low  (0.9  ycal/cm  sec 
Roy  and  Blackwell,  1966)  the  transition  taking  place  across  a 

narrow  zone  50  km  east  of  the  structural  and  physiographic 
boundary.  To  the  east  of  the  Basin  and  Range,  in  the  Colorado 

Plateaus  heat  flow  values  are  normal. 


The  heat  flow  in  the  western  United  States  correlates  very 

well  with  the  P  velocity  and  the  Bouguer  gravity  measurements 

n 

as  shown  in  Figure  49. 
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Seismicity 


We  searched  the  Preliminary  Determination  of  Epicenter 
cards  published  by  the  U.  S.  Coast  and  Geodetic  Survey  for  the 

period  January  1960  through  January  1968  for  events  occurring 
within  the  34  -  38°N,  114  .  119°„  with  focal  depths  less  ^ 

35  km.  We  found  448  events,  but  only  five  small  (mb<4.5) 
events  actually  occurred  within  NTS  over  this  period  that 
could  not  be  attributed  to  a  man-made  event.  However,  Figure 
61  shows  that  there  were  several  events  to  the  northeast  and 
northwest  of  the  test  site.  In  addition,  movement  along  the 
Yucca  Fault  has  been  observed  after  nuclear  tests  (Williams 
and  Hoover,  1962;  Dickey,  1968  )  as  well  as  extensive  fractur¬ 
ing  along  preferred  directions  parallel  to  joint  sets  in  the 
bedrock  unrelated  to  the  Yucca  Fault  (Barosh,  1968).  Further¬ 
more,  following  a  recent  large  explosion  at  the  Nevada  Test 
Site,  Benham,  significant  microearthquake  activity  was  ob¬ 
served,  (Engdahl  et  al,  1969),  (Healey  and  Hamilton,  1969) 
in  the  vicinity  of  the  shot  point.  Aki  et  al,  (1969)  from 
an  analysis  of  Love  waves  generated  by  Benham  suggest  that 
an  earthquake  was  triggered  about  0.5  seconds  after  the  ex¬ 
plosion.  Figure  62  shows  the  aftershock  pattern  produced  by 
the  Benham  event.  This  pattern  clearly  defines  the  prevailing 
north-south  structural  trend. in  the  NTS  area. 

A  Geophysical  Model  for  the  rrust  Beneath  Yucca  F1af 
Test  Site 


Nuclear  explosions  at  the  Nevada  Test  Site  have  been 
detonated  at  several  locations  within  the  test  area  and  in 
various  geologic  media.  Since  the  rocks  exposed  at  NTS  have 
a  wide  range  of  physical  properties  we  will  attempt  to  con¬ 
struct  a  plausible  model  for  the  crust  beneath  Yucca  Flat  and 
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because  the  shallow  structure  beneath  the  flat  is  itself  com¬ 
plicated,  the  surface  layers  of  the  model  must  necessarily  be 
somewhat  idealized.  However,  because  of  the  uncertainties  as 
to  the  gross  structure  of  the  deeper  part  of  the  crust,  differ¬ 
ences  in  the  near  surface  layers  may  be  considered  as  second 
order  effects  for  our  purposes  even  though  they  may  have  a 
profound  influence  on  calculations  of  yields  from  explosions. 

The  near  surface  structure  beneath  Yucca  Flat  is  fairly 
well  known,  from  well-logs,  seismic,  and  gravity  data.  There¬ 
fore  we  have  chosen  to  base  the  upper  layers  of  our  model  on 
the  structure  and  lithology  logged  by  Williams  and  Cole  (Williams 
et  al,  1963  )  at  drill  hole  UE4a,  the  location  of  wh1*  ch  is 
shown  in  Figure  59.  The  strata  encountered  in  drill  hole  UE4a 
consisted  of  approximately  0.17  km  of  alluvium  overlying  0.82 
km  of  tuff  which  rests  on  the  Paleozoic  rocks.  Therefore,  in 
constructing  model  YF111  we  have  adopted  0.2  km  for  the  thick¬ 
ness  of  layer  1  and  0.8  km  for  the  thickness  of  layer  two. 

Compressional  wave  velocities  in  alluvium  at  NTS  range 
from  a  low  of  0.65  km/sec  at  the  surface  to  1.65  km/sec  for 
water  saturated  alluvium  at  a  depth  of  0.2  km  (Hazelwood  et 
al»  1963).  Therefore  for  our  model  we  have  assumed  a  com¬ 
pressional  wave  velocity  of  1.3  km/sec  for  layer  one.  For 
tuff,  Hazlewood  et  aly  (1963)  found  velocities  ranging  from 
2.3  km/sec  to  4.0  km/sec  depending  upon  lithology.  We  have 
adopted  3.0  km/sec  as  the  velocity  in  layer  two  of  our  model 
YFUl.  Sonic  velocities  in  the  Paleozoic  limestones  and  dolo¬ 
mites  which  underlie  the  tuff  in  drill  hole  UE4a  are  typically 
3,9  to  6.5  km/ sec  (Hazlewood  et  al , .  1963),  Data  from  a 
traverse  across  the  buried  Paleozoic  ridge  which  was  down  the 
center  of  Yucca  Flat  (Hazlewood  et'  al,  1963)  suggest  that  the 
range  of  the  seismic  velocities  may  be  somewhat  lower  than  this. 


The  vertical  thickness  of  the  Paleozoic  section  beneath  Yucca 
Flat  is  not  known. 

The  bulk  densities  of  the  various  rock  types  which  occur 
at  the  Nevada  Test  Site  are  shown  in  Table  10.  The  mean  den¬ 
sity  of  the  alluvium  at  Yucca  Flat  is  1.94  gm/cm3  (Healey  and 
Miller,  1963).  Densities  given  by  Healey  and  Miller  (1963) 
for  the  various  tuff  units  range  from  1.94  gm/cm3  in  the  non- 
welded  tuffs  to  2.50  gm/cm3  fcr  the  welded  tuffs.  We  have 
adopted  a  density  of  2.2  gm/cm3  for  layer  two  of  our  model 
YF111.  Densities  in  the  Pre-cambrian  and  Paleozoic  rocks  range 
from  2.49  for  quartzite  to  2.85  for  dolomite  and  average  2.67 
gm/cm  (Healey  and  Miller,  1963). 

Although  the  structure  of  the  surficial  layers  at  NTS 
is  well  documented,  the  structure  of  the  deeper  parts  of  the 
crust  is  only  poorly  known.  We  have  been  unable  to  find  any 
reports  of  a  deep  seismic  refraction  study  through  the  Yucca 
Flat  area.  However,  four  of  the  seismic  refraction  profiles 
shown  in  Figure  53  have  used  locations  within  NTS  as  shot 
points.  Hence,  we  can  estimate  the  deep  crustal  structure  be¬ 
neath  NTS  by  extrapolation  if  we  assume  no  rapid  variations 
in  crustal  structure  within  a  radius  of  approximately  50  km 
from  the  mean  shot  point  (37°  10'  N,  116°  05'  W).  Extrapo¬ 
lating  the  depth  to  the  Moho  beneath  NTS  from  the  data  obtained 
by  Hill  and  Pakiser  (1966),  gives  a  total  crustal  thickness 
of  30  km.  Diment  et  al,  (1961)  also  retort  that  the  crust 
is  about  30  km  thick  at  NTS.  Data  from  the  unreversed  pro¬ 
file  east  from  NTS  give  a  crustal  thickness  of  about  25  km 
at  NTS  (Ryall  and  Stuart,  1963).  The  thickness  of  the 
crust  beneath  NTS  is  34  km  calculated  from  Model  II  of  Gibbs 
and  Roller  (1966).  The  range  of  their  determinations,  25  - 
34  km  suggests  that  the  total  crustal  thickness  may  change 
within  50  km  of  the  mean  sho.  point.  Therefore,  we  have 
adopted  30  km  as  the  depth  of  the  Moho  in  model  YF111. 
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The  attitude  and  thickness  of  the  intermediate  layer  be¬ 
neath  NTS  is  even  more  uncertain.  Based  on  the  data  of  Hill 
and  Pakiser  (1966)  an  intermediate  crustal  layer  with  a 
velocity  of  6.7  km/sec  occurs  at  a  depth  of  20  km  beneath 
NTS  which  is  consistent  with  the  conclusions  of  Prodehl 
(1969).  Clearly,  the  6.7  km/sec  layer  in  the  profile  east 
of  NTS  is  too  thin  (3.8  km  thickness),  to  appear  as  a  first 
arrival,  and  Ryall  and  Stuart  (1963)  only  identify  it  from 
questionable  secondary  arrivals.  For  the  NTS-Ludlow  reversed 
profile,  Gibbs  and  Roller  (1966)  found  first  arrivals  with 
apparent  velocities  of  6.8  km/sec  from  shots  fired  at  Ludlow 
but  not  from  events  at  NTS.  Nevertheless,  Model  II  of  Gibbs 
and  Roller  (1966)  has  a  20  km  thick  intermediate  with  a  veloc¬ 
ity  6.8  km/sec  at  the  NTS  end  of  the  profile.  However,  some 
sort  of  intermediate  layer  probably  exists  to  the  south  of 
NTS  because  Johnson  (1965)  found  a  8  km  thick  layer  with  a 
velocity  of  7 . 1  km/sec  at  25  km  depth  at  the  point  where  his 
profile  intersects  that  of  Gibbs  and  Roller  (1966).  In ^ con¬ 
structing  model  YF111  we  have  adopted  one  km  for  the  thickness 
of  layer  three,  the  near  surface  (M.5  km/ sec)  layer.  We 
have  divided  the  remainder  of  the  crust  into  a  18  km  thick 
"granitic”  layer  with  a  P-wave  velocity  of  6.7  km/sec.  We  have 
assigned  a  PR  velocity  of  7.8  km/sec,  in  agreement  with  the 
results  of  Healey  and  Ploger  (1965). 

The  shear  wave  velocity  distribution  in  the  Basin  and 
Range  Province  is  also  poorly  known.  Model  35CM2  (Figure  57) 
gives  the  distribution  of  shear-wave  velocity  with  depth  which 
best  fits  the  surface  wave  dispel sion  data  from  the  Basin  and 
Range  Province  (Alexander,  1963).  The  shear-wave  velocity  in 
the  near  surface  layer  is  2.2  km/sec,  which  gives  a  Poisson’s 
ratio  of  0.22  which  we  have  used  to  calculate  the  shear-wave 
velocity  in  layers  one,  two  and  three  of  our  model  YF111.  For 
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6.1  km/sec  or  "granitic"  layer,  Alexander  (1963)  obtained  a 
shear  wave  velocity  of  3. 6  km/ sec  which  corresponds  to  a 
Poisson's  ratio  of  0.23.  We  have  adopted  3.6  km/sec  as  the  S- 
wave  velocity  for  layer  three  in  our  model. 

Model  35MC2  has  an  intermediate  section  25  km  thick  which 
occurs  at  a  depth  of  25  km.  In  this  section  the  shear  velocity 
is  4.1  km/sec  whereas  the  P-wave  velocity  is  7.7  km/sec  which 
gives  a  Poisson's  ratio  of  C.3.  To  make  the  S-wave  portion 
of  model  YF111  compatible  with  the  P-wave  portion  and  the 
seismic  refraction  data,  we  have  assigned  a  shear  velocity  of 
4.0  km/sec  to  layer  5.  This  corresponds  to  a  Poisson’s 
ratio  of  0.23  which  is  equivalent  to  Poisson's  ratio  in  the 
"granitic"  layer.  We  have  adopted  4.5  km/sec  for  Sn  which  is 
in  agreement  with  the  model  STAN- 3  (Kovach  and  Robinson,  1969) 
and  close  to  the  Sn  velocity  given  by  Anderson  and  Julian  (1969). 

To  complete  our  model  we  need  to  know  the  density  for 
layers  four  and  five  of  our  model  and  in  the  upper  mantle. 

For  the  "granitic"  layer,  Alexander  (1963)  found  that  a 
density  of  2.78  gm/cm  was  required  in  this  layer  to  fit  his 
dispersion  data.  This  is  close  to  the  value  of  2.80  gm/cm3 
for  the  density  of  crust  below  5  km  depth  adopted  by  Thompson 
and  Taiwan!  (1964).  Archambeau  and  Davis  (1969)  give  the 
density  of  the  6.7  km/sec  layer  as  2.95  gm/cm  .  We  have 
adopted  this  value  for  our  layer  five.  The  density  of  the 
upper  mantle  given  by  several  authors  is  very  close  to  3.30 

3 

gm/cm  and  so  we  have  incorporated  this  value  in  model  YF111. 
Discussion 


Geological  and  geophysical  investigations  at  the  Nevada 
Test  Sit  ?  over  the  past  decade  show  that  Yucca  Flat  is  an 
intermontane  basin  typical  of  the  Basin  and  Range  physiographic 
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province.  The  flat  itself  consists  of  an  alluvium  filled 
trough  formed  by  a  block  of  Paleozoic  rocks  which  have  been 
let  down  by  north  to  northeast  trending  faults  between  the 
eastwards  tilted  blocks  of  Paleozoic  strata  which  form  the 
mountains  to  the  east  and  west.  Gravity,  seismic  refraction, 
and  bore  hole  surveys  showed  that  the  alluvium  at  Yucca  Flat 
attains  a  maximum  thickness  of  0.61  km.  These  surveys  also 
revealed  the  presence  of  a  north  to  northeast  trending  ridge 
of  Palezoic  rock  in  the  floor  of  the  valley  which  protrudes 
through  the  Cenozoic  volcanics.  This  ridge  together  with  the 
known  fault  (Yucca  Fault)  which  extends  down  the  center  of 
Yucca  llat  indicate  that  the  structure  is  considerably  more 
complex  than  a  simple  graben  and  suggest  that  subsurface 
trough  is  probably  a  vested  graben  or  "graben-in-graben"  struc¬ 
ture  similar  to  that  found  beneath  Dixie  Valley  (Thompson  et  al, 
1967)  . 


In  addition  to  the  seismic  and  gravity  surveys,  aeromag- 
netic  surveys  have  also  been  conducted  at  NTS.  Tne  latter  revealed 
that  the  two  positive  anomalies  to  the  north  of  Yucca  Flat 
correlate  with  two  granitic  intrusive  bodies  of  Mesozoic  age. 

The  negative  anomalies  found  in  the  vicinity  of  Yucca  Flat 
correlate  with  outcrops  of  the  Anozoic  volcanics  and  probably 
result  from  reversely  polarized  flows  within  this  sequence. 

1 

Although  the  near-surface  structure  at  Yucca  Flat  is  1 

known  in  great  detail,  we  have  been  unable  to  find  any  ref¬ 
erences  pertaining  to  the  deeper  structure.  Hence,  we  do  not 


know  the  actual  thickness  of  the  Paleozoic  rocks  beneath 
Yucca  Flat,  for  example.  Therefore,  in  constructing  the  crustal 
models  shown  in  Table  11,  we  have  assigned  them  a  thickness 


of  1  km.  This  choice  is  not  quite  as  arbitrary  as  it  may 
seem,  since  Hill  and  Pakiser  (1966)  indicate  1  km  of  low- 
velocity  "sediments”  to  the  north  of  NTS  and  Alexander  (1963) 
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gives  2  -  3  km  as  the  thickness  for  the  "sedimentary"  layer 
in  the  Basin  and  Range.  At  depths  greater  than  3  km  the 
Paleozoic  rocks,  should  they  persist,  are  probably  indisting¬ 
uishable  from  the  "granitic"  layer  in  terms  of  velocity. 

The  structure  of  the  remainder  of  the  crust  at  NTS  itself 
has  not  been  determined  directly.  None  of  the  seismic  refrac¬ 
tion  profiles  discussed  above  actually  traverses  the  test  site 
and  hence  the  thickness  of  the  crust,  the  number  of  layers  and 
their  velocities  have  to  be  inferred  from  structures  determined 
at  points  50  km  away  or  more.  Projecting  the  structures  back 
to  the  mean  shot  point  at  NTS  may  not  be  valid  if  faults  of 
the  type  proposed  by  Hill  and  Pakiser  (1966)  are  present.  How¬ 
ever,  because  of  the  lack  of  data  we  were  forced  to  make  such 
a  projection  to  construct  our  models.  Projecting  the  structures 
back  to  the  mean  shot  point  at  NTS,  we  obtained  estimates  of 
crustal  thickness  of  NTS  ranging  from  25  km  to  34  km.  The 
most  probable  thickness  is  30  km,  the  two  extreme  values  rep¬ 
resenting  thicknesses  estimated  from  an  unreversed  profile 
(Ryall  and  Stuart,  1963)  and  a  model  in  which  a  20  km  thick 
intermediate  layer  was  assumed  (Gibbs  and  Roller,  1966),  re¬ 
spectively  . 

The  existence  of  an  intermediate  crustal  layer  in  the 
Basin  and  Range  has  still  not  been  conclusively  established. 

This  is  because  identification  of  the  intermediate  layer  relies 
for  the  most  part  on  later  arrivals  on  the  seismograms,  either 
refractions  or,  in  many  instances,  reflections.  In  the  case 
of  a  plane  layered  crust,  composed  of  two  layers  with  thick¬ 
nesses  h1  and  h2,  and  velocities  and  V2,  over  a  half-space 
with  velocity  V3,  the  condition  for  a  refraction  from  the  in¬ 
termediate  layer  (V2)  to  appear  as  a  first-arrival  is 

V3J  -  V22  )**  (V2*V,”  (v32  -  vi2)ls 

h2 Tpvprp-  >  hi  -w1(v,,-vp- Tpv3-vp 
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If  we  consider  the  case  where  =  6.1  km/sec,  =  6.7  km/sec, 
and  V  =  7.8  km  sec,  which  are  velocities  commonly  reported 
for  the  "granitic"  layer,  the  intermediate  layer,  and  ,  we 
see  that  for  a  refraction  from  the  intermediate  layer  to 
appear  as  a  first  arrival 

h^  >  0.53  h^ 

Thus,  for  the  profiles  considered  above,  refractions  from 
the  intermediate  layer  theoretically  should  not  appear  as 
first  arrivals,  with  the  exception  of  the  profile  studied  by 
Gibbs  and  Roller  (1964). 

From  the  crustal  model  for  NTS  obtained  by  Gibbs  and 
Roller  (1964),  and  extending  the  formula  gi''en  above  to  the 
three-layer  case,  it  is  clear  that  refractions  from  the  inter¬ 
mediate  layer  should  occur  on  the  seismograms  from  the  NTS  - 
Ludlow  profile.  The  fact  that  Gibbs  and  Roller  (1964)  do  not 
report  first  arriva) s  from  the  intermediate  layer  suggests 
that  the  20  km  thickness  assigned  to  this  layer  is  excessive. 
However,  Hill  and  Pakiser  (1966)  observe  that  the  upper  boundary 
of  the  intermediate  layer  is  much  more  clearly  defined  from 
recordings  of  chemical  explosions  than  nuclear  explosions, 
which  may  explain  why  Gibbs  and  Roller  (1964)  do  not  find 
6.8  km/sec  arrivals  from  events  at  NTS. 

Although  presence  of  the  intermediate  layer  with  the 
velocity  and  thickness  indicated  in  our  model  YF111  cannot 
be  established  from  first  arrival  data,  there  is  an  increasing 
body  of  evidence  from  secondary  arrivals  which  justifies  its 
inclusion  in  this  model.  In  addition,  the  surface-wave  data 
of  Alexander  (1963)  suggest  an  intermediate  layer  is  present 
although  it  is  difficult  to  reconcile  Alexander's  model 
35CM2  with  the  body-wave  models. 


LOCATION  OF  EVENTS  AT  NTS,  TRAVEL-TIMES  AND 
UPPER  MANTLE  STRUCTURE  FOR  THE  BASIN  AND  RANGE  PROVINCE 

Iri  this  section  of  the  report  we  deal  with  the  location  of 
events  at  NTS  using  the  methods  outlined  in  Section  III  of  the 
report.  We  have  extended  the  work  of  Chiburis  (1968)  so  as  to 
incorporate  all  available  LRSM  stations  and  Vela  Seismological 
' )t ^er vatories  to  facilitate  a  more  direct  comparison  between  the 
network's  performance  locating  events  at  Lake  Superior  and  locat¬ 
ing  events  at  NTS.  We  shall  discuss  the  results  of  previous 
investigations  of  the  travel  times  from  events  at  NTS  and  combine 
these  with  results  from  other  studies  of  the  properties  of  the 
mantle  beneath  the  western  United  States,  in  an  attempt  to  derive 
a  model  for  the  structure  beneath  the  Basin  and  Range  province 
which  is  consistent  with  the  surface  wave  dispersion  data. 


Location  of  NTS  Using  the  LRSM  and  Observatory  Network 


Chiburis  (1968)  has  located  19  events  at  NTS,  with  and 
without  travel- time  anomaly  corrections,  using  the  program 
oHIFi  and  only  LRSM  stations  at  distances  greater  than  17  de¬ 
grees.  However,  the  Herrin-68  travel-times  were  not  available 
at  the  time  this  work  was  completed.  We  have  recomputed  the 
epicenters  of  the  19  events  studied  by  Chiburis,  using  his 
travex-time  data  and  the  Herrin-68  travel-time  tables.  Fig¬ 
ure  64  shows  the  results.  The  travel-time  anomalies  which  we 
used  as  anomaly  corrections  *.re  shown. in  Table  12  .  The  re¬ 
sults  (Table  13  )  show  that  the  mean  location  error  for  all 
19  events  is  21.2  km  without  anomaly  corrections  and  2.0  km 
lor  17  of  the  events  located  using  anomaly  corrections. 


Figure  64  shows  the  computed  locations.  There  is  considerable 
reduction  in  the  area  of  the  95%  confidence  ellipses  and  the 
standard  deviation  of  the  residuals  after  travel-time  anomaly 
corrections  are  made. 

In  locating  the  shots  at  Lake  Superior,  we  found  that 
using  the  whole  LRSM  and  Observatory  network  of  stations  we 
were  able  to  obtain  satisfactory  results*  Therefore,  in  order 
to  afford  a  more  direct  comparison,  we  relocated  the  same  17 
of  the  19  events  at  NTS  discussed  above  using  the  whole  net¬ 
work.  We  took  the  observed  travel-times  from  the  shoe  reports 
covering  these  events  and  rejected  those  readings  designated 
questionable  by  the  original  analyst.  Figure  65  shows  the 
recomputed  epicenters  with  and  without  anomaly  corrections. 

The  mean  location  error  without  anomaly  corrections  is  10.9  km 
(see  Table  14)  or  approximately  half  of  that  when  only  sta¬ 
tions  for  which  A  >17°  were  used.  The  mean  location  error  with 
anomaly  corrections  (given  in  Table  12) is  3.2  km  or  0.9  km 
greater  than  for  the  more  distant  group  of  stations.  This  in¬ 
crease  probably  results  from  possible  instability  of  the  anoma¬ 
lies  at  the  shorter  distances  due  to  the  influence  of  the  com¬ 
plex  velocity  structure  (Archambeau  et  al,  1969)  of  the  crust 
and  upper  mantle  under  the  western  United  States.  The  mean 
area  of  the  95%  confidence  ellipse  is  significantly  less  than 
that  computed  using  the  more  distant  stations  both  with  and 
without  anomaly  corrections.  However,  both  with  and  without 
anomaly  corrections,  the  standard  deviation  of  the  residuals 
is  larger  for  the  whole  network  than  when  only  stations  where 
A  > 17°  were  used.  This  is  probably  due  to  reading  error; 
Chiburis  read  his  own  data,  whereas  we  used  the  shot  report 
travel-times . 

The  azimuthal  coverage  of  the  LRSM  network  (about  130°)  is 
insufficient  to  allow  us  to  apply  Spence  and  Alexander's  method 
to  the  19  events  at  NTS  discussed  above.  Instead  we  simply 


show  (Table  15  )  Spence  and  Alexander’s  results  (Spence 
and  Alexander,  1969)  for  the  locations  of  three  large 
events  at  NTS  relative  to  DUMONT  which  they  computed  using 
data  from  the  WWSSS  network.  For  comparison,  we  have  located 
the  same  three  events  using  the  method  of  Chiburis  (1968)  with 
the  same  travel-time  data.  The  results  are  interesting  in 
that  although  the  mean  location  errors  are  about  the  same, 
Spence  and  Alexander’s  method  yields  poor  results  in  terms 
of  the  azimuth  for  the  computed  location  of  BOXCAR,  whereas 
Chiburis ’s  method  gives  a  poor  location  for  FAULTLESS.  We 
do  not  know  the  reason  for  this  discrepancy. 

Travel-Times  and  Upper  Mantle  Structure  for  the  Basin  and 
Range  Province 

P-wave,  travel-times,  amplitudes,  phase  velocities,  and 
dT/d A  in  the  western  United  States  have  been  intensively 
studied  by  Niazi  and  Anderson  (1965),  Johnson  (1967),  Julian 
and  Anderson  (1968),  and  Archambeau  et  al,  (1966,1967,1969). 
Velocity  models  of  the  upper  mantle  derived  by  these  authors 
show  rapid  increases  in  velocity  in  the  depth  ranges  375  - 
650  km  and  650  to  700  km.  Anderson  (1967)  and  Fujisawa  (1968) 
have  discussed  some  of  the  implications  of  these  rapid  in¬ 
creases  in  velocity  in  terms  of  composition  of  the  upper 
mantle.  However,  despite  their  similarities  at  depth,  the 
models  show  some  variations,  particularly  in  the  region  of 
the  low  velocity  zone.  In  a  comprehensive  study  of  travel- 
times,  amplitudes,  and  phase  velocities  of  P-waves  from 
the  Bilby  and  Shoal  explosions  and  the  Fallon,  Nevada,  earth¬ 
quake,  Archambeau  et  al,  (1969)  show  that  the  various  physio¬ 
graphic  provinces  correlate  with  differences  in  structure 
within  the  upper  mantle.  Archambeau  et  al,  (1969)  distinguish 
two  basic  types  of  upper  mantle  structures  in  the  tectonic 
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region  of  western  North  America.  These  are  a  Basin  and 
Range  or  depth  zone  mantle  structure  and  a  plateau-mountain 
mantle  structure.  As  more  data  becomes  available, 
further  subdivisions  of  the  latter  structure  may  be  nec- 
ess-ry.  Since  the  Nevada  Test  Site  is  located  in  the  Basin 
am  Range  province,  we  are  concerned  chiefly  with  Archambeau, 
et  a]/s  model  CIT  111P  (Figure  46),  although  we  will  briefly 
discuss  the  other  models.  Figure  45  shows  the  location  of  the 
profiles  which  Archambeau,  et  a]  used  to  study  the  travel- 
times,  amplitudes,  and  phase  velocities  from  these  events. 

The  Basin  and  Range  model  CIT  111P  (Figure  46)  was 
derived  from  data  recorded  along  the  Shoal-rallon  southeast 
profile  (Figure  45);  Shoal  was  a  nuclear  explosion  detonated 
at  the  site  of  the  Fallon,  Nevada,  earthquake  of  20  July  1962, 
approximately  500  km  northwest  of  NTS.  Figure  66  shows  the 
travel- time  data  for  this  profile  together  with  the  theo¬ 
retical  travel- time  curve  computed  from  model  CIT  111P. 

As  Archambeau,  et  alj  (1969)  point  out,  the  fit  of  the  first 
arrival  branches  of  the  travel-time  curve  is  good.  In  addition, 
the  later  arrivals  correlate  with  the  appropriate  extensions 
of  the  first  arrival  branches  of  the  theoretical  curves.  This 
is  important,  particularly  in  the  neighborhood  of  cusps  where 
more  than  one  phase  is  arriving,  because  the  correlation  of 
phases  between  records  was  achieved  independently  of  any 
travel-time  curve.  (For  details  of  the  digital  processes 
used  see  Archambeau  et  al,  1969). 

The  P  velocity  in  model  CIT  111P  (Figure  46)  is  7.7 
n 

km/sec  which  is  consistent  with  the  Pn  velocity  of  7.7  - 
7.8  km/sec  in  the  Basin  and  Range  obtained  from  the  refraction 
surveys  discussed  above.  Furthermore,  the  P-wave  delays 
through  the  upper  mantle  calculated  from  model  CIT  111P  are 
consistent  with  the  observed  P-wave  delays  for  the  Basin  and 


Range  Province  of  0.6  -  0.7  seconds  with  respect  to  the 
Jef freys-Bullen  travel-times  (Figure  67)  found  by  Cleary 
and  Hales  (1966). 

In  addition  to  satisfying  the  travel-time  data,  model 
CIT  111P  also  quantitatively  explains  the  amplitude  observations. 
Figure  68  shows  the  observed  spectral  amplitudes  at  1.5  Hz  ob¬ 
tained  by  Archambcau,  et  al,  (1969)  from  the  Shoal  northeast 
profile  together  with  the  theoretical  amplitudes  calculated 
trom  the  model  CIT  111P.  (The  spectral  amplitudes  from  the 
Fallon  earthquake  were  not  considered  because  of  uncertainty 
as  to  the  nature  of  the  source  spectrum  and  radiation  pattern). 
The  upper  part  of  Figure  68  shows  the  "inferred"  curves  of 
amplitude  versus  distance,  so  called  because  they  represent 
thv?  mean  variation  of  amplitude  with  distance  for  each  branch. 

The  large  deviations  in  the  observed  amplitudes  from  these 
mean  curves  are  attributed  to  differences  in  the  crustal 
transfer  function  at  the  various  recording  stations  by 
Archambeau,  et  al,  (1969),  and  to  -’nterference  effects 
where  two  phases  arrive  nearly  simultaneously.  However,  the 
mean  trend  of  the  inferred  curves  agrees  well  with  the  theoret¬ 
ical  amplitude-distance  curves  computed  from  model  CIV  111P, 
as  shown  in  the  lower  part  of  Figure  68.  Because  of  this  fit, 
Archambeau,  et  al,  (1969)  conclude  that  the  velocity  gradients 
in  model  CIT  111P  are  indicative  of  those  which  exist  in  the 
upper  mantle  and  that  rapid  velocity  variations  near  150,  375 
and  6b0  km  depths  arc  required  to  explain  the  observations. 

Archambeau  et  al,  (1969)  have  also  derived  a  Q-structure 
for  the  Basin  and  Range  province.  They  determined  the  ob¬ 
served  spectral  amplitudes  from  Shoal  end  Bilby  at  three  fre¬ 
quencies,  0.5,  1.0  and  1.5  Hz  and  observed  that  the  high- 
frequency  amplitudes  decreased  most  rapidly  with  increasing 
distance.  To  explain  this  observation,  Archambeau  et  al. 
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(1969)  constructed  the  Q-model  for  a  frequency  of  1.5  Hz 
shown  in  Figure  69,  for  the  lower  frequencies,  they  in¬ 
ferred  a  Q-model  from  the  long-period  surface  wave  data  pub¬ 
lished  by  Anderson  et  al,  (1965).  In  view  of  the  uncer¬ 
tainties  involved  in  determining  Q  either  from  body  wave  or 
surface  wave  data,  we  must  consider  the  two  Q-structuros 
shown  in  Figure  69  as  placing  rough  bounds  on  the  possible 
range  of  variation  of  Q  with  depth  and  only  suggest  that  Q 
at  high  frequencies  is  greater  than  at  low  frequencies. 
Nevertheless,  the  combined  amplitude  data  from  Bilby  and 
Shoal  are  satisfactorily  explained  by  Qcrust  s  1000, 

^mantle  r  •,0°  at  f  *  1<2S  Hz  and  Qcrust  =  1000>  Qmantle  = 

3CU  at  f  =  0.875  Hz. 

[n  addition  to  the  Fallon-Shoal  southeast  profile, 
Arch^mbeau  et  al,  (1969)  also  analyzed  travel-time  and  amplitude 
data  from  P-waves  recorded  alon,*  three  additional  profiles 
(Figure  M5).  From  the  Shoal-Fa ilon  northeast  profile,  they 
derived  a  model  for  the  upper  mantle  beneath  the  Snake  River 
Plains  -  Northern  Rockies  province,  model  CIT  112P;  from  the 
Bilby  southeast  profile,  model  CIT  HOP  for  the  Colorado  Plateau 
Rocky  Mountains;  from  the  Bilby  northeast  profile,  model  CIT 
109P  for  the  eastern  Basin  and  Range  -  Northern  Rockies.  Fig¬ 
ure  70  shows  all  three  of  these  models,  together  with  CIT  111P. 
Below  200  km  depth  all  four  models  are  quite  similar,  the  major 
difference  being  in  the  upper  200  km  of  the  mantle.  For  the 
three  plateau-mountain  models,  there  is  a  high  velocity  (8.0 
km/sec)  lid  to  the  upper  mantle  30  to  60  km  thick  below  which 
the  P-wave  velocity  decreases  rapidly  reaching  a  minimum  at 
a  depth  of  90  km.  The  velocity  then  increases  with  depth  to 
the  bottom  of  the  low  velocity  zone  at  about  45  km  depth 
where  there  is  a  region  of  large  velocity  gradients  after  which 
the  velocity  again  increases  with  depth.  The  Basin  and  Range 
model  (CIT  111P)  on  the  other  hand  has  no  high  velocity  lid, 


but  rather  the  velocity  remains  nearly  constant  at  7,7  km/ 
sec  from  the  base  of  the  crust  to  a  depth  of  135  km  where 
again  there  is  a  region  of  large  velocity  gradients. 

The  results  obtained  by  Archambeau  et  al#  (1969)  are 
in  good  agreement  with  the  earlier  results  of  Johnson  (1967) 
who  studied  dT/dA  from  short-period  P-waves  in  the  distance 
range  3° £  A <  30°  recorded  at  the  extended  TFO  array  in 
Arizona.  The  upper  mantle  model  obtained  by  Johnson  (1967), 
model  CIT  2 0M  in  Figure  71,  is  essentially  an  average  of  the 
models  obtained  by  Archambeau  et  al,  (1969).  This  can  be 
seen  by  considering  Figure  72  which  illustrates  the  excellent 
fit  of  the  theoretical  dT/dA  curves  calculated  from  Archambeau 
et  al’s  two  extreme  models  CIT  111P  and  CIT  109P  with  Johnson's 
(1967)  data.  Model  CIT  109P  is  also  similar  to  the  model 
NTS-NE19  (Figure  71)  obtained  by  Julian  and  Anderson  (1968)  from 
tnvel-times  recorded  along  a  profile  northeast  from  NTS. 

• 

Kovach  and  Robinson  (1969)  studied  dT/dA  measurements 
for  long-period  S-waves  from  earthquakes  in  the  distance 
range  1M°<A<M0;>  recorded  at  the  extended  TFO  array  in  Arizona, 
Since  these  results  were  obtained  in  a  manner  analogous 
to  that  used  by  Johnson  (1967)  to  study  the  P-wave  velocity 
structure  of  the  upper  mantle  under  the  western  United  States, 
Kovach  and  Robinson's  model  should  be  consistent  with  model 
CIT  20m .  Figure  M7  shows  the  preferred  STAN 3  model  obtained 
by  Kovach  and  Robinson  (1969)  together  with  a  simpler  model 
STAN2  which,  although  it  fits  the  observed  data  fairly  well, 
results  in  poorer  agreement  with  the  observed  travel-time  re¬ 
siduals  with  respect  to  the  Jeff reys-Bullen  travel-time 
tables  reported  by  Doyle  and  Hales  >1967)  for  the  distance 
range  28°<A<82°.  The  main  features  of  model  STAN3  are  a  thin 
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high  velocity  6.5  km/sec  lid  capping  a  broad  low  velocity 
zone  centered  at  100  km  depth,  and  regions  of  rapid  increases 
in  velocity  at  depths  of  360  km  and  620  km.  An  important 
difference  between  STAN3  and  CIT  100  series  of  models  (Archambeau, 
et  al ,  1969)  is  the  absence  of  the  zone  of  rapid  increase  in 
shear  velocity  at  a  depth  of  165  km  the  depth  at  which  there  is 
an  increase  in  P-wave  velocity  in  the  latter  models  and  the 
model  CIT  204  (Johnson,  1967).  However,  the  velocity  gradient 
in  model  STAN 3  does  change  at  this  depth  (Kovach  and  Robinson, 
1969). 

Julian  and  Anderson  (1968)  also  present  a  model  CIT  19B 
for  both  P-  and  S-wave  velocity  structure  in  the  western 
United  States  based  on  Love  wave  dispersion  data.  The  P-wave 
velocities  were  derived  from  the  S-wave  data  using  the  Poisson's 
ratio  distribution  appropriate  to  Gutenberg's  model.  Figure 
75  shows  model  CIT  19B.  Whereas  the  overall  S-wave  velocity 
distribution  of  model  CIT  19B  is  similar  to  that  of  STAN3 , 
the  velocity  gradients  below  the  low  velocity  none  in  CIT  19B 
are  srall  except  in  the  regions  350  -  450  km  and  650  -  750  km. 
Also,  the  regions  of  high  velocity  gradients  are  higher  for 
the  surface  wave  model,  and  in  the  case  of  linear  region  of 
high  velocity  gradients,  deeper  than  those  for  the  body  wave 

model  STAN 3. 
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COMPARISON  OF  THE  LAKE  SUPERIOR  AND 
NEVADA  TEST  SITE  SOURCE  REGIONS 


In  preceding  sections  of  this  Report  we  have  discussed 
the  diverse  data  pertaining  to  Lake  Superior  and  the  Nevada 
Test  Site  as  seismic  source  regions  separately.  In  this  sec¬ 
tion  of  the  report  we  compare  and  contrast  the  structure  of 
the  crust  and  upper  mantle  at  Lake  Superior  as  it  is  known, 
with  the  structure  at  NTS.  This  section  will  also  serve  as 
an  overall  summary  to  the  preceding  sections. 

At  Lake  Superior  we  are  dealing  with  a  geologically  old 
and  stable  regime,  whereas  at  NTS  we  are  concerned  with  a 
young,  tectonically  active  region.  The  structure  of  the  crust 
at  Lake  Superior  has  not  changed  radically  over  the  ast  .500 
million  years.  On  the  other  hand,  the  crust  at  NTS  is  still 
being  subjected  to  tectonic  deformation.  This  continuing  re¬ 
lease  of  stress  at  NTS  is  manifest  in  the  high  level  of  seis¬ 
mic  activity  in  the  Basin  and  Range  province  and  at  NTS  itself 
(Figures  61  and  62).  The  Lake  Superior  region  is  essentially 
aseismic. 

The  crustal  structure  at  Lake  Superior  is  complex  despite 
its  antiquity.  All  interpretations  based  on  data  collected 
during  1963  Lake  Superior  Seismic  Experiment  indicate  consid¬ 
erable  variation  (50%  to  100%^  in  crustal  thickness  in  the 
region.  Crustal  thicknesses  exceeding  50  km  were  reported 
indicating  that  Lake  Superior  is  not  a  typical  portion  of  the 
Canadian  Shield.  Unusually  high  velocities,  6.6  km/sec  to  6.9 
km/sec  at  depths  of  only  10  km,  were  found,  further  emphasising 
the  abnormal  nature  of  the  crust  beneath  the  Lake.  These  high 
velocities  are  thought  to  be  due  to  metamorphism  of  the 
Keweenawan  lavas  through  the  amphibolite  facies  in  the  higher 
parts  of  the  crust.  The  cause  of  the  high  velocities  at  depth 


is  not  known,  but  is  probably  related  to  the  emplacement  of 
the  lavas . 

The  shallow  structure  of  the  Lake  Superior  syncline  sup¬ 
ports  this  hypothesis,  at  least  in  part.  In  the  western  half 
of  the  Lake,  the  Keweenawan  lavas  are  upthrust  to  form  a  horst 
which  appears  to  be  genetically  related  to  the  Mid-Continental 
Gravity  High;  Keweenawan  lavas  having  been  reported  from  bore¬ 
holes  along  the  axis  of  this  feature  together  with  velocities 
similar  to  those  found  at  Lake  Superior.  The  shallow  layers  of 
the  crust  at  Lake  Superior  form  a  broad  syncline,  the  axis  of 
which  swings  around  the  Keweenaw  Peninsula.  Whether  or  not  the 
syncline  dies  out  at  the  southeast  shore  of  Lake  Superior  is 
uncertain. 

The  shallow  structure  at  NTS  is  much  more  complex  than  at 
Lake  Superior.  Whereas  the  wavelength  of  the  near  surface  fea¬ 
tures  at  Lake  Superior  are  of  the  order  of  150  km,  the  wavelength 
3t  NTS  is  approximately  20  km.  The  near  surface  layers  at  NTS 
are  repeated  by  thrust  faulting  and  pervaded  by  Cenozoic  and 
Recent  normal  faults.  Tilting  of  the  faulted  blocks  has  resulted 
in  a  series  of  northerly  trending  mountain  ranges  separated  by 
deep  intermontane  valleys  filled  with  sediments.  Movemem 
along  the  faults  is  still  taking  place. 

The  crust  at  NTS  is  thin  compared  to  that  at  Lake  Superior. 
The  crustal  thickness  at  NTS  is  approximately  30  km.  In  con¬ 
trast  to  the  Lake  Superior  region,  the  crustal  thickness  is 
relatively  uniform  throughout  the  Basin  and  Range  province. 

Seismic  refraction  surveys  in  the  Basin  and  Range  province  re¬ 
veal  a  thick  "granitic”  layer  with  a  P-wave  velocity  of  about 

6.1  km/sec.  The  layer  occurs  at  depths  of  2  to  5  km  and  has  a 

mean  thickness  of  approximately  20  km.  The  exact  thickness  of 
this  layer  is  somewhat  uncertain  because  of  the  indications  of 
an  "intermediate”  layer  with  a  velocity  of  6.7  km/ sec  occurring 


at  +he  base  of  this  crustal  layer.  However,  this  intermedi¬ 
ate  layer  is  quite  variable  in  thickness  and  character  over 
the  basin  and  Range  province  and  its  exact  nature  is  some¬ 
what  uncertain.  Nevertheless,  it  is  clear  that  the  mean  crus¬ 
tal  velocity  at  NTS  is  considerably  less  than  at  Lake  Superior. 
Figure  73  illustrates  this  point. 

There  are  considerable  differences  between  the  properties 
of  the  upper  mantle  at  Lake  Superior  and  those  at  NTS.  Bouguer 
anomalies  over  the  Basin  and  Range  average  about  -120  mgal.  At 
Lake  Superior  they  are  sO,  and  values  in  the  range  -40  to  -60 
mgals  are  typical  of  the  Canadian  Shield  as  a  whole.  In  the 
Basin  and  Range  province  heat  flow  is  high  (>2. 0ycal/cm2-sec) 
whereas  at  Lake  Superior  it  is  low  (0.9pcal/cm2  sec).  The  high 
heat  flow  in  the  Basin  and  Range  is  not  due  to  excessive  heat 
generation  in  the  crust.  Rather,  the  work  of  Roy  et  al  (1968) 
has  demonstrated  that  it  is  due  to  anomalously  high  heat  flow  in 
the  upper  mantle.  This  high  mantle  heat  flow  correlates  with 
high  electrical  conductivities  at  shallow  depths  in  the  upper 
mantle.  Together  these  observations  suggest  partial  melting 
in  the  upper  mantle  bent ath  NTS  and  the  Basin  and  Range,  but 

not  at  Lake  Superior.  jhis  is  consistent  with  the  seismic  data 
from  both  areas. 

Regional  differences  in  P~ delays  and  P^-velocities  in  North 
America  have  been  known  for  some  time.  Figure  67  [Figure  3  of 
Archambeau  et  al  (1969)]  shows  the  P-wave  travel- time  anomalies 
with  respect  to  the  Jeff reys-Bullen  tables  (Jeffreys  and  Bullen, 
1940)  at  North  American  stations  obtained  by  Cleary  and  Hales 
(1966).  The  P  travel-time  anomalies  observed  at  NTS  range  from 
+0.2  to  +0.6  sec  indicating  the  arrivals  are  consistently  late. 

At  Lake  Superior,  arrivals  are  early,  the  travel-time  anomalies 
ranging  from  -0.6  sec  to  -1.0  sec.  The  differences  in  P-wave 
delays  are  reflected  in  the  P^  velocities  (Figure  49).  In  the 
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Archambcau  et  a),  (1969)  have  lucidly  demonstrated  the 

relationships  between  P-wave  velocity  structure  in  the  upper 
mantle  and  the  physiographic  provinces  of  the  western  United 
States.  The  travel-time  data  from  Project  Early  Rise,  when 
inverted  to  yield  velocity  structures  for  the  upper  mantle 
suggests  that  similar  relationships  exist  for  the  eastern 
United  States  and  the  Canadian  Shield.  Figure  74  shows  the 
model  CIT  IIIP  for  the  Basin  and  Range  province  (Archambcau 

j 

et  al j  1969)  together  with  the  Canadian  Shield  model  of  Barr 
(1967),  the  model  for  the  mantle  province  west  of  Lake  Superior 
(Lewis  and  Meyer,  1968)  and  the  model  for  the  central  United 
States  (Green  and  Hales,  1968). 

The  most  important  difference  between  the  P-wave  models 
shown  in  Figure  74  is  the  structure  of  the  upper  mantle  between 
50  km  and  150  km.  The  Basin  and  Range  model  (CIT  IIIP) 
has  a  zone  of  low-velocities  in  this  range  with  a  region  of 
rapid  increase  in  velocity  from  145  to  150  km.  None  of  the 
Early  Rise  models  show  this.  Instead  all  three  models  indi¬ 
cate  sharp  but  relatively  small  increases  in  velocity  at 
depths  of  less  than  100  km.  This  Lewis  and  Meyer  model  is  the 
only  model  in  which  the  velocity  decreases  with  depth  in  the 
range  140  to  160  km.  Below  300  km,  the  Lewis  and  Meyer  model 
differs  significantly  from  the  other  models.  This  model  has 
a  thin  region  of  high  velocity  gradients  centered  at  340  km, 
with  a  second  region  with  much  higher  gradients  at  445  km. 

The  Barr  model  has  a  region  of  high  velocity  gradients  start-  j 

ing  at  370  km,  as  do  the  two  models  for  the  Basin  and  Range.  | 

(Note  that  only  the  upper  260  km  of  the  Green  and  Hales  model  I 

is  based  on  Early  Rise  data;  the  deeper  part  of  this  model  was 
based  on  travel-times  from  events  at  NTS). 


Thus  the  main  difference  in  the  P-wave  models  for  the  two 
regions  is  the  apparent  absence  of  a  low  velocity  zone  in  the 
upper  150  km  of  the  mantle  beneath  the  Canadian  Shield.  It 
is  not  clear  at  precisely  what  depth  the  "MOO  km”  discontinuity 
occurs  for  the  Canadian  Shield.  On  the  basis  of  the  Lewis 
and  Heyer  model,  which  is  for  a  path  across  the  shield  and  the 
northern  Rockies,  this  discontinuity  is  50  km  deeper  than  in 
the  Basin  and  Range  mantle  province. 

We  can  also  compare  the  shear-wjve  models  for  the  Canadian 
Shield  and  the  western  United  States  although  the  S-wave  struc¬ 
tures  are  less  well  known.  Figure  75  shows  the  S-wave  structure 
obtained  by  Lewis  and  Heyer  (1968)  from  body-waves  and  the  S- 
wave  portion  of  CANSD  together  with  STAN- 3 (Kovach  and  Robinscn, 
19S9)  and  the  S-wave  portion  of  Alexander’s  model  32  CM2.  The 
most  important  difference  between  the  structures  is  the  pre¬ 
sence  of  the  deep  low  velocity  zone  for  S-waves  which  occurs 
between  125  and  180  km  depth  in  the  Canadian  Shield.  Such  a 
feature  is  not  present  in  the  Kovach  and  Robinson  model,  which 
indicates  that  the  S-wave  velocity  increases  smoothly  through 
the  region.  The  only  other  S-wave  model  with  a  low-velocity 
zone  similar  to  this,  is  that  obtained  by  Ibrahim  and  Nuttli 
(1967)  from  travel-times  to  stations  distributed  over  the 
entire  United  States.  However,  there  is  not  a  low-velocity 
zone  at  this  depth  in  the  model  US-26  derived  by  Anderson  and 
Julian  (1969)  from  a  re-^xamination  of  this  same  data.  Another 
difference  between  the  Lewis  and  Meyer  model  and  STAN- 3  is 
the  velocity  structure  around  600  km.  STAN- 3  which  was  made 
to  be  compatible  with  the  P-wave  model  of  Johnson  (1967), 
has  a  region  of  large  velocity  gradients  centered  at  400  km 
(as  does  US  26),  whereas  the  Lewis  and  Heyer  S-model  has  a 
region  of  large  velocity  gradients  at  450  km,  also  chosen  to  be 
compatible  with  the  corresponding  P-wave  model. 


The  S-models  from  the  Canadian  Shield  (CANSD)  and  the 
Basin  and  Range  <35  CH2)  derived  from  surface-wave  dispension 
studies  offer  only  a  very  limited  basis  for  comparison.  CANSD 
models  suggest  that  the  S-wave  velocity  in  the  mantle  beneath 
the  Canadian  Shield  stays  constant  nearly  to  a  depth  of  135 
km  below  which  there  is  a  broad  low  velocity  zone  extending 

•  to  a  depth  of  315  km.  Model  35  CM2  suggests  that  the  S-wave 
velocity  decreases  beneath  the  Basin  and  Range  between  the 

* 

depths  of  70  and  125  km. 

The  differences  in  the  properties  of  the  upper  mantle 
between  Lake  Superior  are  conveniently  summarized  in  Figure 
76.  This  figure,  based  primarly  on  the  CIT  100  series  of 
models  traces  the  low-velocity  zone  from  the  Basin  and  Range, 
where  it  is  thick  and  extends  to  the  base  of  the  crust, 
north-eastwards  across  the  United  States  to  the  Superior 
Highlands,  where  it  is  either  very  thin  or  does  not  exist. 

A  similar  section  to  that  shown  in  Figure  76  was  derived 
independently  by  Green  and  Hales  (1968). 


IMPLICATIONS 


The  seismological  differences,  defined  by  the  study, 
between  the  Lake  Superior  Region  and  the  Nevada  Test  Site, 
underline  the  need  to  calibrate  individually  all  of  the 
earth’s  geological  provinces  where  test  sites  may  be 
located.  Every  major  geophysical  parameter  affecting  seis¬ 
mic  wave  propagation  appears  to  differ  between  the  two 
regions;  these  parameters  include  seismic  velocities, 
attenuation,  and  regional  stress  patterns.  The  differences 
between  the  two  regions  extend  from  the  near  surface  layers 
to  at  least  200  km  depth  in  the  upper  mantle  (and  possibly 
even  deeper)  so  that  seismic  observations  at  all  distances 
can  be  expected  to  differ  substantially  for  identical 
sources.  Consequently  it  is  extremely  unlikely  that  the 
results*  for  NTS  events,  tnat  depend  on  structure  parameters 

can  be  generalized  to  apply  to  similar  sources  occuring 
in  regions  as  dissimilar  as  Lake  Superior. 

These  differences  could  have  been  predicted  on  an  a 
priori  basis  from  the  models  derived  from  surface  geological 
studies  which  place  Lake  Superior  on  an  ancient  shield  and 
NTS  in  a  relatively  modern  mountain  building  region.  However, 
in  view  of  the  work  of  Archambeau  et  al  (1969)  in  the  Western 
United  States,  it  is  important  that  the  ’’fine  structure"  of 
the  upper  mantle  as  well  as  the  crust  be  determined  precisely 
because  of  their  influence  on  such  important  identification 
criteria  as  source  depth  and  location,  Mg  vs.  m^,  and  com¬ 
plexity.  This  detail  is  still  largely  lacking  for  the  shield 
area  containing  Lake  Superior.  We  can  infer  however  that  the 
seemingly  low  intrinsic  absorption  in  the  shield  area  would 
tend  to  sharpen  the  contrast  between  earthquakes  and  explosions 
relative  to  the  NTS  area,  because  spectral  differences  at  the 
higher  frequencies  could  be  observed  at  greater  teleseismic 


distances.  Moreover,  the  signal  levels  at  any  given  distance 
should  be  uniformly  higher  for  shield  events  of  a  given  size. 

In  addition,  the  lack  of  any  strong  regional  stress  in  the 
shield  would  preclude  significant  low  frequency  contributions 
to  the  explosion  signal  due  to  stress  release,  such  as  seem 
to  occur  at  the  Nevada  Test  Site. 

We  believe  that  the  promising  new  computational  techniques 
now  under  envelopment  for  wave  propagation  in  media  with  compli¬ 
cated  lateral  and  vertical  variations  in  physical  properties 
can  be  used  to  predict  qualitatively  the  seismic  signals  for 
earthquakes  and  explosions  in  dissimilar  source  regions,  pro¬ 
vided  significant  detailed  knowledge  of  the  crust  and  upper 
mantle  structure  is  available  .  However,  the  present  study 
demonstrates  how  difficult  it  is  to  obtain  reliable  determina¬ 
tions  of  the  "fine"  structure  in  a  region.  Studies  at  least 
as  comprehensive  as  those  conducted  around  Lake  Superior  and 
NTS  are  required  to  document  a  region  of  interest. 

Despite  thn  gross  differences  in  crustal  and  mantle 
structure  we  found  that  using  the  LRSM  and  observatory  net¬ 
work  the  accuracy  of  absolute  event  locations  in  the  regions 
is  similar,  namely  10  to  20  km.  The  aret  s  of  the  95%  confi¬ 
dence  ellipses  for  both  regions  are  large ,  those  for  NTS 
being  larger  by  a  factor  of  two.  Similarly,  the  accuracy  of 
the  relative  event  locations  are  of  the  same  order  for  both 
regions.  However,  using  the  master  event  technique  the  error 
in  location  is  only  2  to  3  km  and  the  areas  of  the  95%  con¬ 
fidence  ellipses  are  correspondingly  smaller.  Judging  by 
these  results  and  other  experiments  such  as  LONG  SHOT,  these 
location  accuracies  appear  to  be  rather  typical  of  what  to 
expect  for  different  source  regions. 

This  study  clearly  demonstrated  the  necessity  of  a 
"regionalized"  approach  to  classifying  seismic  events.  Each 
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source  region  of  interest  must  be  documented  individually  to 
provide  the  best  chance  of  identifying  a  particular  event  from 
that  region.  There  are  at  least  two  practical  ways  of  doing 
this:  (1/  conduct  detailed  and  comprehensive  geophysical  studies 
of  the  region  similar  to  those  for  Lake  Superior  and  NTS 
discussed  in  this  report,  and  (2)  use  an  empirical  reference 
event  approach  whereby  all  events  from  the  region  are  referred 
to  a  "standard"  event  for  that  region  which  can  be  located  and 
identified  with  certainty.  The  latter  is  considerably  less 
costly,  hut  depends  on  having  enough  large  events  to  cover  the 
changes  in  the  structure  throughout  the  region  of  interest. 

We  believe  that  present  technology  and  computation  capability 
make  it  feasible  to  document  each  source  region  of  interest 
individually  along  these  lines  such  that  rapid  and  complete 
evaluation  of  any  event  can  be  routinely  accomplished.  The 
data  brought  together  in  this  report  can  be  used  as  the  point 
of  departure  to  develop  such  regionalized  event  evaluation 
procedures. 
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TABLE  I 


Members  of  the  North  American  Seismic  Group  participating  in  th 3 
1963  Lake  Superior  Seismic  Experiment. 


INSTITUTION 


PRINCIPAL  CONTRIBUTORS 


Carnegie  Institution  of  Washington 

Dominion  Observatory,  Canada 

Graduate  Research  Center  of  the 
Southwest 

Pennsylvania  State  University 
Princeton  University 
University  of  Alberta 
University  of  Manitobe 
University  of  Michigan 

University  of  Minnesota 
Ihiversity  of  Western  Ontario 
University  of  Wisconsin 
U.S.  Geological  Survey 

Vela  Seismological  Center, 

AFTAC 


L. T.  Aldrich,  I.S.  Sacks,  T.J. 
Smith,  J.S.  Steinhart,  Z. Suzuki 

A. M.  Bancroft,  G.W.  Sander 
K.  Graham,  A.L.  Hales 

B. F.  Howell,  Jr. 

M.  Dunne 

G.  Cummings 
D.  Hall 

D.E.  Wil.lis,  C.G.  Bufe,  G.E. 
Frantti 

P.  Farnham,  T.  Fitch,  H.  Mooney 
R.F.  Mereu 

R.P.  Meyer,  T.J.  Cohen 

J.M.  Healy,  W.H.  Jackson, 

L.C.  Pakiser,  J.C.  Roller, 

B.L.  Tibbetts 

R.H.  Mansfield 


TABLE  2 


ORGANIZATION 


PROriLES  OCCUPILD 


/FTAC  East  Liverpool,  Ohio,  to  South 

Mills,  North  Carolina. 

Chapleau,  Ontario,  to  Glace  Bay, 
_ Nova  Scotia. _ 

Southwest  Center  for  Advanced  Lake  Superior,  Michigan,  to  E] 
Studies  Dorado,  Arkansas. 

Lake  Superior,  Michigan,  to 
Charleston,  West  Virginia. 

Colorado  Springs,  Colorado,  to 

_  Charleston,  West  Virginia. _ 

Stanford  Research  Institute  Chapleau,  Ontario,  to  S 2heff erville 

_  Quebec , _ 

University  of  Alberta  Yorkton,  Saskatchewan  to 

Whitecourt,  Alberta. 

Bov/  Island,  Alberta,  to 

Cold  Lake,  Alberta. _ 

University  of  Manitoba  Lake  of  the  Woods,  Ontario,  to 

_  Lake  Manitoba,  Manitoba. _ 

University  of  Michigan  Bear  Lake,  Michigan,  to 

_ Potters  Hill,  North  Carolina. 

University  of  Toronto _ Southeastern  Ontario. _ 

University  of  Western  Toronto  Lake  Nipigon,  to 

Smooth  Rock  Falls,  Ontario. 

Port  Arthur,  Ontario,  to 

_ Barlzon  Lake,  Manitoba. _ 

University  of  Wisconsin  Bemidji,  Minnesota,  to 

_ _ . _ Vancouver,  Washington 

U.  S.  Geological  Survey  Saint  Cloud,  Minnesota,  to 

Cedar  City,  Utah. 

Bemidji,  Minnesota,  to 
Miles  City,  Montana 

Ledve,  Alberta,  to  Mile  Post  750,  B 

Flin  Tlon,  Manitoba,  to 

_  Yellowknife,  Northwest  Territories . 
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TABLE  4 

Plausible  crustal  models  for  the  EARLY  RISE  shot  point 


LAYER 

THICKNESS 

a 

B 

| 

- 

No. 

km 

km/sec 

km/sec 

gm/cm^ 

t 

TABLE  5 


Travel-time  anomalies  for  the  October  1964  series  of  shots 
in  western  Lake  Superior  calculated  from  shot  49A. 


Station 

Anomaly 

Station 

Anomaly 

GP-MN 

+  3.02 

NL-AZ 

+  5.50 

WF-MN 

+  4.02 

WO-AZ 

+  6.07 

RK-ON 

+  3.42 

HR-AZ 

+  5.09 

VO-IO 

+  4.09 

GE-AZ 

+  5.00 

RY-ND 

+  3.46 

JR-AZ 

+  6.62 

BR-PA 

+  2.63 

LG-AZ 

+  6.13 

BL-WV 

+  1.04 

SG-AZ 

+  5.91 

DH-NY 

+  1.81 

SN-AZ 

+  5.75 

CP008 

+  2.10 

RT-NM 

+  3.15 

LS-NY 

+  2.21 

UB010 

+  3.15 

NP-NT 

+  4.66 

HL2ID 

+  4.21 

EU-AL 

+  2.02 

BMQ03 

+  3.54 

HN-ME 

+  0.05 

LC-NM 

+  5.38 

JE-LA 

+  4.41 

KN-UT 

+  5.50 

SV2QB 

+  0.47 

EK-NV 

+  6.09 

DR-CO 

+  2.60 

TF06  0, 1 

+  5.54 

WM006  Reference  Station  MN-NV 

(Residual  -3.62  seconds) 

+  5.58 
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TABLE  10. 

Description  of  rocks  outcropping  in  vicinity  of  Yucca  Flat 

Nevada  Test  Site  ' 


TABLE  11 


A  plausible  Crustal  Model  for  Yucca  Flat  Area,  NTS 


Layer 

Thickness 

.  ot 

0 

1 

0.2 

1.3 

0.78 

1.96 

•  2 

0.8 

3,0 

1.80 

2.20 

3 

1.0 

4.5 

2.68 

2.67 

4 

18.0 

6.1 

3.60 

2.80 

5 

10.0 

6.7 

4.00 

2.95 

6 

7.8 

4.50 

3.30 

TABLE  12 

Travel-Time  Anomalies  (in  seconds)  for  NTS  Computed  Using 
Herrin  1968  Travel-Time  Tables 


Station 

Anomaly 

Source*  Station 

Anomaly 

Source* 

IV I  -  N  V 

-0.12 

A 

AD  -  IS 

+  0.10 

B 

MV-CL 

-0.41 

A 

AX2AL 

+  1.68 

A 

HL-ID 

+  0.42 

A 

BE-FL 

+  1.16 

A 

HI  2  ID 

+  0.42 

E 

BL-WV 

=  0.56 

A 

AZ-TX 

+  4.80 

A 

BR-PA 

+  0.71 

A 

TK-WA 

+  3.03 

A 

CP008 

+  9.82 

A 

SK-TX 

+  4.62 

A 

DH-NY 

-9.02 

A 

GI-MA 

-0.01 

A 

EB-MT 

+  0.09 

A 

TD-NM 

+  0.32 

A 

EN-MO 

+  0.60 

B 

FR-MA 

+  0.36 

A 

EU-AL 

+  1.98 

A  • 

BX-UT 

-9.67 

A 

GG-GR 

+  0.96 

A 

CP-CL 

-0.79 

A 

HN-ME 

+  0.56 

A 

LAO  10 

+  0.05 

A 

KC-MO 

+  1.70 

D 

WN-SD 

+  1.46 

A 

LZ-BV 

+  0.33 

A 

SE-MA 

+  3.95 

B 

NP-NT 

+  1.00 

A 

RG-SD 

+  0.69 

C 

00 -NW 

-0.04 

A 

IJB010 

+  3.12 

A 

PG-BC 

+  2.69 

F 

DR-CO 

-0.42 

A 

SI-BC 

+  2.78 

A 

* 

BMO03 

+  2.63 

A 

WH2YK 

+  0.94 

F 

LC-NM 

+  2.56 

A 

SV2QB 

-0.41 

A 

RT-NM 

+  0.41 

A 

SV3QB 

-0.41 

E 

KN-UT 

-0.98 

A 

PZ-PR 

+  1.22 

A 

WM006 

+  2.97 

A 

LV-LA 

+  1.95 

A 

TF006 

+  0.55 

A 

RK-ON 

0.0 

Reference  Sta. 

MN-NV 

-0.95 

A 

(Residual  =  +0.70  Seconds) 

*Anomalies  were 

determined 

from  events  as 

indicated  below 

A  -  Tan, 

Bilby 

E  - 

Same  as  previous  station 

B  -  Corduroy,  Buff,  Bronze 

F  - 

Commodore,  Agile,  Nash 

C  -  Piledriver, 

Chartreuse, 

Dumont 

D  -  Greeley 
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+  Confidence  ellipse  does  not  contain  shot  point. 

°  Excluding  BOURBON,  COMMODORE,  and  SCOTCH. 
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TABLE  15 


Location  of  3  Large  Events  at  NTS  Relative 
to  DUMONT  (After  Spence  and  Alexander,  1969) 


EVENT 

Actual 

Location 

Computed 

Location 

Location  Error 

Distance 

Km 

Azimuth 

Degree 

Distance 
.  Km 

Azimuth 

Degree 

Distance 

Km 

Boxcar 

40.8 

300.2 

40.7 

291.5 

6.0 

Halfbeax 

31.2 

316.8 

29.4 

317.5 

1.8 

Faultless 

169.6 

355.4 

170.4 

355.3 

0.7 

Mean 
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bPAP  PROVINCE 
!  B  r  k  p  1  a  i  r. 

II  Coppermine  beltei  plain 

SLAVE  PROVINCE 
III  Bathurat  fold  belt 

CHURCHILL  PROVINCE 
IV  Eaat  am  fold  belt 
V  The Ion  plain 
VI  Athabaaca  plain 
VII  Belcher  fold  belt 
VIII  Car*  Smith  fold  belt 
IX  Labrador  trouqh 
X  Hiataaami  belted  plain 


EIPEPKP  PROVINCE 
XI  Croaa  Lake  telt 
XII  Fr.qliah  B.ver  belt 
X  II I  Qiirrt  ;co  belt 
XIV  Sip;  4or  pi  a  ; 

XV  Cobalt  pla.r, 

SOUTHERN  PROVINCE 
XVI  perckean  fold  belt 
XVII  Port  Arthur  belted  plain 
XVIII  Lake  Superior  baa  n 

3PEWILLE  PROVINCE 
XIX  Lake  Melville  plain 
XA  S^mar.ville  IMd  belt 
XX!  Naakaupt  fold  belt 
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Figure  3.  Structural  map  of  the  Canadian  Shield.  (Adopted 
from  Stockwel1  1962;  1964.) 


4’  1 


VELOCITY  (km/sec.) 


Figure  4.  The  model  CANSD  for  the  crust  and  upper  mantle 

beneath  the  Canadian  shield  derived  from  surface 
wave  dispersion  data.  (Brune  and  Dorman,  1963.) 
Dashed  line  indicates  the  P-wave  model  obtained 
by  Barr  (1967)  from  body  wave  data. 
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Figure  5.  Locations  of  shot  points  and  recording 
stations  for  the  1965  Hudson  Bay  crus¬ 
tal  experiment.  (After  Hobson,  1967.) 
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Figure  11.  t  -x  plot  for  data  recorded  at  Otter 
Cove,  Woodsman  and  Tom.  Theoretical 
travel- times  are  for  reflections  from 
interfaces  showing  in  model.  (After 
O'Brien,  1968.) 
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Figure  12.  Locations  of  stations  and  shot  points  for  the 

University  of  Wisconsin  shallow  seismic  profiling 
study.  (Courtesy  of  R.P.  Meyer,  University  of 
Wisconsin  Geophysical  and  Polar  Research  Center.) 
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Figure  13 


Figure  15.  Contour  map  showing  probable  depth  to  the  Moho  in 
the  Lake  Superior  region.  Data  for  crust  beneath 
Lake  Superior  are  taken  from  O'Brien  (1968). 
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Residual  total  magnetic  intensity  contour 
map  of  eastern  Lake  Superior  showing  la¬ 
belled  anomalies  referred  to  in  the  text. 
(After  Hinze  et  al . ,  1966.) 
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Figure  23.  Earthquakes  of  the  United  States  through  1963 
exclusive  of  California  and  western  Nevada. 
(After  Eppley,  1965.) 


Locations  of  the  38  shots  detonated  during  Project 
Early  Rise.  Arrows  indicate  directions  to  the  shot 
points  beyon  1  kin;  values  in  parentheses  indicate 
the  distance. 
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Figure  32.  Assumed  geometry  at  source  in  Spence  and  Alexander's 
method  for  relative  event  location. 
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REDUCED  TIME,(T 


Reduced  travel- times  obtained  during  Project  Early 
Rise.  Solid  curve  represents  theoretical  travel- 
times  calculated  from  model  NCER  1.  (After  Iyer 
et  al. ,  1969.) 


Figure  35 


Figure  39.  Model  for  the  crust  and  upper  mantle  to  the  west 

SuPerior  obtained  by  Lewis  and  Meyer 
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Figure  41, 


Reduced  travel- times  from  Early  Rise  profile  5 
showing  four  branches  of  travel-time  curve  iden¬ 
tified  by  Willis  (1968).  Figure  courtesy  of  D. 
E.  Willis,  University  of  Michigan. 
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Figure  46.  Model  CIT  1J.1-P  for  P-wave 
structure  in  the  Basin  and 


Range  mantle  province. 
(After  Archambeau  et  al., 
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Tertiary  volcanic  racks 
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Generalized  geologic  mao  of 
the  Yucca  Flat  area,  Nevada 
Test  Site.  (After  Hinrichs, 
1968.) 
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DISTRIBUTION 


MALFPINT  PLATE 

PRESENT  OUTCNOP  AND  INTERNED  SUBSURFACE 
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Structural  map  of  the  Yucca  Flat  area,  Nevada  Test 
Site  showing  major  thrust  plates.  (Adopted  from 
Barnes  et  al.,  1963.) 


Compilation  of  results  from  seismic  refraction 
surveys  in  the  western  United  States  plotted  as 
a  fence  diagram. 
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Figure  54. 


Fundamental  and  higher  mode  Love  wave 
group  velocity  dispersion  data  across 
the  Basin  and  Range  province  obtained 
by  Alexander ,  1963.  Dotted  lines  rep¬ 
resent  theoretical ,  group  velocity 
dispersion  curves  computed  from  model 
35  CM  2.  (After  Alexander,  1963.) 
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Figure  55. 


Fundamental  and  higher  mode  Rayleigh 
wave  group  velocity  dispersion  acoss 
the  Basin  and  Range  province  from  two 
events  in  Utah.  Solid  line  represents 
theoretical  group  velocity  curves  com 
puted  from  model  35  CM2.  (After 
Alexander ,  1963.) 
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Figure  56. 


Fundamental  and  higher  mode  Love  and 
Rayleigh  wave  group  velocity  disper¬ 
sion  across  the  Basin  and  Range  pro¬ 
vince  from  Montana  earthquake  after¬ 
shock  sequence.  Solid  line  represent 
theoretical  curves  computed  from 
modal  35  CM2. (After  Alexander,  1963.) 
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Figure  58 


Regional  Bouguer  gravity  fold  in 
the  vicinity  of  NTS.  (Adapted 
from  Miscellaneous  Geologic  In¬ 
vestigations  Map  1-532-B,  pub¬ 
lished  by  U.S.  Geological  Sur¬ 
vey,  1968.) 
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Figure  59.  Detailed  Bouguer  gravity  anomaly 
map  of  Yucca  Flat,  Nevada  Test 
Site.  Note  that  all  values  are 
negative.  (After  Healey,  1968.) 
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Total  magnetic  intensity  in  the 
vicinity  of  the  Nevada  Test 
(Adapted  from  Miscellaneous  Geo¬ 
logic  Investigations  Map  1-532-A 
published  by  U.S.  Geological  Sur¬ 
vey,  1968.) 


Epicenters  of  shallow  events  in  the  vicinity  of 
the  Nevada  Test  Site  for  the  o^riod  January  1960 
through  January  1968  as  dete  *  ;r  Vned  by  the  U.  S. 
Coast  and  Geodetic  Survey.  (1  distinction  as  to 
number  of  events  occurring  at  jiven  location.) 


15,000  FT 


Natural  faults  •  •  •  Earthquakes  sine*  explosion 
Fractures  produced  by  explosion 


Figure  62.  Locations  of  aftert,hocks  follow¬ 
ing  the  detonation  of  underground 
nuclear  explosions  Benham  on  Dec. 
19,  1968  at  NTS. 
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Figure  66, 


Observed  travel-times  from  the  Shoal-Fallon  south¬ 
east  profile  and  theoretical  travel-time  curve 
(solid  line)  calculated  from  model  CIT  111P.  (After 
Archambeau  et  al. ,  1969.) 
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Figure  69.  Q  structure  for  CIT  100  series 
of  models.  (After  Archambeau 
et  al.,  1969. ) 
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Figure  70.  The  CIT  100  series  of  models 
showing  range  of  upper  mantle 
structures  in  the  western 
United  States.  (After 
Archambeau  et  al..  1969.) 
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Figure  71.  Upper  mantle  P-wave  velocity 
distributions  for  model  CIT 
204  (Johnson  1967)  and  NTS- 
NE19.  (Julian  and  Anderson, 
1968.)  Figure  from  Julian 
and  Anderson,  1968. 
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Figure  74.  P-wave  models  from  Early  Rise  compared  with  CIT111P 
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APPENDIX  I  -  BASIC  OBSERVATION  DATA 


The  purpose  of  this  appendix  is  to  bring  together  the 
first-arrival  data  from  the  LRSM  stations  and  Vela  Observa¬ 
tories  obtained  during  this  study  it  SDL  by  standard  filter 
analysis  procedures  for  the  October,  1964  series  of  shots 
in  western  Lake  Superior.  The  travel-time  data  contained 
m  these  tables  were  used  in  the  location  calculations  of 
Section  III..  Corresponding  data  for  the  NTS  events  (Section 
IV)  are  published  elsewhere  (Chiburis,  1968;  and  the  SDL 
"Shot  Reports"  covering  these  events.  ) 


Origin  times,  relative  and  absolute  locations  of  the  38 
shots  fired  during  Project  EARLY  RISE  are  also  given. 
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1.33 

160.8 

10 

30:00.0 

14 

47 

33 

12. ON 

88 

56 

2  8 . 8W 

0.40 

339.6 

09 

30:00.0 

15 

47 

32 

46. 8N 

88 

56 

31. 8W 

0.45 

206.2 

10 

30:00.0 

16 

47 

33 

8.4N 

88 

56 

2  5. 8W 

0.27 

343.8 

09 

30: 00.0 

17 

47 

32 

58. 2N 

88 

56 

16. 8W 

0.13 

116.2 

10 

30:00.0 

18 

47 

32 

58. 2N 

88 

56 

9 .  OW 

0.28 

101.4 

10 

00:00.0 

19 

47 

33 

1.2N 

88 

56 

18. 6W 

0.08 

063.8 

11 

00:00.0 

20 

47 

33 

4. 8N 

88 

56 

48.  OW 

0.56 

285.4 

09 

30:00.0 

21 

47 

32 

52. 8N 

88 

56 

47. 4W 

0.57 

247.1 

10 

30:00.0 

22 

47 

32 

5  6. 4N 

88 

56 

OW 

0.48 

103.5 

09 

30: 00.0 

23 

47 

32 

55. 2N 

88 

55 

55. 8W 

0.57 

105.0 

10 

30:00.0 

24 

47 

33 

12. ON 

88 

56 

29. 4W 

0.40 

337.9 

09 

30:00.0 

25 

47 

33 

16. 8N 

88 

56 

19. 8W 

0.52 

005.5 

10 

30: 00.0 

26 

47 

33 

17. 4N 

88 

56 

17. 4W 

0.55 

010.6 

09 

30: 00.0 

27 

47 

33 

17. 4N 

88 

56 

19. 2W 

0.54 

006.7 

10 

30: 00.0 

28 

47 

33 

4. 2N 

88 

56 

12. 6W 

0.24 

057.1 

09 

30: 00.0 

29 

47 

33 

3.  ON 

88 

56 

7.8W 

0.31(5) 

032.9 

10 

30: 00.0 

30 

47 

33 

7. 2N 

88 

56 

40. 2W 

0.44 

300.6 

09 

30: 00.0 

31 

47 

33 

11. 4N 

88 

56 

3  5. 4W 

0.45 

321.9 

10 

30: 00.0 

32 

47 

33 

13. 3N 

88 

56 

36. OW 

0.51 

325.9 

09 

30:00.0 

33 

47 

33 

9.6N  , 

88 

56 

1.  2W 

0.53 

056.0 

10 

30:00.0 

35 

47 

32 

5  3. 4N 

88 

56 

16. 2W 

0.26 

150.0 

10 

30:00.0 

36 

47 

33 

13. 8N 

88 

57 

6 .  OW 

1.01 

295.0 

09 

30:00.0 

37 

47 

33 

9.6N 

88 

56 

4  3. 8W 

0.54 

303.3 

10 

30: 00.0 

38 

47 

32 

54. 6N 

88 

56 

8.4W 

0.33 

120.0 

09 

30:00.0 

39 

47 

32 

46. 2N 

88 

56 

4.8W 

0,56 

139.5 

10 

30:00.0 

*  From  mean  shot  point,  47.55°N,  88.94°W 


TABLE  A- 1 

Origin  Times  and  Locations  of  EARLY  RISE  shots 
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TABLE  A- 3. 

First  arrival  data  for  shot  49A 


STATION 


Gl’-MN 


VO-IO 

RY-ND 

BR-PA 

BL-WV 

DII-NY 


LS-NII 


EU-AL 

HN-ME 

RT-NM 

JE-LA 

UBO-IO 

SV-2-QB 

DR- CO 


NL-AZ 

BMO-03 

LC-NM 

KN-UT 

WO-AZ 

HR-AZ 

EK-NV 

GE-AZ 

TFO-60 

JR-AZ 

LG-AZ 

SG-AZ 

RN-AZ 

t»iN-NV 

NP-NT 


DISTANCE 

(KM) 


273.5 

441.9 
474.3’ 

598.9 
873.8 

1213.5 

1273.5 
1313.1 


CP0-08  1356.2 


1435.4 


WMO-06  1582.2. 


1629.3 

1672.5 

1678.2 

1740.2 

1759.5 
1807.7 

1831.6 


HL-2-ID  1956.3 


2066.8 

2127.3 

2182.9 
2204.0 

2214.3 

2240.3 

2270.2 

2296.9 
2304.6 
2308.0 

2311.4 

2332.2 

2363.9 
2495,0 

3484.5 


AZIMUTH 

DEGHEES 


278.1 
208.0 
325.7 
198.0 

279.6 
129.0 

143.5 

110.2 
163.1 
097.4 

210.5 

173.3 
086.5 

230.5 

186.7 

250.8 
054.2 

239.6 

266.6 
239.  1 

272.6 

226.4 

246.9 

238.7 

238.5 

255.9 

236.4 

238.6 

240.9 
239.3 

244.5 

238.5 

256.9 

346.9 


ARRIVAL 
TIME  (Z) 


08:30:41.0 
08:31:02.6 
08:31:06.0 
08:31: 22.0 
08:31:55.2 
08 :  32 :-35 . 9 
08:32:41.6 
08:32:47.2 
08:32:52.7 
08:33:02.4 
08:33:17.9 
08:33:25.5 
08:33:28.7 
08:33:32.2 
08:33:41.0 
08:33:42.2 
08:33:45.0 
08:33:50.1 
08:34:05.8 
08:34:19.2 
08:34:23.6 
08:34:31.3 
08:34:33.6 
08:34:35.2 
08:34:36.9 
08:34:40.9 
08:34:42.4 
08:34:43. 7 
08:34:45.1 
08:34:45.0 
08:34:46.8 
08:34:49.6 
08:35:01.8 
08:36:22.5 


TRAVEL 
TIME  (SEC) 


40.8 
62.4 

65.8 

82.8 
115.0 

155.7 

161.4 
167.0 

172.5 
182.2 

197.7 

205.3 

208.5 

212.0 

220.8 
222.0 

224.8 

229.9 

266.6 
259.0 

263.4 

271.1 

273.4 
275.0 
276.7 

255.9 

282.2 
283*5  . 

284.9 

239.3 

244.5 

289.4 

256.9 
382.3 


AMPLITUDE  PERIOD 


41.8 


30.4 


13.9 

43.1 


15.5 

11.8 


15.2 


13.7 


1.1(5) 


4.74 


4.84 
5.86 
5.18 
5.69 
5.47 
5 ,  ?Q 

5.15 
3.51 
4.39 
3.  58 
4.29 

4.34 
3.44 
3.65 

3.16 
3.41 
3.55 
3.22 
3.28 

3.34 

3.85 
3.92 

3.67 
3.49 

3.68 

3.72 
3.83 
3.26 

3.73 
3.61 
4.21 


t. 


;  *  *  ■  * 


•  * 


TABLE  A -4. 

First  arrival  data  for  shot  45E 


STATION 

DISTANCE 

(KM) 

AZIMUTH 

DEGREES 

GP-MN 

243.6 

280.4 

Wl-MN 

422.4 

203.3 

RK-ON 

|  463.1 

329.1 

VO- 10 

583.4 

195.0 

RY-ND 

844.0 

280.1 

BR-PA 

1234.0 

127.5 

BL-IW 

1287.2 

141.9 

DH-NY 

1340.5 

109.  2 

CP0-08 

1359.3 

161.5 

LS-NH 

j  1465.8 

096.7 

WMO-06 

1560.6 

209.3 

BU-AL 

j  1626.6 

171.9 

HN-ME 

:  1704.4 

086.1 

RT-NM 

1649.7 

229.7 

JE-LA 

1730.1 

185.5 

UBO-10 

1727.6 

250.3 

SV-2-QB 

1837.2 

054.3 

DR -CO 

1801.1 

239.0 

HL-2-ID 

1924.5  . 

266.4 

NL-AZ 

2036.4 

238.5 

BMO-03 

2096.2 

272.5 

LC-NM 

2155.6 

225.7 

KN-UT 

2172.4 

246.5 

WO-AZ 

HR-AZ 

2210.0 

238.0 

EK-NV 

2238.0 

255.6 

GE-AZ 

2267.0 

235.8 

TFO-60 

2274.2 

238.0 

JR-AZ 

2277.3 

240.3 

LG-AZ 

2280.9 

238.8 

SG-AZ 

2300.9 

244.1 

SN-AZ 

2333.6 

237.9 

.N-NV 

2462.8 

256.6 

NP-NT 

-  -  -  - 

— 

37.2 
00.3 
04.6 
19.6 

51.5 

38.2 

43.3 

50.6 

53.2 
05.5 

15.4 
25.1 

32.3 

28.7 

38.8 

32.3 

02.2 

46.9 

02.2 

17.5 
20.7 

28.5 

30.3 


11:34:34.1 

11:34:37.4 

11:34:39.7 

11:34:40.8 

11:34:41.6 

11:34:47.0 

11:34:43.7 

11:34:46.9 

11:34:58.3 


TRAVEL 
TIME  (SEC) 


36.9 

60.0 


64.3 

79.3 

111.2 

151.9 
163.0 
170.3 

172.9 
185.2 


208.4 

218.5 

202.9 

241.9 

226.6 
241.9 

257.2 
260.4 

I 

268.2 
270.0 


273.8 

277.1 

279.4 

280.5 
2813 
281.7 
283.4 

286.6 
298.0 


AMPLITUDE 

mu 

27. 

5 

31. 

4 

10. 

4 

19. 

1 

92. 

7 

10. 

1 

(SEC) 


5.12 
4.64 
5.20 
6.19 
5.60 
5.58 
5.39 
5.08 
4.98 
4.72 
4.57 
4.45 
3.68 
4.09 
3.32 
2.78 
3.44 
2.88 

3.12 
2.85 
3.31 
3.15 


3.76 

3.46 
3.19 
3.65 
2.98 
3.63 
2.74 

3.47 
3.39 


A 


TAI5LK  A- 5. 

First  arrival  data  for  idiot  46A 


STATION 


GP-MN 


BR-PA 

BL-IW 

DH-NY 

CP0-08 

LS-NH 

WMO-06 

EU-AL 

HN-ME 

RT-NM 

JE-LA 

UBO-IO 

SV-2-QB 

DR -CO 

HL-2-ID 

NL-AZ 

BMO-03 

LC-NM 

KN-UT 

WO-AZ 

HR-AZ 

EK-NV 

GE-AZ 

TFO-60 

JR-AZ 

LG-AZ 

SG-AZ 

SN-AZ 

..N-NV 

NP-NT 


DISTANCE 

(KM) 

AZIMUIH 

DEGREES 

25C.8 

280.6 

423.7 

204.3 

468.6 

328.6 

583.8 

195.8 

851.6 

280.2 

1228.1 

127.  7 

1282.0 

142.2 

1334.2 

!  109.3 

1355.8 

161.8 

1459.8 

096.8 

1562.8 

209.6 

1624.3 

172.2 

1698.9 

086.1 

1654.3 

229.9 

1734.2 

250.5 

1806.8 

239.2 

1932.3  . 

266.5  j 

2042.0 

238.7 

3104.4 

272.6 

2159.7 

225.9 

2178.8 

246.6 

2189.5 

238.3 

2215.6  ’ 

238.1 

2245.2  • 

255.7 

2272.3 

235.9 

2279.8 

238.2 

2286.5 

238.9 

2307.1 

244.2 

2407.1 

256.7 

...  1 

— 

ARRIVAL 
TIME  (Z) 


11:30:37.9 

11:31:00.4 

11:31:05.4 

11:31:20.0 

11:31:52.0 

11:32:37.6 

11:32:42.7 

11:32:50.4 

11:32:52.7 

11: 33:04.8 

11:33:15.4 

11:33:25.3 

11:33:32.0 

11:33:28.3 


11:33:38.7 


11:33:47.0 

11:34:02.8 

11:34:17.5 

11:34:21.4 

11:34:28.8 

11:34:30.6 

11:34:32.3 

11:34:34.3 

11:34:37.6 

11:34:39.9 

11:34:41.1 


11:34:42.5 
11: 34:44.0 


11:34:58.7 


TRAVEL  |  AMPLITUDE 
TIME  (SEC) 


PERIOD 

(SEC) 


37.7 
60.2 
65.2 

79.8 

111.8 

157.4 

162.5 

170.2 

172.5 

184.6 

195.2 

205.1 

211.8. 

208.1 


218.5 


226.8 

242.6 

257.3 
261.2 

225.9 

270.4 

272.1 

274.1 

277.4 

279.7 

280.9 


282.3 

283.8 


298.5 


43.8 
14.6 

5.3 
6.  5 

79.8 

10.2 

4.2 
2.6 
3.0 

3.2 

7.4 
10.3 

2.5 

2.8 


O'.  4 
0.5 
0.3 
0.4 
0.4 
0.5 
0.3 
0.4 
0.3(5 
j  .  5 
0.8 
0.4 
0.5 
0.6 


4.56 
4.35 
4.61 
6.30 
5.51 
5.24 
4.92 
5.03 
4.81 

4.57 
4.71 
4.00 
3.96 


3.29 


2.81 
2.88 
2.79 
2.55 
3.04 
3.01 
3.44 
3.62 
2  c  96 
2.68 
3.58 


3.63 

3.00 


3.40 


*1 

f,  .t.)  .ift 


IvV':’..  *  - 

■  1 


STATION 
C,T-MN~ 
WI’-MN 
RK-ON 
VO-  IQ 
.  RY-NI) 

br-pa 
BL-IW 
DII-NY 
CP0-08 
LS-NJI 
WMO-06 
EU-AL 
HN-ME 
RT-NM 
UBO-IO 
SV-2-QB 
DR -CO 
ML- 2- ID 
NL-AZ 
BMO-03 
LC-NM 
KN-UT 
WO-AZ 
HR-AZ 
EK-NV 
GE-AZ 
TFO-60 
JR-AZ  • 
LG-AZ 
SG-AZ 
SN-AZ 
MN-NV 
P-NT 


DISTANCE 

(KM) 

|  AZIMUTH 
|  DEGREES 

244.3 

281.1 

419.8 

203.5 

465.7 

329.3 

580. 7 

195.1 

844.  7 

280.  3 

1232.1 

127.4 

1284.8 

141.9 

1339.3 
1356.5 

1465.3 
1558.2 

1623.8 

1704.4 
1648.0 

1726.8 

1799.8 

1924.5 
2034.1 

2096.5 

2153. 7 

2171.4 

2208.7 

2237.5 

2265.5 
2272.9 

2279.6 

2299.8 

2462.3 


109.1 
161.  5 
096.6 
209.4 

171.9 
086.0 

229.9 


TABLE  A- 
First  arrival  data 

muth  arriv'Z 

REES _ Tlf.tr  (Z) 

1.1  14:00:37.2 

3.5  14:01:00.0 

9.3  14:01:05.1  , 

5.1  14:01:19.5 

0.3  14:01:51.4 

7.4  14:02:38.0 

1.9  14: 02:42.8 

?.l  14:02:50.4 

U  5  14:02:52.8 

>.6  14:03:04.6 

>.4  14:03:15.0 

.9  14:03:24.8 

.0  14:03:31.6 

.9  14:03:28.6 


6. 

for  shot  4  51) 


239.0 
266.  5 

238.6 

272.6 
225.  7 
266.5 

238.0 

255.7 

235.8 

238.1 

238.8 

244. 1 


14:03:46.5 
14:04:02.1 
14:04:17.2 
14: 04: 20.7 
14:04:27.6 
14:04:30.0 


14:04:33.7 

14:04:37.2 

14:04:39.4 

14:04:40.5 


14: 04: 41; 7 
14:04:43.8 


TRAVEL 
TIME  (SEC) 

37.0 

59.8 

64.9 
79.3 


(AMPLITUDE  I  PERIOD 


218.5  I  14:03:38.7  I  250.4 


256.7  I  14:04:58.2 


226.3 
241.9 
257.0 
260.5 

267.4 
269.8 

273.5 
277.0 

279.2 

280.3 

281.5 

283.6 


298.0 


mu 

24.6 

18.1 

7.3 

16.9 


(SEC) 

0.3 

0.3 

0.3 

0.5 


111.2 

36.1 

0.3 

157.8 

8.8 

0.4 

162.6 

7.2 

0.4 

170.2 

>— • 
O 

• 

00 

0.4 

172.6 

1.8 

0.4 

184.4 

4.9 

0.6 

194.8 

7.7 

0.8 

206.6 

9.2 

0.5 

211.4 

1,6(5) 

0.4 

208.4 

00 

• 

o 

0.4 

0.6  0.6 

0.3  0.7 

1.7(5)  0.8 


4.88 

4.48 

4.93 

6.08 

5.55 

5.36 

5.53 

4.75 

4.91 

4.59 
4.57 

3.92 

3.60 

3.02 

2.51 

2.88 

2.79 

2.55 

3.16 

3.18 

3.67 
3.04 

2.68 

3.34 

3.40- 

3.06 


3.69 


TABLE  A- 7. 

First  arrival  data  for  shot  47A 


STATION 

DISTANCE 

(KM) 

AZIMUTH 

DEGREES 

ARRIVAL 

TIME  (Z) 

TRAVEL 
TIME  (SEC) 

AMPLITUDE 

mu 

PERIOD 

(SEC) 

Mb 

GP-MN 

258.9 

|  280.6 

14:30: 39.5 

39.2 

10.0 

0.4 

— 

WF-MN 

426.0 

205.4 

14:31:00.6 

60.3 

13.7 

0.5 

4.54 

RK-ON 

473.8 

327.9 

14:31:06.1 

65.8 

2.0 

0.3 

3.92 

VO- 10 

584.9 

196.6 

14:31 : 20. 8 

80.5 

3.  2 

0.3 

3.43 

RY-.ND 

— 

— 

- 

- 

_  - 

BR-PA 

1221.1 

128.0 

14:32:36.0 

155.7 

1.8 

0.3 

4.95 

BL-1W 

1276.2 

142.5 

14:32:41.7 

161.4 

8,5 

0.4 

5.53 

DH-NY 

1326. 2 

109.4 

14:32:49.0 

tl  6  8 .  7 

v  #  4 

0.4 

5.23 

CPO-08 

1352.3 

162.2 

14:32:51.8 

171.5 

1.7 

0.4 

4.74 

LS-NH 

1451.6 

096.8 

14:33:03.3 

183.0 

2.2 

0.5 

4.65 

WMO-36 

1565.9 

210.0 

14:33:15.4 

195.1 

2.5 

0.5 

4.29  : 

EU-AL 

— 

— 

— 

- 

- 

-  - 

HN-ME 

1690.9 

086.1 

14:33:31.0 

210.7 

2.1 

0.3 

4.15 

RT-NM 

1659.8 

230.2 

14:33:29.6 

209.3 

0.5 

0.4 

3.43 

JE-LA 

— 

— 

-  -  - 

- 

• 

I 

UBO-IO 

1741.5 

250.7 

14:33:40.5 

220.2 

0.6(5) 

0.5 

3.02 

SV-2-QB 

— 

— 

— 

- 

- 

-  - 

DR -CO 

1813.1 

239.4 

14:33:47.8 

227.5 

0.2 

0.4 

2.51 

HL- 2- ID 

— 

— 

— 

m 

•.  - 

NL-AZ 

— 

— 

— 

- 

• 

•  - 

BMO-03 

2112.5 

272.7 

14:34:22.0 

261.7 

0.2 

0.5 

2.55 

LC-NM 

2166.8 

226.1 

14:34:29.0 

268.7 

0.7(5) 

0.7 

3.03 

KN-UT 

2185.7 

246.8 

14:34:31.5 

271.2 

0.5 

0.6 

2.95 

WO-AZ 

— 

— 

— 

- 

• 

•  - 

HR-AZ 

2221.8 

238.3 

14:34: 35.0 

274.7 

3.5 

0.6 

3.77 

EK-NV 

2252.8 

255.9 

14:34:38.8 

278.5 

0.4 

0.6 

2.82 

GE-AZ 

2278.3 

236.1 

14:34:40.5 

280.2 

0.3  1 

0.6 

2.67 

TFO-60 

2286.1 

-  238.4 

14:34:41.9 

281.6  .  | 

1.4 

0.6 

3.46 

JR-AZ 

— 

— 

-  i 

- 

- 

-  - 

LG-AZ 

2292.9 

239.1 

14:34:43.0 

282.7  i 

2.0  i 

0.7 

3.55 

SG-AZ 

2313.9 

244.4 

14:34:45.4 

285.1 

0.1  i 

0.6 

2.60 

SN-AZ 

2345.4 

238.3 

14:34:47.9 

287.6  j 

0.6 

0.7 

3.06 

.N-NV 

2477.7 

256.9 

14:34:59.8 

299.5 

0.4 

0.6 

2.98 

NP-NT 

— 

— 

-  1 

-  ; 

- 

Mean 


3.7 


STATION 


GP-MN 

wf-mn’ 

RK-ON 
VO-IO 
RY-N1) 
BR-PA 
BL-WV 
I  DH-NY 
CP0-08 
I  LS-NII 
I  WMO-06 
I  EU-AL 
HN-ME 
RT-NM 
|'JE-LA 
I  UBO-IO 
I  SV-2-QB 
I  DR- CO 
HL- 2- ID 
I  NL-AZ 
I  BMO-03 
I  LC-NM 
KN-UT 
WO-AZ 
HR-AZ 
EK-NV 
GE-AZ 
!  TF0-60 
JR-AZ 
LG-AZ 
SG-AZ 
SN-AZ  | 
.N-NV 
•NP-NT 


DISTANCE 

(KM) 


44.2 

16.9 

67.8 


577.8 


1231.0 

1283.0 

1339.1 

1354.1 
1456. 7 
1555.4 

1621.1 


|  1705.3 

1725.2 

1797.8 

1923.6 
2033.0 

2095.9 

2151.2 

2169.6 

2180.5 

2206.6 

2236.1 

2263.4 
2270.8  . 

2277.5 
2298,0 

2330.1 
2461.0 


TABLE  A- 8. 

I 

First  arrival  data  for  shot  45C 


AZIMUTH 

DEGREES 


281.8 

203.5 

329.5 
195.1 


157.8 

141.7 
109.0 

161.4 
096.5 

209.4 

171.8 

085.9 

250.5 

239.1 
266.  5 

238.6 

272.6 

225.8 

246.6 
238.3 

238.1 

255.7 

235.9 

238.1 
_ 

238.9 

244.2 
238.0 

256.7 


ARRIVAL 
TIME  (Z) 


TRAVEL 
TIME  (SEC) 


37.1 
59.4 
65.  2 
78.9 


157.8 

162.6 


204.0 

211.9 

218.3 

226.1 

241.8 

256.3 
260.6 

267.7 

269.4 

271.4 

273.1 

277.1 
279.0 

280.2  . 

281.5 
283.0 
286.3 

297.8 


AMPLITUDE  I  PERIOD 
|  mu 


15.9 


4.41 

4.34 

4.72 


5.39 
5.13 

5.40 
4.74 
4.83 
4.57 
4.24 


4.04 


3.02 

2.51 
2.58 
2.49 
2.55 
3.31 
2.95 
3.29 
3.67 
3.12 
.  3.03 
3.27 


3.45 
2.42 

3.46 
3.51 


TABLE  A- 9. 

First  arrival  data  for  shot  73A 


STATION 

DISTANCE 

(KM) 

AZIMUTH 

DEGREES 

ARRIVAL 
TIME  (Z) 

1  travel 

1  TIME  (SEC) 

[AMPLITUDE  1  PERIOD 
rnu  (SEC) 

i  VI 

GP-MN 

» 

264.1 

2  78.4 

11:30:40.0 

39.7 

34.9 

0.5 

-  1 

WF-MN 

437.2 

205.8 

11:31:01.9 

61.6 

9.2 

0.5 

4.37 

RK-ON 

469.4 

326.7 

11:31:05.6 

65.3 

2.4 

0.4 

3.87 

VO-IO 

595.6 

197.1 

11:31:21.2 

80.9 

3.4 

0.4 

4.  33 

|| 

864.5 

279.6 

11:31:54.0 

113.7 

71.9 

0.5 

6. 16 

BR-PA 

1220.7 

128.6 

11:32:36.7 

156.4 

1.2(5 

;  0.5 

4.  70 

BL-WV 

1278.8 

143.1 

1  11:32:42.5 

162.2 

6.8 

0.4 

5.26 

1322.0 

110.0 

11:32:48.1 

167.8 

13.3 

0.4 

5.62 

CP0-08 

1358.6 

'  162.7 

11:32:53.3 

173.1 

3.9 

0.5 

5.00 

LS-mi 

097.3 

11:33:02.5 

182.2 

1.9 

0.4 

4.68 

WMO-06 

1576.9 

210.1 

11:33:17.1 

196.8 

1.3 

0.5 

4.02 

EU-AL 

f 

— 

HN-ME 

1682.1 

086.4 

11:33:29.8 

209.5 

2.9 

0.6 

3.98 

RT-NM 

1670.5 

230.2 

11:33:31.3 

211.0 

1.2 

0.5 

3.66 

JE-LA 

— 

— 

... 

. 

UBO-10 

1750.3 

250.6 

11:33:40.6 

220.3 

0.7 

0.5 

3.03 

SV-2-QB 

— 

— 

•»  « 

L-El 

— 

— 

HL-2-ID 

•  1946.7 

266.5 

11:34:04.6 

244.3 

0.3 

0.7 

2.58 

NL-AZ 

2058.4  • 

238.9 

11:34:19.2 

258.9 

1.0 

0.5 

3.31 

BMO-03 

2117.7 

272.5 

11:34:22.8 

262.5 

0.2 

0.5 

2.55 

LC-NM 

KN-UT 

2175.7 

226.2 

11:34:30.6 

'  270.3 

1.4 

0.5 

3.43 

WO-AZ 

— 

— 

HR-AZ 

EK-NV 

2231.9 

238.3 

11:34:36.0 

275.7 

2.4 

0.5 

3.68 

GE-AZ 

— 

— 

TF0-6& 

2296..  2 

238.4 

11:34:43.2 

282.9 

0.6 

0.5 

3.18 

JR-AZ 

-  ! 

— 

LG-AZ 

2302.9 

239.1 

11:34:43.8 

283.5 

2.0 

0.7 

3.55 

SG-AZ 

*  —  • 

— 

_ 

SN-AZ 

— 

N-NV 

NP-NT 

— 

- 

Mean 


4.0 


TABLE  A- 10. 

First  arrival  data  for  shot  45P 


TABLE  A- 11 


First  arrival  data  for  shot  48A 


* 


V 


STATION 

DISTANCE 

IKM) 

AZIMUTH 

DEGREES 

ARRIVAL 

TIME  (Z) 

TRAVEL 
TIME  (SEC) 

AMPLITUDE 

mu 

PERIOD 

(SEC) 

_ _Mb_ 

GP-MN 

264.0 

279.6 

14:00:39.9 

39.5 

11.8 

0.4 

WF-MN 

432.2 

206.0 

14:01:01.4 

61. 0 

5.  4 

0.5 

4.14 

RK-ON 

473.2 

327.1 

14:01:06.6 

66.2 

2.2 

0.3 

3.97 

■VC-10 

590.4 

197.1 

14:01:21.0 

80.6 

5.9  ! 

0.4 

4.57 

RY-ND 

864.5 

2  79.9 

14:01:54.1 

113.7 

33.  3 

0.4 

5.92 

BR-PA 

JPP 

128.4 

14:02:36.5 

156.1 

1.9 

0.4 

4, 8  7 

BL-WV 

142.9 

14:02:41.4 

161.0 

3.3 

0.5 

5.02 

109.7 

14:02:48. 5 

168.1 

11.8 

0.4 

5.57 

CP0-08 

162.5 

14:02:52.6 

172. 2 

1.2 

0.4 

4.57 

LS-NH 

1445.1 

097.1 

14:03:03.0 

182.6 

2.2 

0.5 

4.65 

WMO-06 

1572.0 

210.2 

14:03:16.6 

196.2 

1.2 

0.4 

4.10 

EU-AL 

— 

— 

— 

- 

- 

HN-ME 

1683.4 

086.3 

14:03:29.9 

209.5 

2.2 

0.4 

4.04 

RT-NM 

1666. 5 

230.  3 

14:03:30.8 

210.4 

0.5 

0.5 

3.09 

UBO-IO 

1747.7 

250.7 

14:03:40.9 

720.5 

0.  2 

0.4 

2.60 

SV-2-QB 

— 

— 

— 

- 

- 

DR- CO 

— 

— 

—  —  • 

- 

- 

HL- 2- ID 

.1945.5 

266.6 

14:04:04.5 

244.1 

0.1 

0.6 

2.26 

NL-AZ 

2054,9 

239.0 

14:04:19.6 

258.6 

0.4(5) 

0.6 

2.78 

BMO-03 

2117.1 

272.6 

14:04:22.5 

262.1 

0.2 

0.5 

2.55 

LC-NM 

— 

— 

.  - 

- 

- 

KN-UT 

2192.1 

246.9 

14:04:32.4 

272.0 

0.2 

0.4 

2.71 

WO-AZ 

— 

— 

*  - 

- 

HR-AZ 

2228.5 

238.4 

14:04:35.8 

275.4 

0.6 

0.4 

3.19 

EK-NV 

2258.6 

255.9 

14:04:39.5 

279.1 

0.2 

0.5 

2.65 

GE-AZ 

— 

— 

— 

- 

- 

TFO-60 

2292.7 

238.5 

14:04:42.7 

282.3 

0.6 

0.5 

3.17 

JR-AZ  • 

— 

— 

— 

- 

-  - 

LG-AZ 

2299.5 

239.2 

14:04:43.6 

283.2 

0.3 

0.4 

2.99 

,  SG-AZ 

2320 . 3 

244.4 

14:04:45.5 

285.1 

0.2 

0.6 

2.70 

SN-AZ 

2352.0 

238.4 

14:04:48.2 

287.8 

0.4 

0.6 

2.92 

MN-NV 

2483.5 

256.9 

14:05:00.3 

299.9 

0.4(5) 

0.7 

3.01 

NP-NT 

•  *  “ 

•  •  “ 

• 

Mean 

3.7 

TABLE  A- 12. 

First  arrival  data  for  shot  4  5A 


STATION 

DISTANCE 

(KM) 

AZIMUTH 

DEGREES 

ARRIVAL  ’ 
TIME  (Z; 

TRAVEL 
TIME  (SEC) 

A 

MPLITUDE 

mu 

PERIOD 

(SEC) 

Mb 

Gl'-MN 

243.9 

281.7 

12:00:37.4 

37.0 

3.5 

0.4 

W1;-MN 

417.1 

203.  5 

12:00:59.7 

59.3 

6.2 

0.5 

4.19 

RK-ON 

467.4 

329.5 

12:01:05.4 

65.0 

2.8 

0.3 

4.06 

VO-IO 

578„  1 

195.1 

12:01:19.4 

79.0 

5.0 

0.  4 

4.50 

RY-ND 

844.  2 

280.4 

12:01:52.3 

111.9 

12.9 

0.4 

5.  51 

BR-PA 

1231.4 

127.3 

12:02:37.9 

157.5 

3.8 

0.4 

5.17 

BL-WV 

1283.4 

141.7 

12:02:42.7 

162.3 

,  3.3 

0.5 

4.92 

1339.4 

109.0 

12:02:50.5 

170.1 

6.0 

5.18 

C?0-08 

1354.4 

161.4 

12:02:53. 1 

172.7 

!  !•» 

|  4.75 

LS-M1 

1465.9 

096.  5 

12:03:05.7 

185.3 

i  2.2  j 

4.65 

WMO-06 

1555.5 

709.4 

12:03:14.7 

194.3 

0.8 

^9 

4.22 

EU-AL 

— 

— 

1 

- 

-  - 

HN-ME 

1705.5 

C85.9 

12:02:31.5 

211.2  ( 

2.2  | 

|| 

3*  86 

RT-NM 

1645.7 

229.8 

12:03:28.8 

208.4 

0.5 

3.42 

UBO-IO 

1725.1 

250.5 

12:03:38.5 

218.2 

0.4  ; 

7.90 

SV-2-QB 

— 

-  -  - 

-  ! 

...  1 

- 

- 

•  • 

1797. 7  j 

239.1 

12:03:46.4  ! 

0.2 

0.4 

2.52 

ML-  2-ID 

- - 

— 

_ 

•> 

- 

-- 

NL-AZ 

'  2032.9 

238.6 

12:04:16.5 

256.1 

0.2 

2.54 

BMO  03 

2095.7 

272.6 

12:04:21.4 

261.0 

0.2 

2.55  ; 

J.C-NM 

2151.3 

225.8 

12:04:27.8 

267.4 

0.7 

3.16 

KN-UT 

2169.5 

246.6 

12:04:30.1 

269.7 

0.3 

2.65 

WO-AZ 

-  -  - 

— 

— 

- 

- 

-- 

HR-AZ 

2206.5 

238.1 

12:04:33. 5 

273.1 

1.6 

3.43 

EK-NV 

2235.9 

255.7 

12:04:37.0 

276.6 

0.2(5) 

2.63 

GE-AZ 

2263.3 

235.9 

12:04: 39.2 

278.8 

0.1 

2.39 

TFO-60 

2270.8 

238.1 

12:04:40.5 

280.1 

0.6 

3.03 

JR-AZ 

«•  -  - 

— 

- 

- 

-  - 

LG-AZ 

3 ] 77. 5 

238.9 

12:04:41.2 

0.6 

3.19 

SG-AZ 

2297.9 

244.4 

12:04:43.4 

283.0 

0.4 

0.5 

SN-AZ 

2330.  ; 

238.0 

12:04:46.6 

286.2 

1 

0.7 

Q I 

3.  74 

MN-NV 

2460.8 

256.7 

12:04:58.3 

297.9 

0.4 

C9 

2.99 

T-NT 

•  -  - 

— 

— 

- 

- 

-- 

| 

Mean 

3.6 
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